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Laboratory  research  activities  (under  Grant  AFOSR  91-0170)  for  the  one-year  period  ending  14 

February  1992.  Among  our  research  results  described  in  detail  in  the  cited  references  (Section  5). 

Perhaps  the  most  notweworthy  are  the  development  of: 

R 1  rational  correction  factors  to  account  for  the  effects  of  suspended  particle  morphology  on  convective- 

diffusion  mass  deposition  rates 

R2  quantitative  criteria  for  influence  of  particle  ihermophoresis  on  the  structure  of  two-phase 

laminar  counterflow  diffusion  flames  (potentially  useful  to  predict  IR  radiation  from  such 
flames  and  optimize  particle  properties  in  synthesis  applications) 

R3  quantitative  methods  for  predicting/coirelating  the  effects  of  particle  inertia  on  theimophoretic  deposition 

across  laminar  boundary  layers  on  targets  with  streamwise  curvature  (experimentally  verified  by  our  seeded 
micro-combustor  experiments  on  concave  ribbon  targets) 

15  presentations  and  2  PhD  dissertations  have  resulted  from  this  research  program.  Copies  of  3 

reprints  appearing  during  this  period  are  included  in  the  Appendices  (Section  6)  of  this  report. 
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TRANSPORT  PHENOMENA  AND  INTERFACIAI.  KINETICS 
IN  MULTIPHASE  COMBUSTION  SYSTEMS 


1.  INTRODUCTION 

The  performance  of  ramjets  burning  slurry  fuels  Grading  to  condensed  oxide  aerosols  and 
liquid  film  deposits),  gas  turbine  engines  in  dusty  atmospheres,  or  when  using  fuels  from  non- 
traditional  sources  ,  depends  upon  the  formation  and  transport  of  small  particles  across  non- 
isothermal  combustion  gas  boundary  layers  (BLs).  Even  airbreathing  engines  burning  "clean" 
hydrocarbon  fuels  can  experience  soot  formation/deposition  problems  (e.g.,  combustor  liner 
burnout,  accelerated  turbine  blade  erosion  and  "hot"  corrosion).  Moreover,  particle  formation  and 
transport  are  important  in  many  chemical  reactors  used  to  synthesize  or  process  aerospace  materials 
(turbine  blade  coatings,  optical  waveguides,  ceramic  precursor  powders,...).  Accordingly,  our 
research  is  directed  toward  providing  chemical  propulsion  systems  engineers  and  materials- 
oriented  engineers  with  new  techniques  and  quantitative  information  on  important  particle-  and 
vapor-mass  transport  mechanisms  and  rates. 

The  purpose  of  this  report  is  to  summarize  our  research  methods  and  accomplishments 
under  AFOSR  Grant  91-0170  (Technical  Monitor:  J.M.Tishkoff)  during  the  l-year  period:  15 
February  '91-14  February  '92.  Readers  interested  in  greater  detail  than  contained  in  Section  2  are 
advised  to  consult  the  published  papers  explicitly  cit^  in  Sections  2  and  5.  Copies  of  any  of  these 
published  papers  (Section  5.2)  or  preprints  (Section  5.3)  can  be  obtained  by  writing  to  the  PI: 
Prof.  Daniel  E.  Rosner,  at  the  Department  of  Chemical  Engineering,  Yale  University,  Box  2159 
Yale  Station,  New  Haven  CT  06520-2159  USA.  Comments  on,  or  examples  of,  the  applications 
of  our  research  (Section  3.4.)  will  be  especially  welcome. 

An  interactive  experimental/theoretici  approach  has  been  used  to  gain  understanding  of 
performance-limiting  chemical-,  and  mass/energy  transfer-phenomena  at  or  near  interfaces.  This 
included  the  development  and  exploitation  of  seeded  laboratory  flat  flame  burners  (Section  2.1), 
and  new  optical  diagnostic/spectroscopic  techniques.  Resulting  experimental  rate  data,  together 
with  the  predictions  of  asymptotic  theories  (Section  2),  were  used  as  the  basis  for  proposing  and 
verifying  simple  viewpoints  and  effective  engineering  correlations  for  future  design/optimization 
studies. 

2  RESEARCH  ACCOMPLISHMENTS  AND  PUBLICATIONS 

Most  of  the  results  we  have  obtained  under  Grant  AFOSR  91-0170  can  be  divided  into  the 
subsections  below: 

2. 1 .  TRANSPORT  OF  AGGREGATED  PARTICLES 

The  Brownian  diffusion-,  inertial-,  and  optical-properties  of  aggregated  particles,  as  formed  in 
sooting  diffusion  flames,  are  quite  sensitive  to  size  (eg.  number  N  of  "primary"  particles)  and 
morphology  (configuration  of  the  N  particles).  In  Fig.  1  (eg.)  we  summarize  the  orientation- 
averaged  Brownian  diffusivity  of  aggregates  in  the  continuum  (high  presure)  limit.  Of  interest  are 
meth^s  for  anticipating  coagulation  and  deposition  rates  of  populations  of  such  particles, 
especially  in  non-isothermal  flow  systems.  We  predict  that  capture  rates  from  "coagulation-aged" 
polydispersed  aerosols  of  a  particular  moiphology  by  the  mechanisms  of  convective-diffusion  and, 
especially,  thermophoresis  will  be  close  to  those  calculated  if  all  such  particles  had  the  mean 

particle  volume  <|)p/Np  ,where  4)p  is  the  particle  volume  fraction  and  Np  the  aggregate  number 
density.  Correction  factors  (Fig.2)  for  convective-diffusion  are  found  to  be  ca.  0.9  for  a 
surprisingly  wide  range  of  gas-dynamic/environmental  conditions,  and  particle  morphologies. 
However,  clearly,  if  one  does  not  know  the  suspended  particle  morphology  (distribution),  a 
knowledge  of  volume-  (or  mass-)  fraction  alone  is  not  enough  to  accurately  preset  convective- 
diffusion  deposition  rates,  with  (Section  2.2)  or  without  simultaneous  inertial  effects. 
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The  ability  to  reliably  measure  and  predict  the  thermophoretic  properties  of  aggregated 
flame-generated  particles  (carbonaceous  soot,  AI2O3,  ...)  is  important  to  many  technologies, 
including  chemical  propulsion  and  refractory  materials  fabrication.  We  are  now  completing  a 
manuscript  describing  our  measurements  of  the  size-  and  morphology  insensitive  thermophoretic 
diffusivity  of  flame-generated  submicron  TiOafs)  "soot"  particles  using  TiCl4(g)-seeded  low 
strain-rate  counterflow  laminar  diffusion  flame  (CDF-)  techniques  (Fig.  3,  Gomez  and  Rosner, 
1992)).  One  of  our  important  conclusions  is  that  engineering  c^culations  of  thermophoretically 
dominated  soot  deposition  rates  will  usually  not  require  details  of  the  aggregate  size-  and 
morphology-  distribution.  Depending  on  the  sharpness  of  the  gas  temperature  and  velocity 
gradients  in  counterflow  diffusion  flames  there  are  one  or  more  particle  stagnation  pIanes.(PSPs, 
Fig.4).  A  knowledge  of  the  relative  positions  of  the  gas  and  particle  stagnation  planes  (shown 
(Fig.5)  vs.  flame  stoichiometry),  and  the  associated  chemical  environments,  can  be  used  to  control 
the  composition  and  morphology  of  flame-synthesized  particles.  These  factors  should  also 
influence  particle  production  and  radiation  from  turbulent  non-premixed  "sooting"  flames,  as  will 
be  discussed  further  in  Gomez  and  Rosner,  1992 


2.2.  MULTIPHASE  BOUNDARY  LAYER  THEORY:  THERMOPHORESIS  AND  INERTIA 

The  results  of  Seeded  micro-combustor  experiments,  and  ancillary  theoretical  calculations 
on  the  interesting  competition  between  particle  inertia  and  particle  thermophoresis  for  the  case  of 
laminar  gaseous  boundary  layers  on  surfaces  with  streamwise  curvature  {eg.,  turbine  blades)  are 
shown  in  Fig.6  (from  Konstandopoulos  and  Rosner,  1992;  Rosner  et.  al.  1991).  The  combined 
effects  of  thermophoretic-  and  inertial-particle  drift  (here  for  compact  aggregates  of  MgO(s)  with 

volume  mean  diameter  of  ca.  0.8  |J.m)  toward  a  concave,  "cold"  target  are  seen  to  substantially 
augment  local  particle  deposition  rates,  in  excellent  accord  with  the  predictions  of  two-phase 
laminar  bound^  layer  theory.  As  implied  by  Fig.l,  morphology  will  influerre  inertially  modified 
deposition  rates,  reported  here  in  terms  of  the  dimesionless  Stanton  number  (Rosner,1986). 
Indeed,  if  these  MgO(s)  particles  were  "open"  agglomerates  of  equal  mass  the  relevant 
dimensionless  stopping  time  parameter:  Stk*Re^/2  (where  Stk  is  the  ordinary  Stokes  number 
(Rosner, 1986)  and  Re  the  prevailing  length-Reynolds  number  evaluated  at  the  deposition  location) 
would  have  been  too  small  to  have  a  measureable  effect  for  the  cases  shown.  On  convex  cooled 
surfaces  {eg.  the  "suction"  side  of  a  GT  blade)  our  prediction-correlation  methods  predict  the 
competing  effects  of  thermophoretically-augmented  capture  and  inertial  drift  away  from  the  target. 
Instructive  results  for  the  angular  distribution  of  mass  transfer  to  cooled  cylinders  and  spheres  are 
displayed  in  Fig.  7.  Our  recommended  correlation/prediction  techniques  for  high  Schmidt  number 
combined  inertia/phoretic  deposition  will  be  the  subject  of  a  follow-on  paper. 

2.3.  HEAT  TRANSFER  EFFECTS  ON  COAGULATION  DYNAMICS 

Our  continuing  studies  of  submicron  particle  motion  in  host-gas  temperature  gradients  and 
radiation  fields  (see,  eg.,  Mackowski  et.al.\992,  Rosner,  et.  a/.,1992,  Mackowski,  D.W.,1990) 
have  demonstrated  that  even  spherical  particles  thermally  out  of  equilibrium  with  their  local  host 
gas  will  coagulate  with  particle-particle  encounter  rate  constants  quite  different  from  their  usual 
"isothermal"  Brownian  values  (see,  eg.,  the  "correction"  factors  of  Fig.8  for  unequal-sized 
carbonaceous  particles  in  a  black-body  radiation  environment  with  the  stated  temperatures).  This 
leads  to  unusual  population  dynamics,  including  the  strong  tendency  of  radiatively  cooled  large 
particles  to  grow  still  larger  by  rapidly  "scavenging"  smaller  ones.  Localized  radiative  cooling  can 
also  thermophoretically  concentrate  the  radiating  particles,  with  interesting  stability  consequences 
(Mackowski  et.al.  1992 

2.4.  KINETICS  AND  MORPHOLCXJY  OF  CVD-MATERIALS  IN  MULTI-PHASE  ENVIRONMENTS 

A  small  impinging  jet  (stagnation  flow)  reactor  is  being  used  to  study  the  CVD-rates  of 
refractory  layers  on  inductively  heated  substrates,  including  those  deliberately  "pre-loaded"  with 
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a  TiKnnophortiically-iIaininaicd  deposition  rate  ratios  (closer  to 
unity)  not  reported  here;  present  values  are  for  negligible 
transport  by  non-convcctive 

mechanisms  other  than  Brownian  diffusion  (numerical  columns 
1-4)  or  inertia  (column  S) 

b  log-normal  orientation  averaged  populations  with  Og(fm)»2.46 
and  ag(c)  >2.30  (sec.r.x.  Rosner  it  Tassopoulos(1989;) 
d  Reference  ease;  sintered  single  spheres  with  abovementioned  Og- 
values  for  BSD 

c  For  continuum  cases  not  strictly  a  power-law  (Fig  1), 
however,  values  stated  ate  for  approximate  power- 
law  representation  for  <Dn(N» 


Fig.l.  Predicted  nonnalized  orientation-averaged  Brownian 
diffusivity  of  aggregated  particles  containing  N-primary 
spheres  in  the  near-continuum  limit;  Normalization: 
Brownian  diffusion  coefficient,  D],  of  primary  sphere  in  the 
same  local  environment  (after  Rosner  etMl.  1991) 


Fig.  2  Predicted  values  of  (•m"p/(-m"p(v)]  for  isothermal 
deposition  from  coagulation-aged  aggregated  particle 
peculations  (after  Rosner  etMl.,1992) 
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Fig.3  a).  Axisymmetric  counterflow  diffusion  flame  (CDF) 
burner  configuration  used  for  experimental  determination  of 
panicle  thcrmophoretic  diffusivity  based  on  observed  axial 
position  of  "particle  stagnation  planes"  (PSPs);  b) 
Appearance  of  axially  mounted  125  pm  diam.  filament 
revealing  existence  of  dust-free  region  straddling  the  flat 
diffusion  flame  (after  Gome/  and  Rosner,1992) 


Fjg.4  Axial  gas  velocity  and  particle  thermophoretic  velocity 
profiles  in  near-stoichiometric  counterflow  laminar  diffusion 
flames;  a)  "Sharp  flame-modest  strain"  case  with  three 
eligible  panicle  stagnation  planes  (PSPs);  b)  Broad  flame- 
high  strain"  case  with  one  eligible  particle  stagnation  plane 
(after  Gomez  and  Rosner,1992) 
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Fig.5  Locus  of  dimensionless  axial  positions  of  particle 
stagnation  plane(s)  (PSP)  and  gas  stagnation  plane  (GSP) 
(relative  to  the  flWe  location,  zf)  for  laminar  counterflow, 
diffusion  flames  with  fuel  velocity  Up  as  a  function  of  the 
effective  fuel  vapor  equivalence  ratio',  ^^ggr.Case  a)  "Sharp 
flame-modest  strain  rate  a”.  Case  b  "Broad  flame-high  strain 
rate  a”  (after  Gomez  and  Rosner,  1992) 


Inertiol  Effects  on  Thermophoretic  Deposition 
to  Wolls  with  Streomwise  Curvature 
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Fig.6  Experimental  and  theoretical  dimensionless  particle 
deposition  rate  coefficient  (Stanton  number)  StmRcx*^  vs. 
surface-to-gas  temperature  ratio:  Tw/Te  and  the  relevant 
particle  relaxation  time  parameter  (effective  Stokes  number ) 
Stl^ex^^ ;  laminar  boundary  layer  flow,  concave  circular  arc 
solid  target  (afto-  Konstandopoulos  and  Rosner,  1992) 


F^.  7  Predicted  inertial  effects  on  thermophoretic-convective 

mass  transfer  rates  to  cooled  convex  objects  at  subcritical  Fig.  8.  Multiplicative  "correction"  to  the  Brownian 
Stokes  numbers;  angular  distribution  of  mass  transfer  (Smoluchowski  isothermal)  coagulation  rate  constant  due  to 
coefficient  for  spheres  and  cylinders  at  high  Reynolds  energy  exchange  to/from  spherical  particles  of  unequal  radius 
numbers  (Normalization:  Stagnation  point  value  at  the  (after  Rosner  cr.a/.,1992,  Mackowski  ef.n/.  1992) 
prevailing  value  of  Stk)  (after  Konstandqx>ulos(1991)) 


porous  particulate  deposits.  These  measurements  and  associated  "chemical  vapor  infiltration" 
theory  will  be  useful  to  understand  deposition  rates  and  deposit  microstructuies  that  have  been 
observed  in  particle-containing  CVD  environments.  Analogous  questions  arise  in  anticipating  the 
vapor  scavenging-  (or  sorption-)  ability  of  coagulation-aged  populations  of  suspended  aggregated 
particles,  on  which  we  are  also  working. 

3|c3|c4c4t 

In  our  OSR-sponsored  Yale  HTCRE  Lab  research  during  the  past  year,  briefly  reviewed 
here,  we  have  shown  that  new  methods  for  rapidly  measuring  particle-mass  transfer  rates 
combined  with  recent  advances  in  transport  theory,  provide  useful  means  to  identify  and 
incorporate  important,  but  previously  neglected,  mass  transport  phenomena  in  many  propulsion 
engineering  and  materials  engineering  design/optimization  calculations. 

Despite  the  formidable  complexities  to  be  overcome  in  the  design  and  operation  of  air- 
breathing  propulsion  power  plants  utilizing  a  broad  spectrum  energetic  fuels  tiiese  particular 
techniques  and  results  are  incticative  of  the  potentially  useful  simplifications  and  generalizations 
which  have  emerged  from  our  present  fundamental  AFOSR-funded  research  studies  of 
combustion-generated  particle  transport  mechanisms.  It  is  hoped  that  this  presentation  and  its 
supporting  (cited)  papers  will  facilitate  the  retinement  and/or  incorporation  of  some  of  the  present 
ideas  into  engineering  design  procedures  of  much  greater  generality  and  reliability.  This  work  has 
alredy  helped  identify  new  directions  where  research  results  would  have  a  significant  impact  on 
engineering  practice. 


3.  ADMINISTRATIVE  INFORMATION: 

PERSONNEL,  PRESENTATIONS,  APPLICATIONS, 

"COUPLING'  ACTIVITIES 

The  following  sections  summarize  some  pertinent  'non-technical'  facets  of  the 
abovementioned  Yale  HTCRE  Lab/AFOSR  research  program: 

3.1  Personnel 

The  present  results  (Sections  2,4  and  5)  are  due  to  the  contributions  of  the  personnel  listed 
in  Table  3.1-1,  which  also  indicates  the  role  of  each  researcher  and  the  relevant  time  interval  of  the 
activity.  It  will  be  noted  that,  in  addition  to  the  results  themselves,  this  program  has  simultaneously 
contributed  to  the  research  training  of  a  number  of  students  and  fresh  PhDs,  who  will  now  be  in  an 
excellent  position  to  make  future  contributions  to  technologies  oriented  toward  air-breathing 
chemical  propulsion,  and  high-tech  materials  processing. 

3.2  Cooperation  with  US  Industry 

The  research  summarized  here  was  supported  by  AFOSR  under  Grant  91-0170  (2/15/91- 
2/14/92).  The  Yale  HTCRE  Laboratory  has  also  been  the  beneficiary  of  continuing  smaller  grants 
from  U.S.  industrial  corporations,  including  GE-Schenectady,  DuPont,  and  Shell  as  well  as  the 
feedback  and  advice  of  principal  scientists/engineers  from  each  of  these  corporations  and 
Combustion  Engineering-ABB.  We  appreciate  this  level  of  collaboration,  and  expect  that  it  will 
accelerate  inevitable  applications  of  our  results  in  areas  relevant  to  their  technological  objectives 
(see,  also.  Section  3.4,  below). 

3.3  Presentations  and  Research  Training 

Apart  from  the  publications  itemized  in  Section  4  and  our  verbal  presentations  (of  progress) 
at  regular  AFOSR  Contractors  Meetings,  our  results  have  also  been  presented  at  annual  or  topical 
conferences  of  the  following  professional  organizations: 

AAAR  Combustion  Inst.(Central  States  Section) 

AIChE  ASME-Engrg.  Fdm. 

MRS 

Indeed,  our  paper  on  the  subject  matter  of  Section  2.3  was  awarded  first  prize  at  the  April  '91 
annual  meeting  of  the  Central  States  Section  of  the  Combustion  Inst. 


Table  3.1-1  Summary  of  Research  Participants^  on  AFOSR  Grant ; 

TRANSPORT  PHENOMENA  AND  INTERFACIAL  KINETICS 
IN  MULTIPHASE  COMBUSTION  SYSTEMS 


Name 

Princioal  Research  Activitvb 

Castillo,  J. 

vs 

3/15-4/15 

interface  stability 

Collins,!. 

GRA 

'91, '92 

CVD  of  ceramic  coatings 

Garcia- Ybarra,  P. 

VS 

3/15-4/15 

thermophoresis  (kin.  theory) 

Gomez,A. 

I-Asst.Prof. 

'91, '92 

meas.  of  particle  transp.  props.(CDFs) 

Konstandopoulos,A-j-. 

GRA 

'91 

thermophoresis/inertia  coupling 

Kho,T. 

GRA 

'92 

chemical  vapor  infiltration  (coatings) 

Labowsky,M. 

VS 

'91,'92 

method  of  images  calculations 

PapadopoulosJD 

GRA 

'92 

transport  phenomena  in  CVD  reactors 

Rosner,  D.E. 

PI 

'91-'92 

program  direction-dep.theory/exp 

Tassopoulos.M.f 

GRA 

'91 

aggregate  particle  deposition 

Tandon,  P. 

GRA 

'91-'92 

transport  phenomena  in  BLs  and  CDFs 

Tsang,J. 

UGRA 

Summer  '91 

CVD  of  Ti02(s)  Films 

a  PDRA=Post-doctoral  Research  Asst  (arrival  delayed  to  ’92)  GRA=  Graduate  Research  Assistant 
PI  =  Principal  Investigator...,VS  =  Visiting  Scholar 
b  See  Section  5  for  specific  references  cited  in  text  (Section  2);  +  PhD  work  completed  during  '91 


In  addition,  during  the  period:  2/15/91  -2/14/92,  the  PI  presented  seminars  at  the  following 


Universities: 

U.  Wisconsin 

Waterloo 

Trondheim-SINTEF 

U.  Oslo 

KTH-Stockholm 

U.Limoges 

U.Toulouse 

Technion  (Haifa) 

U.  Paris-Nord 

NTU-Athens 

In  all,  a  total  of  about  15  talks  were  given  by  Yale-HTCRE  people  based  in  part  on  the 
results  of  this  AFOSR  research  program. 

This  program  involved  the  dissertation  research  of  two  graduate  students 
(A.G.Konstandopoulos  and  M.  Tassopoulos;  cf.Table  3.1-1)  who  completed  their  PhD  degree 
requirements  in  1991. 

3.4  Known  Applications  of  Yale-HTCRE  Lab  Research  Results 

It  has  been  particularly  gratifying  to  see  direct  applications  of  some  of  this  AFOSR- 
supported  particle  and  vapor  mass  transfer  research  in  more  applications-oriented  investigations 
reported  in  recent  years.  Indeed,  the  writer  would  appreciate  it  if  further  examples  known  to  the 
reader  can  be  brought  to  his  attention. 

Explicit  examples  are  provided  in  ongoing  work  at  MIT,  and  Sandia  CRF,  both  groups 
having  incorporated  our  rational  correlation  of  inertial  particle  impaction  (e.g.  a  cylinder  in  cross- 
flow)  in  terms  of  an  effective  Stokes  number. 

This  PI  was  also  pleased  to  confirm  applications  of  our  AFOSR  and  DOE-supported 
resepch  (on  the  correlation  of  inertial  impaction  by  cylinders  in  crossflow)  by  the  National 
Engineering  Laboratory  (NEL)  of  Glasgow  Scotland  (Contact:  Dr.  Andrew  Jenkins).  NEL  is 
apparently  collaborating  with  Marchwood  Labs-CEGB  on  developing  mass-transfer  prediction 
methods  applicable  to  waste-heat  recovery  systems  in  incinerators,  as  well  as  pulverized  coal-fired 
boilers.  These  applications  are  somewhat  similar  to  those  reported  by  the  Combustion  Lab  R&D 
group  at  MIT 

In  the  area  of  alkali  sulfate  vapor  deposition  in  combustion  systems  additional  applications 
of  our  predictive  methods  (for  "chemically  frozen"  and  LTCE  multicomponent  laminar  boundary 


layers)  continue  to  be  made  by  British  Coal  Corporation-Power  Generation  Branch  (I.  Fantom, 
contact)  in  connection  with  their  topping  cycles  which  run  gas  turbines  on  the  products  of  fluidized 
bed  coal  combustors/gasifiers.  Our  Poster  in  the  1988  Combustion  Symposium  (Seattle)  appears  to 
have  motivated  Prof.  Takeno  (Nagoya)  to  initiate  experiments  using  a  laminar  counterflow 

diffusion  flame  to  estimate  the  thermophoretic  diffusivity  of  LDV  seed  particles.  Our  paper  on  the 
experimental  determination  of  smaller  aggregated  particle  thermophoretic  diffusivities  using 
CDF/LDVA^aser  Light  Scattering  techniques  is  nearing  completion,  and  will  be  submitted  to 
Combustion  Science  and  Technology  early  in  our  next  year’s  program.  Also,  in  combustion 
research  many  groups  (eg.  Dobbinbs  et.al.  (Brown  U.),  Faeth  et.al.  (U.  Mich.),  Katz  et  al. 
(J.Hopkins  U.))  are  now  utilizing  "thermophoretic  sampling"  techniques  to  exploit  the  size-  and 
morphology-insensitive  capture  efficiency  characteristics  that  we  have  proven  in  our  AFOSR 
research  (Section  2.1). 

Explicit  use  of  our  studies  of  self-regulated  "capture"  of  incident  impacting  particles 
(Rosner  and  Nagaragan,1987)  is  being  made  in  current  work  on  impact  separators  and  ceramic 
heat  exchangers  for  coal-fired  turbine  systems  in  high  performance  stationary  power  plants.  Other 
potential  applications  arise  in  connection  with  "candle  filters"  used  to  remove  fines  (sorbent 
particles,...)  upstream  of  the  turbines.  A  useful  summary  of  work  in  these  interrelated  areas  (Solar 
Turbines,  Textron  Defense  Systems,  Hague  International,...)  was  presented  at  the  Engineering 
Foundation  Conference  (Palm  Coast,  FL  for  March  1991):  Inorganic  Tran^ormations  and  Ash 
Deposition  During  Combustion. ,  to  appear  (ASME/EF)  in  1992. 

Clearly,  fruitful  opportunities  for  the  application  of  our  recent  "non-Brownian"  convective 
mass  transfer  research  now  exist  in  many  of  the  programs  currently  supported  by  the  US  Air 
Force. 


4.  CONCLUSIONS 

In  the  OSR-sponsored  Yale  HTCRE  Lab  research  during  the  period:  2/15/91-2/14/92, 
bnetiy  described  above,  we  have  shown  that  new  methods  for  rapidly  measuring  particle-mass 
transfer  rates  ,  combined  with  our  recent  advances  in  mass  transport  theory,  provide  useful  means 
to  identify  and  inco^orate  important,  but  previously  neglected,  mass  transport  phenomena  in  many 
propulsion  engineering  and  materials  engineering  design/optimization  calculations. 

Despite  formidable  complexities  to  be  overcome  in  the  design  and  operation  of  mobile  and 
stationary  power  plants  utilizing  a  broad  spectrum  energetic  fuels  ^e  abovementioned  techniques 
and  results  (Section  2)  are  indicative  of  the  potentially  useful  simplifications  and  generalizations 
which  have  emerged  from  our  present  fundamental  AFOSR-funded  research  studies  of 
combustion-generated  particle  transport  mechanisms  and  interfacial  reactions.  It  is  hoped  that  this 
report  and  its  supporting  (cited)  papers  will  facilitate  the  refinement  and/or  further  incorporation  of 
many  of  the  present  ideas  into  design  and  test  procedures  of  greater  generality  and  reliability.  This 
work  has  also  helped  identify  new  Sections  where  it  is  anticipated  that  research  results  will  have  a 
significant  impact  on  future  engineering  practice. 
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The  Mahiliiv  of  the  planar  interface  of  a  structureless  solid  growing  from  a  depositing  component  dilute  in  a  carrier  fluid  is 
studied  when  the  mam  solute  transport  mechanism  is  thermal  (Sorett  diffusion.  A  linear  stahilitv  analysis,  carried  out  in  the  limit  of 
low  growth  Peclet  number,  leads  to  a  dispersion  relation  which  shows  that  the  planar  front  is  unstable  either  when  the  thermal 
diffusion  lactor  of  the  condensing  component  is  positive  and  the  latent  heat  release  is  small  or  when  the  thermal  diffusion  factor  is 
negative  and  the  vilid  grows  over  a  thermally-insulating  substrate.  Furthermore,  the  influence  of  interfacial  energy  effects  and 
constitutional  supersaturation  in  the  vicinity  of  the  moving  interface  is  analyzed  in  the  limit  of  very  Schmidt  numbers  (small 
solute  Ficliian  diffusioni  The  analysis  is  relevant  to  physical  vapor  deposition  of  very  massive  species  on  cold  surfaces,  as  in  recent 
evpenments  of  organic  wilid  film  growth  under  microgravity  conditions. 


I.  Iniroduclion 

The  growth  of  solid  deposits  from  the  vapor 
phase  is  a  technique  widely  used  for  the  fabrica¬ 
tion  of  thin  layers  of  metals,  insulators,  and  semi¬ 
conductor  materials;  see.  for  instance,  the  review 
papers  by  Stringfellow  (1]  and  Bryant  (21.  Even 
though  crystal  growth  from  melts  may  be  a  more 
"efficient”  process  (usually  giving  growth  rates 
several  orders  of  magnitude  larger  than  vapor 
phase  growth),  the  latter  is  often  preferred  be¬ 
cause  of  film  quality  considerations.  In  particular, 
new  pharmacological  and  electro-optical  applica¬ 
tions  demand  the  production  of  crystalline  films 
of  organic  and  organometallic  macromolecules. 
The  optimal  growth  of  these  materials  from  the 
vapor  phase  is  now  under  intense  experimental 


investigation  internationally,  even  in  microgravity 
environments  to  avoid  crystal  imperfections  in¬ 
duced  by  buoyancy-driven  instabilities,  see  Wal¬ 
ter  (3]. 

Especially  since  the  pioneering  work  of  Mullins 
and  Sekerka  (4.5)  and  the  review  by  Langer  (6).  it 
is  well  known  that  planar  crystal  growth  may  be 
morphologically  unstable  during  the  solidification 
of  a  supercooled  liquid  (i.e..  in  the  absence  of 
capillarity,  a  forward-facing  bulge  in  the  interface 
steepens  the  thermal  gradient  in  the  fluid  ahead 
of  it.  so  that  latent  heat  is  more  rapidly  removed 
from  the  surface  and  the  bulge  amplifies  unsta¬ 
bly),  and  in  the  similar  case  of  isothermal  solidifi¬ 
cation  of  a  liquid  mixture  (respectively,  direc¬ 
tional  solidification  of  a  binary  liquid)  due  to  the 
constitutional  supersaturation  (respectively,  const i- 
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tutional  supercooling)  induced  ahead  of  the  ad¬ 
vancing  solidification  front.  In  the  above-men¬ 
tioned  cases,  the  morphological  instability  arises 
as  a  consequence  of  the  existence  of  a  metastable 
liquid  phase  near  the  fluid-solid  interface. 

A  simplified  criterion  for  the  existence  of  con¬ 
stitutional  supersaturation  was  given  by  Reed  and 
LaFleur  [7]  and  reconsidered  by  Rosenberger  ct 
al.  (Sj.  This  criterion  leads  to  the  existence  of  a 
critical  temperature  gradient  reached  when  the 
concentration  gradient  at  the  interface  just  over¬ 
comes  the  gradient  of  the  saturation  concentra¬ 
tion.  Temperature  and  concentration  gradients 
increase  with  the  growth  velocity;  then,  for  a 
large  enough  crystal  growth  velocity  constitu¬ 
tional  supersaturation  appears  and  the  planar 
growth  may  become  unstable.  The  criterion  may 
easily  be  extended  to  account  for  the  thermal 
diffusion  transport  of  solute  (i.e.,  the  cross  diffu¬ 
sion  of  solute  driven  by  the  temperature  gradient). 
Consider  a  planar  crystal-fluid  interface  with  the 
fluid  consisting  of  an  inert  (carrier)  component 
and  a  dilute  component  which  forms  a  solid  de¬ 
posit  with  total  rejection  of  the  carrier  species. 
The  advance  of  the  interface  is  linked  to  the 
diffusive  arrival  of  depositing  material,  resulting 
in 


/  d<a  io  dT  \ 

(I) 

where  pj  is  the  solid  deposit  density,  the 

interface  advance  velocity,  p  the  fluid  density, 
oj  =  p|/p  the  solute  mass  fraction,  T  the  temper¬ 
ature  and  Y  the  coordinate  normal  to  the  inter- 
The  right-hand  side  in  eq.  (1)  stands  for  the 
mass  diffusion  flux  and  it  is  eomprised  of  the 
Fickian  diffusion  flux  (with  D  being  the  binary 
diffusion  coefficient)  plus  the  thermal  diffusion 
contribution  (with  Oj  being  the  thermal  diffusion 
factor).  In  the  high  dilution  limit,  the  solute  mass 
fraction  relates  to  its  partial  vapor  pressure  p,  by 
(o  =  (Mf/M2)p^./p~  where  A/,  and  Af,  art  the 
molar  masses  of  the  solute  component  and  car¬ 
rier  species,  respectively,  and  p  is  the  total  pres¬ 
sure.  The  equilibrium  partial  vapor  pressure  over 
a  flat  solid  is  temperature  dependent  and  assum¬ 


ing  that  the  total  pressure  is  constant,  we  can 
write 


dw"*  M.  dp"**  dT  AH  dT 

— —  *  *  __  —  ___ 

dT  Af,p  dT  dT”'"  RT'dT’ 


where  the  Clausius-Clapcyron  relation  has  been 
used  and  the  solid  molar  volume  has  been  ne¬ 
glected.  AH  is  the  latent  heat  per  mole  of  de¬ 
positing  species  converted  to  solid  and  R  is  the 
universal  gas  constant.  When  the  mass  fraction 
gradient  at  the  interface  is  equal  to  this  equilib¬ 
rium  mass  fraction  gradient,  we  obtain 


pD- 


AH  dT 

-  -f  Q,  - 

RT  'IdY 
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For  a  constant  temperature  gradient,  an  increase 
of  the  crystal  growth  velocity  beyond  this  value 
will  produce  a  mass  fraction  gradient  at  the  inter¬ 
face  larger  than  the  equilibrium  mass  fraction 
gradient  indicated  in  eq.  (2).  Then,  the  fluid 
ahead  of  the  front  becomes  locally  supersatu¬ 
rated.  This  local  supersaturation  is  a  necessary 
condition  for  the  appearance  of  morphological 
instability.  Thus,  for  any  growth  velocity  slower 
than  the  value  given  by  eq.  (.^)  the  planar  inter¬ 
face  remains  stable,  eq.  O)  shows  that  the  range 
of  stable  growth  velocities  increases  when  the 
thermal  diffusion  factor  takes  on  positive  values. 
In  fact,  is  positive  for  massive  species  in  light 
carrier  gases,  which  is  a  usual  characteristic  of 
fluid  mixtures  in  crystal  growth  from  the  vapor 
phase. 

As  will  be  shown  below,  this  oversimplified 
analysis  provides  only  a  rough  lower  bound  for 
the  instability  threshold.  Dissipative  effects  as 
well  as  the  effects  of  latent  heat  release  and 
surface  tension  arc  stabilizing  mechanisms  which 
may  delay  the  instability  onset  well  beyond  this 
local  supcrsaturation  criterion.  Nevertheless,  al¬ 
though  the  planar  interface  may  remain  morpho¬ 
logically  stable  due  to  these  stabilizing  mecha¬ 
nisms.  the  existing  supersaturation  may  lead  to 
nucleation  phenomena  on  the  fluid  side.  Under 
these  circumstances,  depositing  material  arrives 
to  the  interface  by  diffusion,  as  indicated  in  eq. 
(1).  and  also  by  thermophoretic  transport  of  the 
nucleated  particles.  Even  if  crystalline,  these  par- 
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tides  would  deposit  at  random  orientations  onto 
the  interface  and  reduce  the  quality  of  the  grow¬ 
ing  crystal.  This  combined  transport  of  nucleated 
particles  and  vapors  was  analyzed  by  Castillo  and 
Rosner  [9],  although  in  a  different  physical  prob¬ 
lem;  namely,  the  deposition  of  vapors  on  to  cold 
surfaces  across  laminar  thermal  boundary  layers. 
The  velocity  given  by  eq.  (3)  provides,  in  practice, 
an  upper  limit  below  which  no  nuclcation  can 
tKcur.  In  any  ease,  we  will  not  consider  further 
the  prrssibility  of  homogeneous  nucicaiion  in  the 
fluid  side. 

The  case  of  crystal  growth  from  the  tapor 
phase  presents  some  other  particularly  interesting 
features.  Due  to  the  high  density  difference  be¬ 
tween  gas  and  solid  phases,  the  growth  process 
can  usually  be  taken  as  quasi-stationary.  More¬ 
over.  in  dilute  systems,  the  Stefan  flow  induced 
by  the  advancing  front  is  negligible  compared  to 
the  diffusive  transport  and  fluid-dynamics  often 
plays  a  secondary  role;  .see,  c.g..  Rosner  (10].  For 
condensing  macromolccules  dilute  in  a  light  car¬ 
rier  gas,  an  interesting  limiting  case  is  the  so- 
called  quasi-Lorentzian  gas  mixture,  analyzed  by 
Mason  [llj.  for  which  the  Schmidt  number  Sc 
(ratio  of  momentum-diffusivity  to  mass-diffusiv- 
ity)  and  the  thermal  diffusion  factor  take  on 
very  large  values,  but  their  ratio  remains  of  order 
unity  (e.g.,  «|/Sc  tends  to  3/4  for  hard  spheres 
undergoing  elastic  specular  collisions,  sec  also 
Garcia- Ybarra  and  Rosner  [12]).  Therefore,  when 
such  a  gas  mixture  is  subjected  to  a  temperature 
gradient,  thermal  diffusion  (Sorct  “cross"  mass 
transport)  becomes  the  dominant  mass  trans¬ 
fer  mechanism,  overwhelming  the  more  familiar 
Fickian  (concentration)  diffusion  everywhere  in 
the  bulk  (although  the  latter  must  be  retained  in 


*'  Due  to  the  Onsattcr  recipriKal  relation  (see,  e.g..  Pc 
Gr(X)t  and  Ma/ur  [l.t])  use  of  the  simple  Fourier  law  for 
heat  diflusion  could  apparently  he  questioned,  hut  the 
ratio  of  Dufour  cross  heat  transfer  contribution  to  the 
Fourier  heal  diffusion  is  proportional  not  only  to  but 
also  to  the  mass  fraction  of  the  heavy  component.  Under 
conditions  of  high  dilution,  the  Dufour  effect  can  there¬ 
fore  be  neglected. 


the  vicinity  of  vapor/solid  interfaces  to  fulfill  the 
required  boundary  conditions). 

In  two  recent  papers  [14,15],  the  authors  have 
shown  than  in  cry.slal  grovth  from  quasi- 
Lorentzian  gas  mixtures  (when  Fickian  mass  dif¬ 
fusion  plays  a  minor  role)  there  exists  a  new  type 
of  morphological  instability.  In  these  mixtures, 
the  mass  deposition  rate  is  directly  proportional 
to  the  temperature  gradient  normal  to  the  inter¬ 
face.  When  the  thermal  conductivity  of  the  solid 
is  larger  than  that  of  the  fluid  mixture,  a  defor¬ 
mation  of  the  solid-gas  interface  induces  a  defor¬ 
mation  of  the  gas  phase  i.sotherms  that  increases 
the  local  thermal  gradients  and  therefore  in¬ 
creases  the  local  mass  deposition  rate  at  the 
cicsts  and  reduces  them  in  the  valleys.  Thus,  any 
initial  surface  deformation  is  enhanced  by  the 
modified  deposition  rates.  Planar  growth  is  then 
always  unstable  (only  the  large  wavenumber  dis¬ 
turbances  are  stabilized  by  surface  energy  effects) 
and  the  system  tends  to  present  as  large  an 
interfacial  area  as  po.ssible.  This  is  the  morpho¬ 
logical  instability  induced  by  thermal  mass  diffu¬ 
sion.  In  the  opposite  case  (A/Aj)>  1  (a  rather 
unphysical  situation),  the  planar  front  would  be 
linearly  stable.  However,  latent  heat  effects  were 
not  considered  in  our  previous  analyses.  This 
heat  release  is  proportional  to  the  local  mass 
deposition  rate  and  will  play  a  stabilizing  role  by 
increasing  locally  the  temperature  at  the  crests 
(thus  reducing  the  temperature  gradients  and 
therefore  the  local  deposition  rates)  and  decreas¬ 
ing  the  local  temperature  at  the  valleys  (increas¬ 
ing  temperature  gradients  and  deposition  rates). 

In  this  paper  we  study  morphological  instabili¬ 
ties  that  may  appear  during  planar  growth  of  a 
solid  deposit  generated  from  a  quasi-Lorentzian 
vapor  mixture.  The  condensate  is  assumed  to  be 
comprised  exclusively  of  the  heavier  component, 
which  drifts  towards  the  gas-solid  interface  due 
to  the  thermal  diffusion  induced  by  some  exter¬ 
nally  imposed  thermal  gradient.  Latent  heat  re¬ 
lease  will  be  taken  into  account,  as  well  as  the 
influence  of  surface  energy  and  Fickian  diffusion 
for  large  wavenumber  disturbances.  Hurle  [16] 
presented  a  similar  analysis  of  the  opposite  limit; 
that  is,  the  influence  of  thermal  diffusion  on  the 
morphological  instability  of  a  crystal  growing 
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mainly  by  Fickian  solute  diffusion.  Thus,  this 
paper  was  f(K'used  on  the  influence  of  thermal 
diffusion  on  the  supcrsaturation-induccd  instabil¬ 
ity.  but  its  effects  were  assumed  to  be  rather 
small  and  he  did  not  consider  the  thcrmophorcti- 
cally  induced  instability  emphasized  here. 

This  paper  is  organized  as  follows.  In  section  2 
we  establish  the  governing  equations  for  the  liKal 
fields  of  velocity,  temperature  and  mass  fraction 
of  the  depositing  species  in  the  fluid  phase  as 
well  as  for  the  temperature  field  in  the  .solid 
deposit,  using  a  reference  system  in  which  the 
solid  moves  al  the  constant  velocity  Fj.  In  section 
?  we  detail  the  boundary  conditions  that  must  be 
satisfied  at  the  moving  vapor-solid  interface.  In 
section  4  the  equations  and  boundary  conditions 
are  rendered  dimensionless  and  the  relevant  di¬ 
mensionless  parameters  are  identified.  The  fluid 
Peclet  number  (based  on  the  fluid  velocity  for  a 
planar  solidification  front)  is  shown  to  be  a  small 
I. umber  (proportional  to  the  mass  fraction  of  the 
dilute  depositing  species)  and  an  expansion  in  the 
Peclet  number  is  formally  carried  out.  To  leading 
order,  convective  as  well  as  unsteady  terms  arc 
negligible  in  the  governing  equations  compared 
to  the  diffusive  terms.  Since  the  Schmidt  number 
is  a  large  parameter,  a  matched  asymptotic  ex¬ 
pansion  in  the  inverse  of  the  Schmidt  number  is 
performed.  As  done  by  Goren  [17],  we  consider 
an  outer  region  far  from  the  interface  where 
Fickian  (Brownian)  diffusion  can  be  neglected 
and  an  inner  region  (a  Brownian  sublayer,  see 
Gokoglu  and  Rosner  [18])  adjacent  to  the  solid- 
vapor  interface,  where  concentration  diffusion 
becomes  important,  bringing  the  mass  fraction 
from  the  value  required  to  matcb  with  the  outer 
region  to  the  value  required  by  the  boundary 
conditions.  The  temperature  and  mass  fraction 
fields  (in  both  outer  and  inner  regions)  in  the 
stationary  state  of  planar  growth  are  given  in 
section  5.  The  remainder  of  the  paper  is  devoted 
to  a  study  of  the  stability  of  this  stationary  state 
with  respect  to  small  unavoidable  disturbances.  A 
linear  stability  analysis  is  provided  in  section  6, 
giving  a  dispersion  relation  linking  the  disturbance 
amplification  factor  with  the  corresponding 
wavenumber.  Since  this  analysis  breaks  down 
when  the  disturbance  wavelength  is  comparable 


to  the  Brownian  sublayer  thickne.ss,  it  is  supple¬ 
mented  in  section  7  to  take  into  account  the  ease 
of  very  large  wavenumber  disturbances.  Section  8 
contains  remarks  on  the  implications  of  our  re¬ 
sults.  necessary  extensions,  and  general  conclu¬ 
sions. 


2.  Governing  equations 

We  consider  a  binary  fluid  mixture  with  a 
depositing  species  (density  p,  and  velocity  V,) 
and  a  carrier  component  (density  P2  and  velocity 
F;).  The  total  density  is  p  =  P|  -i-p;  and  the  mix¬ 
ture  barycentric  velocity  V  =  (p,F,  p^F^I/p. The 
depositing  (active)  component  drifts  by  the  com¬ 
bined  effect  of  Fickian  and  thermal  diffusion 
towards  a  solid  surface  where  a  solid  deposit 
forms  with  total  rejection  of  the  carrier  species. 
The  governing  equations  for  the  fluid  mixture, 
considered  incompressible  and  with  constant 
thermophysical  properties,  are  (see.  e.g..  the  book 
by  Rosner  [10]): 


r  •  F  =  0. 

(4) 

ftF 

p—  -t-pF-  rF=  -rp-t-pr-F, 

:5) 

pc,,—  -(-pc,.F-  rT=  AT-T. 

(f ' 

dp, 

— -Hr(p,F,)  =  (). 

(7) 

In  the  above  equations,  p  is  the  dynamic  vis¬ 
cosity,  the  specific  heat  of  the  fluid  mixture  at 
constant  pressure,  and  A  the  Fourier  thermal 
conductivity.  The  mass  flux  of  the  depositing 
species  relates  the  velocity  F,  to  the  barycentric 
velocity  F  and  the  concentration  and  thermal 
gradients.  This  relation,  as  given  by  De  Groot 
and  Mazur  [1.^]  (see  also  Rosenberger  [19]).  is 

rr 

pwF,  =  pwF  —  pDTw  —  paT-Dwf  1  —  ft)) -y- .  (8) 

The  first  term  on  the  right  hand  .side  is  the 
convective  flux,  the  second  term  represents  the 
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direct  Fickian  diffusion  flux  and  the  third  term 
represents  the  thermal  diffusion  (Sorct)  effect. 

For  the  temperature  field  in  the  solid  deposit, 
the  governing  equation  is 

'*7ii 

=  Aj,r'7j.  (d) 

where  the  subscript  "d"  refers  to  the  deposit  and 
all  symbols  have  their  usual  meaning. 

The  above  set  of  equations  must  be  solved 
subject  to  appropriate  boundary  conditions  on 
the  fluid  and  solid  sides  and  at  the  moving 
boundary  between  the  two  phases. 


3.  Boundary  conditions  at  the  solid-fluid  inter¬ 
face 

For  simplicity  we  restrict  ourselves  to  a  two-di¬ 
mensional  geometry  iX.Y)  and  consider  a  refer¬ 
ence  system  in  which  the  solid  deposit  moves  with 
constant  velocity  f  j  directed  along  the  Y  axis. 

The  separation  line  between  the  fluid  mixture 
and  the  solid  deposit  is  defined  by  the  position 
vector  r.: 

r,  =  (A-.i{T.  r)).  (10) 

The  locus  is  an  unknown  of  the  problem,  as  in 
other  problems  with  free  moving  boundaries  (for 
example,  combustion  fronts  in  Clavin  and 
Garcia- Ybarra  [20],  liquid-gas  interfaces  in 
Garcia-Vbarra  ct  al.  [21],  or  solidification  fronts 
in  Langcr  [ft]).  Then,  the  unit  normal  vector 
pointing  towards  the  fluid  side  and  the  unit  tan¬ 
gent  vector  arc 


(II) 

The  normal  component  of  the  local  intcrfacial 
velocity  is  related  to  the  evolution  of  the 


interface  location,  f(X.  t).  via  the  kinematic  con¬ 
dition: 

-=A/K-n.  (12) 

We  as.sume  that  the  inert  component  of  the 
fluid  mixture  docs  not  penetrate  the  deposit,  that 
is,  it  moves  with  the  same  speed  as  the  interface. 
Therefore: 

F,/i  =  K-n.  (13) 

Moreover,  the  no-slip  condition  at  the  interface 
corresponds  to  continuity  of  tangential  velocities, 
that  is: 

F/=Kj/.  (14) 

On  the  other  hand,  continuity  of  mass  flux  across 
the  interface  implies: 

p(y-  K)  n  =  pa>(r,  - »;)  n  =pj(ij-  rj  •«. 

(15) 

By  using  relations  (12).  (13)  and  (K)  in  bi>undary 
condition  (l.S).  we  obtain: 

i'C  p.V  ID  VT\ 

‘•t  ‘  p,,  \\- to  T  I 

(1ft) 

and  from  eqs.  (12)  and  (1.^). 

Pj  I  PkI  \ 

F-nN= —F., •«/)/+  1  - (17) 

P  \  p  I  til 

Note  that  the  unknown  I’,  can  be  removed  by 
ch(X)sing  a  reference  system  in  which  the  solid  is 
at  rest  and  the  solid-fluid  interface  is  moving.  In 
that  ca.se.  Fj  =  0.  but  this  is  not  the  most  conve¬ 
nient  reference  system  because  then  there  is  no 
steady  state  solution.  Instead,  we  use  a  reference 
system  linked  to  the  interface  in  which  the  solid  is 
retreating  with  the  crystal  growth  velocity.  Once 
has  been  specified  by  a  proper  selection  of  the 
reference  system,  cq.  (Ift)  determines  the  evolu¬ 
tion  of  the  fluid-solid  interface  IcKation  and  the 
boundary  conditions  (14)  and  (17)  give  the  fluid 
vcIcKity  (tangential  and  normal  component,  re¬ 
spectively)  at  the  interface. 
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Regarding  the  asst»eiated  thermal  conditions, 
we  consider  that  the  temperature  and  the  energy 
flux  are  continuous  across  the  interface,  so  that; 

T=  7',. 


M, 


where  l„  fi  is  given  by  (12)  and  iAH/M,)  is  the 
latent  heat  per  unit  mass  converted  into  solid. 

Finally,  local  thcrmiKhemical  equilibrium  at 
the  fluid-solid  interface  will  be  assumed,  i.e.: 

w{  X.  t )  =  /',  local  curvature) .  ( Id) 

where  the  functional  dependence  of  is  gov¬ 
erned  by  thermodynamic  considerations.  We  a.s- 
sume  that  to"'*  depends  »»nly  on  the  l(Kal  interface 
curvature  and  temperature,  see,  for  instance,  the 
biHik  by  Frenkel  (22). 


4.  Dimensionless  equations  and  boundary  condi¬ 
tions 

At  this  point  it  is  quite  useful  to  normalize  and 
make  dimensionless  the  variables  of  the  problem, 
by  defining  r  s  ,v  s  X/L,  v  s  Y/L,  i-  = 

F/i;,,.  «=7/7,,,.  (%^T,/T,,.,. 
and  Tjl.v,  7)  =  4{A',  i)/L.  where  the  reference 
values  will  be  selected  later.  Also,  we  define  the 
dimensionless  distance  from  the  interface: 

2  =  y-T}{.v,T)  (20) 

and  rewrite  our  equations  in  a,  z  coordinates. 
Then,  the  governing  equations  on  the  fluid  side 
(z  >  0),  cqs.  (4)  to  (7),  transform  to 

rr  =  0,  (21) 

,  I  (tr  ;>TJ  ftr  \ 

Pr  '  Pc  —  -  ^ - -(- 1-  ■  Tr  =  -  n/  -I-  r-r. 

\  dr  ftr  dz  ) 


I  dW  dri  dW  \ 

Pe  - - -i-r  - 

\  dr  dr  dr  I 


,  /  d/7  d77  d/7  ) 

Pr  '  Pc  ^ - -H'  r/7 

\  dr  dr  dz  / 

,  (24) 

=  —  r-/7  4ar[/7(l  -<o,,,/7)r  in  «), 

where  cq.  (8)  has  been  used  in  cq.  (7)  to  obtain 
cq.  (24). 

Whereas,  in  the  solid  (z  <  0),  from  cq.  (9),  wc 
find; 


Pe  CiliA 
\  Ki 

= 


d(9.|  dT>  df-7.| 

dr  dr  dz 


Wc  have  introduced  the  dimensionless  pressure 
77  ^  and  the  dimensionless  parameters 

Pc  =  F',^.,Z./(A/(pc^,)]  (thermal  Peclet  number), 
the  Prandtl  number  Pr  =  i'/[A/(pr,,)].  the 
Schmidt  number  Sc  =  v/D.  the  parameter  a  = 
a-fD/v  that  determines  the  importance  of  ther- 
mophorcsis  (Sorct  effect),  and  v  is  the  kinematic 
viscosity  (momentum  diffusivity),  p/p. 

Note  that  due  to  our  selection  of  the  x.  z 
coordinate  system,  the  operators  T  and  T-  are. 
respectively,  given  by; 


/  d  drj  d  d  \ 

\  d.v  d.v  dz  *  dz  /  ' 


d-Tj  d  ^dT7  d- 

dA  '  dz  dA  dA'dz 


/dry)-  d- 
+  1  +  —  — 
d.v  dz  - 


Therefore  the  fluid-solid  interface  is  located  at 
z  =  0.  and  the  following  boundary  conditions  must 
be  fulfilled; 


r/  =  r.,/. 


\  P  I  P  ! 


'■*9  P  "refPr 

—  =rjnN-l- - - — Nn 

dr  Pj  Pe 


•  —  r/7-i-/7(l  -w,„/2)rtr  In  (9  .  (29) 


(2.1) 
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«  =  (30) 

Pol  \ 

-  Nn.  ■  V  (-)  —  I  •  nN - I 

p  \  “  dr) 

=  -yNiir«,.  (31) 

/2  =  /2'‘’{6>) -y/C,  (32a) 

with  the  equilib¬ 

rium  mass  fraction  over  a  flat  surface,  K  the  local 
curvature  and  y  a  dimensionless  parameter  that 
takes  into  account  the  interfacial  energy.  Using 
the  Gibbs-Kelvin  relation  (see,  for  instance,  the 
bot>k  by  Frenkel  [22]),  this  parameter  can  be 
written  as 

y  =  2o-T/W,/pj/?rw„,L,  (32b) 

where  (Tj  is  the  interfacial  tension  of  the  solid- 
fluid  interface. 

As  in  the  problem  of  directional  solidification, 
a  stationary  state  can  be  achieved  by  imposing 
the  remaining  boundary  conditions  at  constant 
distances  from  the  fluid-solid  front,  let  us  say,  at 
a  distance  L  on  the  fluid  side  (that,  is  at  y  =  1) 
and  at  a  different  distance  Z.j  on  the  solid  side 
(that  is,  at  y  =  -6,  where  8  s  L^/L).  Therefore, 
we  assume  that  the  mass  fraction  at  y  =  1  is  kept 
fixed  and  the  temperature,  at  y  =  1  and  y  =  -5, 
reaches  some  particular  values  during  the  station¬ 
ary  state.  That  is 

/2{y  =  1)  =f2,  = 6>{y  =  1}  =0,-l-nst, 

"ref 

0{  y  =  -6)  =  0J5  +  nst, 

(33) 

with  the  subscript  e  denoting  values  at  y  =  1.  and 
subscript  8  values  at  y=  -8.  The  notation  nst 
stands  for  nonstationary  terms  corresponding  to 
disturbances  of  the  stationary  state.  For  these  nst, 
Biot  type  boundary  conditions  will  be  used  and 
discussed  later  on  (see  eqs.  (58)  and  (59)). 

As  mentioned  before,  is  only  a  function  of 
the  reference  system.  In  the  case  of  a  planar 
growing  interface,  a  useful  reference  system  is 
that  in  which  the  interface  docs  not  move  at  all 
and  is  always  located  at  y  =  0.  Then  d-q /dr  =  0 


and  Cj  is  found  to  compensate  the  otherwise 
secular  growth  of  q  given  by  boundary  condition 
(29).  In  that  ca.se,  from  eq.  (28)  and  taking  into 
account  that  r  is  of  order  unity,  and  cq.  (29),  it  is 
clear  that  Pc  ~  Pr  In  many  practical  cases 
the  depositing  species  is  very  dilute  in  the  carrier 
fluid,  so  that  is  quite  small  and  therefore 
Pe  ■«  1  provided  that  the  Prandtl  number  Pr  is  of 
order  unity  (as  in  gases).  To  leading  order  in  a 
small  Pecict  number  expansion,  the  physical 
problem  described  by  eqs.  (21)  to  (25)  with 
boundary  conditions  (27)  to  (33)  becomes  a  purely 
diffusive  (negligible  advection)  problem  and  the 
flow  field  in  the  gas  plays  no  role.  In  the  follow¬ 
ing  we  restrict  ourselves  to  this  case  (i.e.,  «*:  1 

and  Pe  <*:  1 ).  Note  that  we  have  incorporated  the 
Peclet  number  into  the  definition  of  the  latent 
heat  parameter  y,  because  we  arc  interested  in 
emphasizing  the  stabilizing  effects  of  latent  heat 
release  on  the  thermophoretically  induced  insta¬ 
bility,  as  mentioned  in  the  introduction. 

We  focus  here  on  cases  in  which  the  main 
transport  mechanism  of  the  depositing  species 
toward  the  solid  is  thermal  (Soret)  diffusion,  with 
Fickian  diffusion  playing  a  minor  role.  This  situa¬ 
tion  corresponds  to  high  values  of  the  Schmidt 
number.  Then  Fickian  diffusion  effects  become 
important  only  in  a  very  thin  region  (the  Brown¬ 
ian  sublayer)  surrounding  the  fluid-solid  inter¬ 
face.  This  sublayer  allows  the  mass  fraction  pro¬ 
file  to  satisfy'  boundary  condition  (32a).  To  deal 
with  this  situation,  we  develop  a  boundary  layer 
analysis,  as  indicated  in  the  book  by  Van  Dyke 
[23].  Thus,  we  define  an  inner  region  near  the 
interface  using  the  stretched  variable; 


^  =  Sc-2.  (34) 

Therefore,  in  the  limit  Pc  ■«  1,  eqs.  (21)  to  (24) 

transform  to 

rT  =  0.  (35) 

r-r  =  r//.  (36) 

r=0  =  O.  (37) 

^r-/?-Har(/}/0)  •  r0  =  o.  (38) 

r-0j  =  ().  (39) 


n: 
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with  the  operators  V  and  T’.  in  the  outer  region 
given  hy  (26a)  and  (26h).  respectively.  Whereas  in 
the  inner  region  (Brownian  sublayer)  they  are 
given  by: 


r- 


<<  ftn  (I 

—  -  Sc—  — 

i'.V  (tX 


r-  = 


(*■  'd'-q 

—  -  Sc— 

ll.v  ■  <*.V  • 


(t  flT)  ft' 

- 2Sc— : - 

HA'  (tA<>^ 


+  Sc- 


(4()a) 


(4()b) 


The  boundary  conditions  at  f  =  0  are  still  given 
by  eqs.  (27)  to  (.'^2)  with  the  terms  in  (2*)) 

being  negligibly  small. 


5.  Stationary  state 


Under  the  above-mentioned  circumstances, 
there  exists  a  steady  solution  with  a  planar  inter¬ 
lace  moving  at  constant  velocity.  We  choo.se  our 
reference  system  such  that  the  planar  interface 
does  ni)t  move  and  remains  at  y  =  0.  Therefore 
7j  =  (I  and  r  =  y. 

Then  the  solution  of  eq.  (.'^7)  and  eq.  with 
boundary  conditUms  (.lO)  and  (.Tib)  may  be  writ¬ 
ten: 


(41) 

(42) 

with  and  (i  and  the  dimension¬ 

less  temperature  gradients  in  the  fluid  and  solid 
phase,  respectively.  At  this  point  we  select  the 
reference  velocity  f'^.,  as  the  fluid  velocity  during 
the  stationary  growth,  then  r"  =  (().  -I)  and  r,, 
=  (p/p,, )/•''.  Therefore,  from  boundary  condi¬ 
tions  (.11)  and  (.1.1c),  we  obtain  the  dimensionless 
stationary  temperature  gradients 


L  dT"  (■),-(■),- {  A /\,) 
T,,,  iiY  1+«(A/A„) 

L  dV  A  f-K  -  (% 

dv  A,  1 -(-«(A/Aj  ■ 


For  the  mass  fraction  in  the  outer  region,  to 
leading  order  in  a  Sc  expansion,  eq.  (.IK)  leads  to 
the  solution: 


/i, 

ii"  =  —  f-r' 
(•)  * 


(44) 


where  the  boundary  condition  (.1.1a)  has  been 
imposed.  In  order  to  satisfy  the  mass  fraction 
boundary  condition  at  the  solid  interface  (32a). 
the  inner  diffusion  layer  must  be  considered.  In 
this  Brownian  sublayer,  eq.  (38)  taking  into  ac¬ 
count  (40)  and  (41)  becomes  (to  leading  order  in 
Sc): 


d(- 


«  d/2:;. 


(-1 


d^ 


(4.S) 


with  (->„  =  f-2,.  - /3,  being  the  temperature  (in  di¬ 
mensionless  form)  of  the  planar  interface. 

The  solution  of  eq.  (4.S)  satisfying  boundary 
condition  (32a)  at  £  =  0.  and  matching  with  the 
outer  solution.  (/2:i{^  -♦  x)  =  /2:,u,{z  =  ()))  is: 


C 

(46) 


//::=  /2:;'--/2,—  exp  -  — 


with  /i:;**  S  /2^''^{(-2„). 

On  the  other  hand,  because  r"  =  (0.  -  1).  from 
the  boundary  conditions  (28)  and  (29)  with 
=  0.  we  find 

Pc=«u,,,Pr/2,(o)B/(-2,).  (47) 

Because  the  velocity  of  the  solid  deposit,  from 
boundary  condition  (14)  is  fj=  it 

turns  out  that 

p  1  dT" 

K,=  -aiDo>  r".  (48a) 

P.I  T- 

with  dT'Vdy  the  constant  temperature  gradient 
in  the  fluid  mixture  indicated  in  eq.  (43).  Note 
that  -  Kj  is  the  velocity  of  the  moving  .solid-fluid 
front  measured  in  a  reference  system  in  which 
the  solid  deposit  remains  at  rest,  that  is  l  ',,r,mih- 
Then,  from  eq.  (48a)  the  product  aj(dT''/d>') 
(and  therefore  ap)  is  positive  for  a  growing  solid 
(and  it  will  be  negative  for  a  subliming  solid). 


J.I..  CaMilli)  cl  at.  /  Morpholosit  al  inslahihn'  ttf  ihcrmiiphiirclically  unm  inn  (Icimimi 


113 


Wc  can  also  imagine  that  in  some  actual  exper¬ 
iments  what  is  fixed  instead  of  the  distances  /, 
and  /-j  is  the  total  distance  I  =  L  +  over  which 
the  temperature  and  mass  fractions  arc  con¬ 
trolled.  In  that  case,  from  cqs.  (4.^a)  and  (48a) 
depending  on  the  crystal  growth  vcliK'ity.  the 
location  of  the  solid-fluid  interface  can  be  ob¬ 
tained;  i.c..  the  relative  value  of  L.  One  finds 


(J//  (tjl)pM^.\  p  ■ 

at;  )  ^ 

A  /  J//  (I  I  lipui^.  \ 


The  thickness  of  the  fluid  layer  may  vary  from  the 
maximum  value  L  =  I  to  its  minimum  value  L  =  0. 
These  two  extreme  values  provide  limiting  crystal 
growth  vcliKilics.  resulting 


p  a  I  Do)  T  —  T 
Pd  '  L 

y  ( /  =  (u  =  Ai  _ 

AH  ayDpio- 

j  + - -  -- 

M,  AT. 


This  result  can  be  used  as  a  quasistcady  ap¬ 
proximation  to  the  nonstationary  growth  when 
the  concentration  and  temperatures  arc  kept  con¬ 
stant  at  fixed  locations.  The  growth  is  assumed  to 
accommodate  at  every  instant  to  the  new  inter¬ 
face  location,  with  a  growing  velocity  varying  from 
=  /}.  at  the  beginning,  when  no  solid  is 
present,  to  the  value  =  0},  when  the 

solid  occupies  the  entire  cavity. 

We  are  now  interested  in  finding  under  what 
conditions  the  solution  given  by  eqs.  (41),  (42). 
(44)  and  (46)  is  stable  and  when  an  initially  in¬ 
finitesimal  disturbance  will  grow  to  a  finite  ampli¬ 
tude. 


6.  Linear  stability  analysis 

Wc  assume  that  stationary  planar  growth  is 
perturbed,  so  that  the  interface  deforms  and  the 
fields  become  H  =  0"  +  0,  =  0^'  +  0^,  = 


+  /ou.  and  When  the  ampli¬ 

tudes  of  the  disturbances  are  infinitesimally  small, 
the  corresponding  governing  equations  and 
boundary  conditions  can  be  linearized  around  the 
initial  undisturbed  state,  resulting  from  eqs.  (M) 
to  (.19)  with  r  and  T’  given  by  cqs.  (26) 


whereas  in  the  Brownian  sublayer,  from  eqs.  (.18) 
and  (40) 


J!_r  + 


0{z  =  o}  w:'„  g /ctg\ 

f-)„  \  i»r  ) ; (t£ 

(52) 


The  boundary  conditions  for  the  temperature  dis¬ 
turbances  at  the  interface  are 


0  =  0^ 

at  2  =  0, 

(-5-^) 

»0  _  Ad 

Hz  p  Hr  A  Hz 

at  2  =  0. 

(54) 

The  matching  condition  between  the  inner  and 
the  outer  mass  fraction  disturbances  corresponds 
to 

=/;„,{  2 -»0).  (55) 

Linearizing  cq.  (29)  with  respect  to  the  distur¬ 
bances  and  taking  into  account  (47),  we  obtain 


drj  p  0^ 

ilT  pj  Si^ap 


a/3  a 
+  — 


at  £  =  2  =  0, 


(56) 
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and  fn)m  boundary  condition  (32a) 

'*’’7 

at  ^  =  0,  (57) 

<>.v ' 


by  expanding  each  variable  in  a  Fourier  series.  As 
the  system  is  linear,  we  can  study  the  behavior  of 
the  stationary  state  (41),  (42),  (44)  and  (46)  via 
each  Fourier  mode.  Therefore,  we  consider 


with 

^  I  df-> 

the  parameter  that  measures  the  variation  of  the 
equilibrium  mass  fraction  with  temperature,  see 
eq.  (2). 

F-inally.  we  must  specify  the  disturbances 
boundary  conditions  at  the  fluid  side,  y  =  1  (i.e., 
c  =  I  -  T/f.v.  t))  and  at  the  solid  side,  y  =  —fi 
(i.e.,  e  =  -<S  -  t;{.v.  t}).  For  the  temperature,  we 
consider  Biot  type  boundary  conditions,  that  is. 
the  local  heat  flux  perturbation  is  taken  to  be 
proportional  to  the  local  temperature  field  distur¬ 
bance.  then,  we  obtain  the  linearized  conditions 

—  =  -  Bi(  -/3t7  +  0)  at  z==  1,  (58) 

d: 

^  =  Bij(-p,T,  +  e,)  at  z=-8.  (59) 

nr 

The  terms  proportional  to  -q  arise  from  the  ap¬ 
parent  displacement  of  the  boundaries  with  re¬ 
spect  to  the  fluid-solid  interface.  The  Biot  num¬ 
ber  at  the  upper  (Bi)  and  lower  (Bij)  boundaries 
arc  defined  positive  and  measure  the  thermal 
resistance  of  the  boundaries.  The  limit  of  a  very 
large  Biot  number  corresponds  to  a  perfectly 
conducting  boundary  where  the  temperature  re¬ 
mains  constant  and  thermal  disturbances  vanish. 
The  opposite  limit  of  very  small  Biot  numbers 
corresponds  to  an  thermally  insulating  boundary 
where  the  heat  flux  is  kept  fixed. 

The  mass  fraction  is  assumed  to  remain  con¬ 
stant  at  y  =  1.  therefore 


F{ jr,  z,  t}  =  F{ 2}  cxp(it»- +  o-t).  (61) 

where  F  stands  for  0,  6^.  r,„  and  tj  (this  last 

one  without  the  2  dependence).  By  using  (61)  in 
eq.  (49).  we  find  that  the  solution  is  given  by 

0  =  ^7} -I-/4  c*-  -i-F  e'*-.  (62) 

From  Eqs.  (50)  and  (51)  we  obtain,  respectively. 

»d  =  i3aV)-h/l,e‘--  +  Bje-‘--.  (63) 

77,  72,  ,  77,  , 

(64) 


Eq.  (52)  may  be  integrated  once,  yielding 


d^  ■"  0„  0„ 

,  Pj  77,0^ 


(65) 


where  the  second  term  on  the  right-hand  side  of 
the  equation  is  an  integration  constant  which  has 
been  obtained  by  using  boundary  condition  (56). 

The  solution  of  eq.  (65)  with  given  by  eq. 
(46)  is 


.  77, 


\  P 


0„ 


+  -  +/!—-- 


1 


1  A: 
+  B\  —  -I-  - 

P 


“  0., 


00  P 
0. 


77, 


at  2  =  1. 


(60) 


where  we  have  used  the  value  of  77'*  given  by  eq. 
(44). 

The  system  of  eqs.  (49)  to  (52)  with  boundary 
and  matching  conditions  (53)  to  (60)  can  be  solved 


+  B\-^+k 


^  }  exp 


a 


0, 


Pi 


(66) 
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The  values  of  A,  B.  A^,  flj,  C,  E  and  t)  must 
satisfy  the  boundary  eonditions  (58)  and  (60)  at 
;  =  1,  (59)  at  z=  -d,  the  matching  condition 
given  by  (55),  and  the  conditions  (53),  (54)  and 
(57)  at  z  =  0.  By  substitution  of  (62),  (63),  (64) 
and  (66)  in  the  boundary  conditions,  the  following 
homogeneous  system  of  algebraic  equations  is 
obtained; 


/((A  +  Bi)  e* +  B(Bi -A)  c'*  =  0,  (67) 

02 

/I  c* +B  e'* +  =0,  (68) 

/lj(Bi,-A)  e-**  +  Bj(Bij  +  A)  e*^  =  0,  (69) 

^  ^  .  P  d 

(70) 

a+b-a.-b.  +  mb-BJ-o,  (71) 

^d  Pd  ^ 

A-B-  -V  =  0,  (72) 

A  A  p  K 


k\n. 


'0  B )  0c- 


0  n  \  jOi  ^*1  ^ 


+  B 


+  V 


n 


-0, 


n. 


<^Pd 


0„ 


-  yk^ 


+  £  =  0. 


(73) 


Note  that  eq.  (73)  is  needed  just  to  specify  the 
value  of  £  which  determines  the  inner  mass 
fraction  disturbance  by  eq.  (66).  This  value  of 
£  is  such  that  the  equilibrium  condition  at  the 
solid  interface  is  fulfilled.  As  a  consequence,  the 
Brownian  sublayer  has  no  influence  at  all  on  the 
stability  of  the  planar  front  as  long  as  the  distur¬ 
bance  wavenumbers  are  not  too  large  (see  next 
section).  No  matter  how  complicated  the  mass 
fraction  boundary  condition  is  at  the  fluid-solid 
interface,  Fickian  diffusion  will  bring  the  mass 
fraction  from  its  value  at  the  outer  edge  region  of 
the  Brownian  sublayer  to  the  required  value  at 
the  interface. 

Thus,  the  value  of  a  given  below  is  not  af¬ 
fected  by  the  boundary  condition  corresponding 


li.s 


to  the  mass  fraction  at  the  interface,  eq.  (32a). 
The  system  (67)  to  (72)  is  a  homogeneous  alge¬ 
braic  system  of  six  equations  with  six  unknowns 
{A,  B,  /4j,  Bj,  C  and  -q).  The  determinant  of  the 
coefficients  must  vanish  in  order  to  have  a  non¬ 
trivial  solution.  Imposing  this  condition,  the  fol¬ 
lowing  dispersion  relation  is  obtained: 


p  p-Bii ,  1 

'  P  ,  V 

a  = 

A 

—  -t-  A  smh  A  -H  Bi  cosh  A 

Pd  B  ' 

.0c  1] 

X 


I  ^ 

{ A  cosh  A  -I-  Bi  sinh  A  -t-  — 


X 


A  sinh  A  -h  Bi  cosh  A ) 


A  cosh  kS  +  Bij  sinh  AS 
A  sinh  AS  -t-  Bij  cosh  AS 


(74) 


Recalling  eq.  (61),  the  sign  of  a  determines 
the  stability  of  the  stationary  state  against  in¬ 
finitesimal  disturbances  of  wavenumber  A.  The 
stationary  state,  with  a  planar  interface,  is  unsta¬ 
ble  when  a  >  0  and  an  initially  small  disturbance 
will  grow  until  the  non-linear  terms  become  im¬ 
portant  and  determine  the  new  shape  of  the 
advancing  front.  When  the  amplification  factor. 
O’,  is  negative  the  planar  front  is  linearly  stable 
against  infinitesimally  small  disturbances,  al¬ 
though  it  may  be  destabilized  by  finite  amplitude 
disturbances. 

The  value  of  cr  in  eq.  (74)  depends  on  the 
density  ratio  p/pj  (only  as  a  multiplicative  factor), 
the  thermal  Biot  numbers,  Bi  and  Bi,j,  the  ther¬ 
mal  conductivity  ratio  A/A^.  the  depth  ratio  8. 
and  the  dimensionless  groups 


B  dr'  p, 

0^  7^  dk'  p 

B  M,T,  A  A  dr'/dY 

(75a) 


for  quasi-Lorentzian  gas  mixtures  {oj  >  0)  all  the 
mentioned  parameters  are  positive  definite,  and 
except  for  B/^e  3’’*^  y/B<  which  depend  on 
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growth  conditions,  all  others  arc  determined  by 
material  properties.  Note  that  V//a  is  the  ratio 
of  latent  heat  release  (due  to  the  mass  deposition 
flu.x)  to  the  diffusive  heat  flux  towards  the  inter¬ 
face  from  the  fluid  side. 

For  positive  values  of  o,  (i.e.  /(>(•),  all  the 
mentioned  parameters  are  positive  and  therefore, 
the  sign  of  <r  is  the  .same  as  the  sign  of  the  group 


,  A 
H  -'a 


(7.Sb) 


When  this  group  is  positive.  <t  given  by  eq.  (74)  is 
a  monotonuolh  incraisiiin  fiinclioii  of  k.  Planar 
growth  is  always  unstable  for  perturbations  of  any 
wavenumber  and  the  system  tends  to  present  as 
large  surface  as  possible.  A  deformation  of  the 
solid-gas  interface  induces  a  deformation  of  the 
isotherms  in  the  gas  phase  that  increases  the 
thermal  gradients  and  therefore  the  local  mass 
deposition  rate  at  the  crests  and  reduces  them  in 
the  valleys.  Thus,  the  initial  surface  deformation 
is  enhanced  by  the  modified  deposition  rates.  In 
the  opposite  case.  ( ^  - )//3 )  <  0.  the  planar 
front  is  linearly  stable.  Therefore,  we  can  con¬ 
clude  that 


for  a ,  >  0.  planar  growth  is  unstable  W'hcn 


(7.Sc) 


Thus,  for  negligible  latent  heat  release  planar 
growth  is  always  unstable  because  Aj/A  >  1  un¬ 
der  usual  conditions.  Latent  heat  plays  a  stabiliz¬ 
ing  role,  as  pointed  out  in  the  introduction  of  this 
paper.  Eq.  (75c).  with  the  value  of  y/p  given  by 
eq.  (75a).  provides  a  minimum  value  of  the  group 
pw^./T^  for  stable  planar  growth  and  shows  that, 
for  sufficiently  low  dilution,  planar  growth  is  un¬ 
stable. 

The  result  indicated  in  eq.  (75c)  shows  that 
when  thcrmophoreticaily  growing  planar  deposits 
arc  stable,  from  eq.  (48c).  it  appears  that 
=  0}  >  =/}•  Therefore,  for  the  case  of 

thermal  boundary  conditions  imposed  at  fixed 
locations  and  quasistcady  growth,  the  .stable  ther- 
mophorctic  planar  growth  corresponds  to  a 
solid-fluid  front  which  starts  moving  at  the  vcloc- 


'•y-  K!«)))ii.  =  =  /)•  and  accelerates  as  the 

solid  deposit  depth  increases.  As  the  growth  vc- 
liKity  increases,  constitutional  supersaturation 
may  build  up  in  the  ncighborh(X)d  of  the  advanc¬ 
ing  interface  (i.c.,  in  the  Brownian  sublayer)  ren¬ 
dering  unstable  this  planar  growth.  To  account 
for  this  possibility  we  have  to  study  very  large 
wavenumber  disturbances  which  arc  strongly  af¬ 
fected  by  the  Brownian  sublayer  structure.  How¬ 
ever.  short  wavenumber  disturbances  arc  not  af¬ 
fected  in  the  linear  stage  by  supersaturation  ef¬ 
fects  and  the  above  results  remain  valid. 

On  the  other  hand,  the  analysis  carried  out 
here  predicts  that  the  planar  front  is  unstable 
when  a,  (and  therefore  )3)  is  negative,  and 
-(p/&\J>Bi,  as  could  be  the  case  in  some 
liquid  mixtures.  There  is  a  singular  value  of  k, 
denoted  by  k  *,  for  which  the  numerator  in  eq. 
(74)  vanishes  (that  is.  k  *  sinh  k*  +  B'\  cosh  k*  = 
-(p/(-)^).  and  thus  a{k*)  =  0).  Thus  for  (/3- 
/3j)/^)<()  (respectively.  (/3  -  ^j)//3)  >  0).  the 
planar  interface  is  unstable  for  wavenumber 
smaller  (re.spectively.  larger)  than  k  *  and  linearly 
stable  otherwise.  Therefore 


for  a,  <  0.  planar  growth  is  unstable  when 


(7.>d) 


that  is,  planar  growth  is  unstable  when  the  upper 
boundary  is  a  thermal  insulator. 

It  should  be  kept  in  mind  that  in  liquid  mix¬ 
tures  Fickian  (concentration)  diffusion  is  usually 
more  important  than  thermal  mass  diffusion  and 
the  present  analysis  must  therefore  be  extended 
to  take  this  into  account.  Ordinarily  the  velocity 
of  planar  deposit  growth  from  liquids  will  be 
determined  principally  by  direct  Fickian  trans¬ 
port.  Nevertheless  thermal  mass  diffusion  can 
play  a  non-negligible  role  in  influencing  the  sta¬ 
bility  properties  of  such  systems;  see.  by  analogy. 
Garcia-Ybarra  ct  al.  [2.^]. 

Eq.  (74)  extends  our  preliminary  results. 
Castillo  et  al.  [14,15],  to  arbitrary  Biot  numbers  at 
both  ends  and  also  incorporates  the  influence  of 
the  latent  heat.  We  proceed  now  to  analyze  the 
most  important  features  of  the  dispersion  rela- 
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lion.  cq.  (74).  The  growth  rate  of  a  planar  distur¬ 
bance  (zero  wavenumber)  is 

ff{Bi,,  *  I).  k  =  0} 

Pa  l  ^  ) 


Bi  + 


A 

1  -1-  Bi  -1-  — 

y' 

Bi  -)-  -  I 

1  -i-^Bij 

Bid 

P 

1 

(76a) 


Note,  however,  that  this  growth  rate  vanishes  in 
the  limit  of  a  perfect  insulating  support,  when  a 
quadratic  dependence  for  .small  wavenumbers  is 
obtained 

crfBij  =  0.  k  -»  0} 


right  hand  side  numerator  is  missing  and  in  the 
denominator,  0\  must  be  replaced  by  0\  + 
A/G,). 

Eq.  (77)  shows  that  when  (p  -  +  p^)> 

0,  the  amplification  factor  <t  increases  linearly 
with  k  and  the  smaller  wavelength  disturbances 
arc  the  most  dangerous  to  the  planar  growth 
because  they  amplify  faster.  However,  in  this 
analysis,  the  value  of  the  disturbance  wavenum¬ 
ber  has  been  assumed  to  be  of  order  unity  in  the 
Schmidt  number  expansion.  This  assumption  fails 
when  the  wavenumber  takes  on  large  values  and 
a  different  scaling  law  is  required  to  improve  the 
result  shown  in  eq.  (74).  This  new  scaling  is 
presented  in  the  next  section. 


7.  Study  of  large  wavenumber  disturbances 
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Pd  \  ^ 


Bi  -(-  p/e. 
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y'i  P 
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In  the  case  of  very  large  wavenumber  distur¬ 
bances  (i.e.  when  k  =  Sc  a.  with  a  =  0(Sc")).  the 
stability  analysis  performed  in  section  6  is  no 
longer  valid.  The  variations  with  x  become  more 
important  and  a  new  scaling  is  needed.  There¬ 
fore.  we  define  the  new  rescaled  coordinate  if 
and  time  t'  via: 


On  the  other  hand,  in  the  opposite  limit  of  very- 
large  wavenumbers,  cq.  (74)  leads  to  a  linear 
behavior 


which  coincides  (for  Pj  =  p)  with  the  Huric  [16] 
dispersion  relation  in  the  limit  of  very  high  Soret 
effect  (once  the  several  misprints  in  Hurle's  equa¬ 
tions  are  corrected:  using  his  notation.  NGy/k 
should  be  replaced  by  A/G,  in  his  eq.  (11),  in  his 
eq.  (2.1)  a  minus  sign  in  the  first  contribution  to 
(■>*  is  forgotten,  and  in  his  dispersion  relation,  eq. 
(24),  a  in  the  left-hand  side  of  the  numerator 
should  be  in  the  denominator,  the  bracket  in  the 


tp  =  Sc  X,  t’ s  Sc  t.  (78) 

It  can  be  shown  that  for  a  given  order  of  the 
mass  fraction  disturbance  amplitude,  the  distur¬ 
bances  in  the  front  shape  and  in  the  temperature 
fields  are  Sc"'  times  smaller.  The  mass  fraction 
field  in  the  vicinity  of  the  interface  is  very  sensi¬ 
tive  to  front  corrugations.  Then,  different  expan¬ 
sions  are  needed  for  the  temperature  fields  in  the 
inner  and  outer  regions.  Therefore  we  take: 

1 

77  =  — 77,  (79a) 

Sc 

^^ou.  =  ()<-•<  1-77. 

(79b) 

=  +  ()<^<x.  (79c) 
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-  f>  -  rj  <  z  <{). 


«d.,n  =  +  »U..„  )  -  *  <  ^ 


n  »  O'*  +  r 

^*oul  '‘oul  *  i)u« 


IA„  =  /C  +  /'.n 


()<r  <  I  -  Tj, 

(79f) 

0<f  <x,  (79g) 


Eqs.  (79)  define  the  tilde-quantities.  Then  the 
operators  T  and  T-  in  the  outer  regions  beeome 

/  (<  'd  (t  \ 

'■=  Sc— (KO) 
^  dif  dz  dz  f 

d-  'd-fj  d  drj  d- 

r-  =  Sc--^  -  Sc—r  —  -  2Sc—  r-7- 

dif'  dif'  dtf  dz  dif 


+  1  + 


dri  \']  0- 


dif  j  ^2  - 


and  in  the  inner  regions 

Id  dri  d  d 

r=  Sc—  -  Sc— —  .Sc— 

y  dtf  dtp  d(  d( 


,  ,  f*  ,  'dii  d- 

V-  =  Sc--^  -  Sc’—;  —  -  2ScM 

dif~  dip~  dif  d^  dip 

Then  to  lowest  order  in  the  Schmidt  number 
expansion,  the  linearized  equations  for  the  distur- 
banees,  from  eqs.  (37)  to  (39)  with  (80)  and  (81 )  in 
the  outer  region,  beeome 

-^-—^=.0,  (84) 

Oif  dip 


-0 

dip^  dip^-^'^ 


whereas,  taking  into  aceount  (82)  and  (83),  in  the 
inner  regions,  we  obtain: 

a’«,„  d'r,  d-e,„ 

irT-ir3^  +  i;7r  =  “-  («7) 

Oif  Pif 

dip-  dip'^'^ 

^  ^  ap  di\ 

d('  dip-  dip-  d£  0„  df 

«  ^1.  diC 

+  (89) 

6)„  di 

The  boundary  conditions  at  the  interface,  from 
eqs.  (53),  (54).  (56)  and  (57)  are 


«.n  =  eu..n 


af  p  dr  A  d^ 

^  +  —r 


at  ^  =  0; 


at  ^  =  0;(91) 


dT  pj  \  a^ 


+  —O' 


- 

=  Ars®.n-  rs-^ 

dip 


at  ^  =  0;(92) 


at  f  =  0;(93) 


with  Xs  -  ^f/Sc  and  =  ySc,  taken  to  be  of 
order  unity. 

As  in  the  ease  studied  in  section  6.  the  system 
of  linear  equations  can  be  analyzed  in  a  Fourier 
series,  by  considering  any  disturbance  in  the  form 

/{v*.  Z,  t’}  =  F{Z}  exp(io^ -t-i^r'),  (94) 

where  F  stands  for  rj,  and  F,  and  Z  for  2 

(for  the  outer  quantities)  or  (  (for  the  inner 
variables). 

Then  the  solution  of  eqs.  (84).  (85),  (87)  and 
(88)  satisfying  boundary  conditions  (90)  and  (91) 
and  the  matching  principle  is 


0. 


(95) 
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IJ'  1  Pj  A  /3-^j 

®,n  >  -*■  - - - JT- 

\P  a  p  Aj  P 


X(1  +A/AJ'  c 


a  solution  different  from  the  trivial  one  when  the 
following  eompatibility  condition  is  fulfilled: 


P-Pa 

a - r 


\P^  Alylp^  p-p. 


p  aAP  a  P  ^  P 


x(l  +A/Aj)''  e“'  . 


Also,  from  eq.  (K6).  we  obtain 


I  \  “  I 

,rs«.  ■>  jll 


(w) 

Note  that,  to  leading  order,  the  boundary  condi¬ 
tions  (58)  to  (6(1)  are  automatically  satisfied,  irre¬ 
spective  of  the  Biot  numbers.  Thus,  the  outer 
regions  act  as  a  thermal  bath  for  the  inner  region 
and  the  boundary  conditions  far  from  the  inter¬ 
face  play  no  role. 

The  solution  to  eq.  (89)  which  vanishes  when  | 
goes  to  infinity  (to  match  with  the  outer  solution 
(99)1  is 

I  /A?  \ 

r.- 


/  i'  Pj  A  P-Pj 

X  1  -t- - 

\  ^  fl  p  Aj  P 


X(1  +A/AJ'‘  e 


I  a0\- 


The  constants  A  and  f)  in  eq.  (1(X))  should 
satisfy  two  different  conditions,  coming  from  eqs. 
(92)  and  (93).  The  resulting  algebraic  system  has 


^  fl  /3  Aj  /?,.  (  a  (9„  I  I 


In  the  limit  of  small  values  of  a.  this  expression 
reduces  to 


y 

1  +  - 


Comparing  this  result  with  the  asymptotic  limit 
of  eq.  (74)  for  very  large  wavenumbers  (which  is 
given  by  eq  (77)).  clearly  shows  the  same  trends. 
Therefore,  both  dispersion  relations,  eq.  (74)  and 
eq.  (102).  match  perfectly  in  the  intermediate 
wavenumber  range  and  a  composite  relation  can 
be  written  down  as  usual,  just  by  adding  the 
right-hand  side  of  both  dispersion  relations  and 
subtracting  the  common  part  given  by  the  right- 
hand  side  of  eq.  (103). 

In  order  to  present  the  results  in  a  simple  way. 
we  define  the  saturation  parameter 


^2^ 


0 1  4e„n;i- 


I.  M 


(104a) 
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with 


then,  cq.  (H)2)  becomes 


H,,  jdn:;, 

de  d( 

Note  that  ./  can  also  be  written 


(KWb)  v«  _  +  (/l  -t-fl*-  -  l)(y'-a*‘) 

\/l  +a*-  -  1  >7/3 

A - - - - - / 

a*  2+y'/p 


(lOS) 


./  = 


A 

r  ^  j 

1  4/27-), 1 

Ys 

A 

1  .  — 
Aj 

1  a/3&{ 

^  r 

(KWe) 


Both  parameters  5  and  ./  have  the  same  sign 
and  they  measure  the  degree  of  saturation  in  the 
vicinity  of  the  growing  deposit.  Note  that  when  ./ 
is  negative  the  Brownian  sublayer  is  undersatu- 
rated  ahead  of  the  moving  fluid-solid  interface 
whereas  positive  values  of  .!/  correspond  to  con¬ 
stitutional  supersaturation. 

On  the  other  hand,  the  parameter 


0  ys<’\.afi(]  +  A/Aj) 


(105) 


provides  a  measure  of  the  thcrmophorctically  in¬ 
duced  instability.  When  .V  is  positive,  planar 
crystal  growth  is  unstable,  as  indicated  in  section 
6. 

In  terms  of  the  parameters  .!/’  and  .'<J,  the 
dispersion  relation  expressed  by  eq.  (102)  can  be 
rewritten  in  a  simplified  manner  by  using  the  new 
normalized  amplification  factor  J*  and 
wavenumber  a* 


Pd  20„  , 


* 


p  afi  ^  ^0/3  Sc  ’ 


and  defining  the  parameter 


,1.-11.- 


(9,a^  rs(l+A/AJ  ' 


(106) 


(107) 


Within  the  range  of  validity  of  the  present 
analysis,  the  sign  of  1*  in  eq.  (108)  is  given  by 
the  numerator  of  the  last  fraction  (occasionally, 
the  denominator  may  vanish  or  become  negative, 
but  this  happens  either  for  unrealistically  large 
values  of  the  surface  tension  parameter  or  for 
very  large  values  of  the  wavenumber  and,  in  both 
cases,  the  expansion  in  the  Schmidt  number  we 
have  performed  breaks  down).  When  the  numer¬ 
ator  is  positive  for  some  wavenumber  a*,  planar 
growth  is  unstable.  Therefore  the  stability  of  the 
front  is  governed  by  only  two  parameters,  the 
thermophoretic  parameter  defined  in  eq.  (lO.S) 
and  the  saturation  parameter  v  defined  in  eq. 
(104).  These  parameters  account  for  the  two  mor¬ 
phological  instabilities  which  may  arise  during 
planar  growth  of  a  solid  deposit  by  thermal  diffu¬ 
sion  mass  transport. 

When  is  positive,  the  amplification  factor 
is  positive,  at  least  for  the  small  values  of  o*. 
This  is  the  thermophoretically-induccd  morpho¬ 
logical  instability  pointed  out  in  the  previous  sec¬ 
tion.  Besides,  eq.  (108)  shows  that  for  very  large 
wavenumbers,  a*,  the  stabilizing  effect  of  surface 
energy,  measured  by  the  parameter  and  repre¬ 
sented  by  tne  last  term  in  -o*'.  is  dominant,  and 
1*  becomes  negative.  Therefore  large  wavenum¬ 
ber  disturbances  are  always  damped  out  by  their 
excess  in  surface  energy. 

On  the  other  hand,  whenever  >  is  positive 
there  exists  a  constitutional  supersaturation  in 
the  Brownian  sublayer.  As  it  was  schematically 
pointed  out  in  the  Introduction  and  shown  by 
Huric  [25]  and  Sekerka  [26].  this  is  a  potentially 
unstable  situation.  This  is  the  constitutional  su¬ 
persaturation-induced  morphological  instability. 
However,  when  .y  is  negative,  interface  protu¬ 
berances  have  more  difficulty  in  developing  into 
the  fluid  region  because  it  is  a  thermodynamically 
stable  phase. 
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Fip,  1.  Stuhllily  diagram  of  Ihcrmophorclically  growing  de¬ 
posits.  Planar  growth  is  unstable  in  regions  I.  II  and  III  and 
stable  in  regains  IV  and  V.  In  each  region  a  schematic 
representation  of  the  amplification  factor  2  *.  as  a  function  of 
the  disturbance  wavenumber  u*.  is  presented.  The  solid  line 
separating  regions  III  and  IV  corresponds  t<i  the  critical 
values  given  by  egs.  (11(1)  whereas  the  dashed  line  is  the 
asymptotic  result  given  by  eq.  (III). 


The  main  features  of  the  dispersion  relation 
( lOS)  are  depicted  in  fig.  1.  There  exist  five  differ¬ 
ent  regions  in  the  parameter  space  labeled  I  to  V 
in  fig.  1.  The  behavior  of  the  linear  amplification 
factor.  1*,  as  a  function  of  the  disturbance 
wavenumber,  a*,  is  schematically  indicated  in 
each  region. 

In  region  I.  .'i  is  positive  and  the  planar  growth 
is  unstable  due  to  the  morphological  instability 
indicated  in  section  b,  but  is  negative,  there¬ 
fore  the  fluid  ahead  of  the  interface  is  undersatu¬ 
rated.  For  large  values  of  a*,  the  stabilizing 
influence  of  the  undersaturation  overcomes  the 
destabilizing  effect  of  thermal  diffusion  within 
the  diffusive  boundary  layer.  The  right-hand  side 
term  of  eq.  (108)  becomes  negative  and  large 
wavenumber  disturbances  arc  damped  out. 
Therefore,  only  disturbances  with  wavelengths 
larger  than  the  diffusive  boundary  layer  thickness 
arc  expected  to  develop.  Thus,  the  undersatura- 
tion  modifies  the  result  given  by  eq.  (77),  which 
predicted  a  monotonic  increase  of  the  amplifica¬ 


tion  factor  with  the  disturbance  wavenumber 
when  Brownian  diffusion  and  surface  energy  ef¬ 
fects  were  not  taken  into  account.  Only  distur¬ 
bances  with  wavenumbers  ranging  from  a*  =  0  to 
the  value  at  which  {a*}  vanishes  again,  arc 
linearly  unstable. 

In  region  11,  both  parameters  ^  and  •/'  arc 
positive.  The  planar  front  is  already  unstable  for 
the  thermal  diffusion-driven  morphological  insta¬ 
bility  indicated  in  the  previous  section,  and  a 
positive  value  of  results  in  even  higher  ampli¬ 
fication  factors  for  intermediate  wavenumber  dis¬ 
turbances.  This  is  the  most  unstable  situation. 
Both  effects  cooperate  to  render  the  planar 
growth  unstable  until  surface  energy  effects  take 
over  and  stabilize  the  very  large  wavenumbers. 

In  region  111,  is  a  negative  quantity  and  the 
planar  front  is  stable  for  thermal  diffusion  insta¬ 
bility,  but  a  high  enough  supersaturation  (that  is. 
a  positive  and  large  value  of  y)  render  the  value 
of  i  *  in  eq.  (108)  positive  for  some  intermediate 
values  of  the  wavenumber,  a*.  Then  a  .supersatu¬ 
ration  induced  morphological  instability  appears. 
The  small  wavenumber  disturbances  arc  stabi¬ 
lized  by  latent  heat  release  but  at  an  intermediate 
range  of  wavenumbers,  supersaturation  over¬ 
comes  these  stabilizing  effects  and  makes  unsta¬ 
ble  the  planar  growth.  Again,  surface  energy  ef¬ 
fects  stabilize  the  very  large  wavenumbers. 

In  region  IV.  is  negative  and  is  positive 
but  is  not  large  enough  to  overcome  the  stabiliz¬ 
ing  effects  of  both  latent  heat  release  and  surface 
energy  excess.  Planar  growth  is  stable. 

Finally,  in  region  V.  both  and  .y  are  nega¬ 
tive  and  the  planar  interface  during  the  growth  of 
a  solid  deposit  by  thermophoretic  transport  is 
linearly  stable.  Disturbances  of  any  wavenumber 
are  damped  out. 

All  these  mentioned  features  are  clearly 
pointed  out  when  the  dispersiem  relation  is  ex¬ 
panded.  around  the  origin  in  the  parameter  space, 
for  small  values  of  the  wavenumber  as 

=  +  (109) 

this  simple  expression  incorporates  the  whole 
qualitative  description  of  the  morphological  insta- 
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bility  in  the  five  distinctive  regions  of  the  parame¬ 
ter  space. 

Over  the  line  separating  region  111  from  region 
IV.  1*  starts  at  a  zero  value  at  fl*=0  and 
decreases  as  a*  increa.ses  due  to  the  stabilizing 
role  of  latent  heat  release.  Further  increases  of 
a*  tend  to  increase  the  value  of  2*  due  to  the 
constitutional  supersaturation,  but  1*  reaches  a 
maximum  value  equal  to  zero  and  comes  down 
back  again  due  to  the  surface  energy  stabilizing 
influence.  Therefore,  the  locus  of  this  critical  line 
may  be  easily  obtained  by  noting  that  2*  van¬ 
ishes  when  it  is  also  a  maximum  respect  to  a*. 
Imposing  these  two  conditions  (2*  =  0  and 
d2*/dfl*  =0),  wc  obtain  the  critical  values  for 
the  saturation  parameter,  critical 

wavenumber,  a*,,,  as  functions  of  the  thermal 
parameter  3 .  It  results  that 


(  ^cril  ) 


I +2\1 +«,,)- 


( 1 10a) 


Fig.  2.  Normalized  critical  wavenumber  u*.  as  a  function  of 
the  thermophoretic  parameter  H'.  Solid  line  represents  the 
exact  result  from  eq.  (IlOb).  whereas  dashed  line  is  the 
asymptotic  result  given  by  eq.  ( 1 1 1 ). 


-“iTU 


I  +«..)■-  Vl+(«cn.)' 


(110b) 


Thus,  given  a  negative  value  of  .■!/,  eq.  (110b) 
provides  the  critical  wavenumber.  and  eq. 
(1 10a)  the  critical  value  of  the  saturation  parame¬ 
ter.  which  corresponds  to  the  separation 

line  between  region  111  and  region  IV  in  fig.  1. 
For  positive  values  of  or  values  of  .y’  larger 
than  the  critical  value  (which  is  a  function  of  .'■5 
only),  planar  growth  is  unstable. 

The  behavior  of  the  critical  values,  for  small 
negative  values  of  can  easily  be  obtained  from 
eqs.  (109)  and  (110),  resulting  in 

=  3(-.¥^)'/\  = 

when  0  <  -y^  •«  1 . 


Also  the  line  corresponding  to  =  0.  and  neg¬ 
ative  values  of  the  saturation  parameter  ^  (i.e. 
y  =  0  and  J'  <  0),  is  a  critical  line  because  it 
separates  stable  and  unstable  regions  of  planar 
growth.  Over  this  line,  the  critical  wavenumber, 
a*i,,  is  always  equal  to  zero. 

When  y  is  negative,  the  critical  value  of  the 
saturation  parameter,  •nn.  .  provided  by  eq.  (109a) 
is  a  monotonically  increasing  function  of  {-S’). 
Therefore,  the  range  of  experimental  operation 
conditions  for  stable  planar  growth  is  enlarged 
beyond  the  limit  growth  velocity  given  by  the 
oversimplified  criterion,  eq.  (3),  due  to  the  stabi¬ 
lizing  influence  of  both  latent  heat  release  and 
interfacial  energy.  Thus,  we  can  conclude  that  for 
negative  values  of  y,  the  maximum  growth  veloc¬ 
ity  for  stable  planar  growth  is 


This  asymptotic  results  is  plotted  in  fig.  1  as  a 
dashed  line. 

The  critical  wavenumber,  is  plotted  in  fig. 
2  as  a  function  of  the  thermal  parameter  y.  The 
solid  line  corresponds  to  the  exact  result  given  by 
eq.  (110)  and  the  dashed  line  gives  the  asymptotic 
limit  for  small  values  of  the  thermal  parameter 
provided  by  eq.  ( 1 1 1 ). 
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where  the  critical  value  .S,,,,,  may  be  obtained 
from  the  definition  (l(l9a)  once  is  known. 
As  is  always  positive  (when  the  planar  growth 
is  thcrmophorctically  stable;  i.e.,  for  negative  val¬ 
ues  of  ^  )  the  limit  velocity  is  larger  than  the 
given  by  the  simplified  criterion.  However,  for 
ptrsitivc  values  of  the  thcrmophoretic  parameter 
planar  growth  is  unstable  for  any  given  veloc¬ 
ity. 


8.  Discussion  of  results  and  conclusions 

We  have  performed  a  linear  stability  analysis 
of  a  solid  deposit  with  a  planar  interface  growing 
from  a  fluid  mixture.  The  depositing  species  is 
considered  to  be  dilute  in  the  carrier  fluid  and 
the  growth  limited  by  thermal  diffusive  (Soret) 
transport,  with  Brownian  diffusion  playing  a  mi¬ 
nor  role  (restricted  to  a  boundary  (sub-)layer  in 
the  immediate  vicinity  of  the  solid-liquid  inter¬ 
face)  and  without  considering  interfacial  kinetic 
limitations.  Due  to  the  diluteness  assumption,  the 
Peclet  number  associated  with  the  advancing  so¬ 
lidification  front  is  quite  small  and  the  unsteady 
term  as  well  as  the  convective  terms  have  been 
neglected  in  the  governing  equations. 

The  resulting  dispersion  relation  between  the 
disturbance  amplification  rate  cr  and  the  distur¬ 
bance  wavenumber  k  is  given  by  eq.  (74),  which 
describes  a  new  type  of  morphological  instability. 
The  driving  mechanism  is  the  thermal  (Soret) 
diffusion  phenomenon.  Sufficient  conditions  for 
the  appearance  of  the  instability  are  explicitly 
given  by  eq.  (75c)  and  eq.  (75d).  Planar  growth  is 
unstable  either  when  the  thermal  diffusion  factor 
of  the  depositing  species  is  positive  (as  is  the 
case  for  very  massive  molecules  in  a  light  carrier 
gas)  and  the  latent  heat  release  is  very  small,  or 
when  Oj  is  negative  and  the  upper  boundary  is  a 
good  thermal  insulator.  To  be  consistent  with  the 
assumption  of  a  comparatively  small  diffusive 
(Fickian)  mass  transport,  Oj  should  be  consid¬ 
ered  positive,  Rosner  [27].  This  situation  corre¬ 
sponds  to  large  depositing  molecules  dilute  in  a 
light  carrier  gas.  Such  particles  drift  by  thermal 
diffusion  towards  cooler  regions,  that  is,  the  de¬ 
posit  should  then  be  colder  than  the  fluid  mix¬ 


ture.  In  that  case,  eq.  (74)  shows  that,  subject  to 
our  present  model  idealizations,  when  the  latent 
heat  release  is  very  small  (and  the  thermal  con¬ 
ductivity  of  the  solid  is  larger  than  the  fluid 
thermal  conductivity,  as  it  usually  happens),  there 
is  no  possibility  of  stable  planar  growth  -  i.e., 
disturbances  of  any  wavenumber  will  be  magni¬ 
fied,  with  the  amplification  factor  being  a  mono- 
tonically  increasing  function  of  wavenumber. 
However,  the  analysis  leading  to  eq.  (74)  breaks 
down  when  the  disturbance  wavelength  is  of  the 
order  of  the  thickness  of  the  diffusive  boundary 
layer  [that  is,  when  k  =  0(Sc)].  Then,  transverse 
diffusive  fluxes  should  be  taken  into  account  and 
eq.  (102)  provides  the  corresponding  dispersion 
relation.  The  implemented  dispersion  relation  in¬ 
corporates  a  second  type  of  morphological  insta¬ 
bility  induced  by  constitutional  supersaturation. 
The  two  kinds  of  instabilities  may  compete  or 
cooperate  in  this  interval  of  large  wavenumber 
disturbances.  Thus,  when  the  diffusion  sublayer  is 
stable  by  itself  (that  is,  when  there  is  no  constitu¬ 
tional  supersaturation  near  the  interface),  this 
stabilization  mechanism  may  provide  a  cutoff  to 
the  instability  induced  by  thermophoresis.  On  the 
other  hand,  when  constitutional  supersaturation 
is  present,  the  band  of  unstable  wavenumber  is 
enlarged  up  to  the  scale  where  interfacial  energy 
phenomena  take  over  and  the  smaller  wave¬ 
lengths  are  damped  out.  When  the  latent  heat 
release  is  able  to  stabilize  the  small  wavenumbers 
disturbances  but  supersaturation  exists,  the  desta¬ 
bilizing  effect  of  supersaturation  may  render  un¬ 
stable  the  planar  growth  for  some  wavenumber  of 
the  order  of  the  Schmidt  number.  Under  these 
circumstances,  a  front  corrugated  on  very  small 
length  scales  is  expected. 

Debe  and  Poirier  [28]  have  recently  reported 
some  experimental  results  on  the  physical  vapor 
deposition  of  dilute  copper  phtalocyanine  under 
microgravity  and  ground-based  conditions.  Al¬ 
though  in  the  latter  case  buoyancy  currents  *' 
appearing  in  the  bulk  may  alter  the  resulting 
crystallization  pattern,  they  observed  in  both  cases 

Gas  “creep"  along  the  non-isothermal  solid  sidewalls  is 

also  expected  to  drive  an  appreciable  convection,  even 

under  micrograviiy  conditions  (Rosner  129]). 
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the  formation  of  a  very  corrugated  front,  leading 
finally  to  the  development  of  non-linear  struc¬ 
tures  of  finger-shaped  form,  with  an  average  size 
of  the  order  of  10-^  A.  Copper  phtalocyanine 
molecules  are  very  large  and  in  the  usual  gas 
carriers  they  exhibit  a  low  Fickian  mobility  (with 
Schmidt  numbers  about  10)  and  a  large  thermal 
diffusion  effect  (with  thermal  diffusion  factors 
about  7).  Thus,  the  solid  growth  was  actually 
controlled  by  thermal  diffusion.  Under  their  ex¬ 
perimental  conditions,  constitutional  supersatura¬ 
tion  was  actually  present  (see,  e.g.  Rosner  and 
Keyes  [30]).  Moreover,  the  copper  phthalocyanine 
was  very  dilute  in  the  gas  mixture  and,  then,  the 
latent  heat  release  rate  was  very  small  and  unable 
to  stabilize  the  small  wavenumber  disturbances. 
Therefore,  planar  growth  was  unstable  against 
both  types  of  morphological  instabilities.  To  avoid 
or  reduce  the  instability  one  might  use  a  carrier 
of  larger  molecular  size  (comparable  to  that  of 
the  depositing  species)  although  the  “price”  paid 
would  be  a  notable  decrease  in  the  molecular 
diffusion-limited  growth  rate.  In  any  case,  our 
approach  predicts  an  important  dependence  of 
surface  morphology  on  the  nature  of  the  carrier 
gas,  which  may  have  interesting  practical  implica¬ 
tions. 

On  the  other  hand,  the  experiments  conducted 
by  Reed  and  LaFleur  [7]  on  the  growth  of  iodine 
solids  from  iodine-helium  mixtures  were  carried 
out  for  heavy  mass-loaded  mbctures  and,  in  spite 
of  the  molar  dilution,  the  mass  diluteness  as¬ 
sumption  can  not  be  applied  and  the  Stefan  flow 
must  play  a  non-negligible  role.  Furthermore,  the 
Schmidt  number  was  of  the  order  of  the  unity 
and  Fickian  diffusion  transport  was  as  important 
as  thermal  diffusion  mass  transfer.  In  any  case, 
the  planar  growth  was  found  to  be  initially  stable, 
becoming  unstable  only  for  sufficiently  large 
growth  velocities.  The  experimentally  observed 
critical  growth  velocity,  as  indicated  by  Rosen- 
berger  et  al.  [8],  was  one  order  of  magnitude 
larger  than  the  value  predicted  for  the  oversim¬ 
plified  criterion,  eq.  (3).  Our  result  (eq.  (112)) 
explains  qualitatively  this  discrepance  because  a.^ 
as  well  as  are  positive  and  a  larger  growth 
velocity  may  be  achieved  with  a  planar  interface 
still  remaining  stable. 


It  is  clear  that  the  additional  mechanisms  of 
Fickian  diffusion  and  Stefan  flow  cannot  be  ne¬ 
glected  in  many  practical  situations,  and  the  pre¬ 
sent  analysis  should  be  extended  to  account  for 
mixtures  in  which  both  diffusive  mass  fluxes 
(Fickian  diffusion  and  thermal  (Soret)  diffusion) 
as  well  as  the  deposition-induced  convective  mass 
transport  have  comparable  importance  through¬ 
out  the  domain.  It  would  be  also  interesting  to 
relate  this  class  of  continuum  models  to  discrete 
computer  simulations  of  deposit  formation  (see, 
e.g.,  Tassopoulos  et  al.  [31]  and  Rosner  et  al.  [32]) 
which  predict  deposit  microstructural  properties 
(surface  roughness,  porosity,  etc.).  Our  further 
studies  will  be  directed  along  these  lines. 
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Nomenclature 


A,  B,  C,  E 
a 

Bi 

fd 

D 


Integration  constants 
Rescaled  disturbance  wavenum¬ 
ber,  a  =  k/Sc 

Biot  numbers  introduced  in  eqs. 
(58)  and  (59) 

Deposit  specific  heat 

Fluid  specific  heat  at  constant 

pressure 

Binary  diffusion  coefficient 
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»■ 

Thermophoretic  parameter  de¬ 

P 

Dimensionless  fluid  temperature 

fined  in  eq.  (105) 

gradient 

AH 

Molar  heat  of  solid  deposit  sub¬ 

r 

Mass  fraction  disturbance 

limation 

y 

Dimensionless  interfacial  energy 

K 

Dimensionless  local  surface  cur¬ 

parameter,  y  = 

vature 

(p^RTLu,„t\  eq.  (32b) 

k 

Dimensionless  disturbance 

7s 

Stretched  y,  y^  =  ySc 

wave-number 

s 

Thickness  ratio,  8  s  Lj/L 

Dimensionless  latent  heat  pa¬ 

Interfacial  deformation,  eq.  (10) 

rameter  4  A/  Pe/(Af,c^T„f) 

V 

Dimensionless  interfacial  defor¬ 

L 

Fluid  layer  thickness 

mation,  TJ  s  (/L 

La 

Deposit  thickness 

0 

Dimensionless  temt)erature 

1 

/  =  L-(-L, 

field,  0  s  T/T„, 

M 

Molar  mass 

e 

Temperature  disturbance 

N 

Normalization  factor  defined  by 

A 

Thermal  conductivity 

eq.  (11) 

A 

Parameter  defined  by  eq.  (107) 

n 

Unit  vector  normal  to  the 

M 

Dynamic  viscosity 

fluid-solid  interface,  eq.  (11) 

V 

Kinematic  viscosity,  v  =  p/p 

Pe 

Fluid  thermal  Peclet  number,  Pe 

Inner  stretched  coordinate,  ^  = 

s  y„,L/lA/(pc^)] 

Sc  z 

Pr 

Prandtl  number. 

n 

Dimensionless  pressure  field. 

Pr  =  v/[A/(pCp)] 

n  =  pL/pV„f 

P 

I-ocal  pressure 

p 

Density 

R 

Universal  gas  constant 

I 

Normalized  disturbance  amplifi¬ 

r 

Exponential  factor  defined  by 

cation  factor,  1  =  o-/Sc 

eq.  (101) 

<T 

Disturbance  amplification  factor 

Interface  position  vector,  eq.  (10) 

try 

Interfacial  tension  of  the  fluid- 

5 

Measure  of  supersaturation  de¬ 

solid  interface 

fined  by  eq. (104b) 

t' 

Stretched  dimensionless  time,  r' 

Supersaturation  parameter  de¬ 

=  T  Sc 

fined  by  eq.  (104a) 

T 

Dimensionless  time,  t  = 

Sc 

Schmidt  number.  Sc  s  p/O 

<t> 

Stretched  horizontal  coordinate. 

T 

Local  temperature 

<^  =  Sc  x 

t 

Time 

X 

X  =  (dfA'Vd©)  at  ©  =  00 

t 

Unit  vector  tangent  to  the 

Xs 

Stretched  Xs  =  a/Sc 

fluid-solid  interface,  eq.  (11) 

n 

Normalized  mass  fraction,  {i  = 

V 

Local  velocity 

V 

Dimensionless  velocity,  r  s 

w 

Depositing  species  mass  frac¬ 

tion,  CD  =  p|/p 

X,Y 

Space  coordinates 

X,  y 

Dimensionless  coordinates 

Subscripts 

z 

Dimensionless  vertical  distance 
from  the  interface, 
z  S  y  -  rjlx,  t) 

Critical  values,  see  eqs.  (110) 

crit 

d 

Depxjsit 

a 

Dimensionless  thermophoretic 

e 

Values  at  y  =  L 

diffusivity,  a  =  OjD/v 

in 

Inner  solution 

«T 

Thermal  diffusion  factor 

out 

Outer  solution 
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ref  Reference  values  used  to  nondi- 

mensionalize  the  variables 
s  Interface 

0  Refers  to  the  value  at  the  inter¬ 

face  during  stationary  growth 

1  Refers  to  the  deposit  forming 
fluid  component 

2  Refers  to  the  inert  (carrier)  fluid 
component 

8  Values  at  y=  -Lj 

Sititerscripls 

eq  Equilibrium  mass  fraction  val¬ 

ues.  eq.  (19) 

St  Stationary  state,  section  5 

*  Normalized  values  defined  by 

eqs.  (106) 

Miscellaneous 

Disturbance  z-dependence  (sec 
eq.  (61)) 

Refers  to  disturbances  defined 
by  eq.  (79) 

{  )  Functional  parentheses 

r.  r*  Gradient  and  Laplacian  opera¬ 

tors,  respectively 
nst  Non-steady  terms 
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Heat  transfer  (associated  with  Fourier  diffusion,  photon  absorption  and/or  emission,  surface  chemical 
reaction,  and/or  evaporation)  can  lead  to  dramatic  changes  in  the  motion  of  small  particles  in  gases, 
modifying  coagulation  rates,  and/or  capture  rates.  Cases  discussed  in  this  review  of  recent  studies 
include  (a)  fhermophoretic  drift  of  spherical  particles  in  a  gas  through  which  energy  is  diffusing, 
(b)  photophoretic  d^  of  radiation-abrorbing  suspended  spheres,  (c)  phoresis  and  apparent  Brownian 
diffusivity  of  asymmetric  particles,  including  spheres  with  nonuniform  accommo^tion  coefficient 
and  unequal  primary  partide  size  aggregates,  (d)  scavenging  of  fine  particles  by  evaporating  droplets, 
and  (e)  coagulation  rate  of  unequal  temperature  spherical  particles.  These  transport  phenomena, 
often  associated  with  the  localized  “overh^ting”  or  ‘Hmdercooling”  of  individual  partides,  have  largely 
been  ignored  by  applied  sdentists  on  the  (usu^y  implidt)  assumption  that  the  more  familiar  effects 
of  ordinary  Brownian  diffusion  would  dominate  in  near-isothermal  laminar  gas  (sub-)regions. 
However,  we  show  that  even  in  regions  free  of  macroscopic  (bulk)  temperature  gradients,  these 
phenomena  can  dominate  ordinary  Brownian  transport  by  orders  of  magnitude,  with  potentially 
important  engineering  consequences.  General  expressions  for  the  relevant  particle  drift  velocities, 
coagulation  rate  constants,  or  effective  particle  diffusivities  are  used  to  carry  out  illustrative  cal- 
oilations  for  multiphase  systems  of  current  and  future  engineering  interest,  including  fine  particle 
synthesis/ materials  processing,  and  combustion  gas  cleaning. 


1.  Introduction,  Motivation,  and  Goals 

Technologies  that  exploit  or  inadvertently  generate 
particle-laden  gases  are  often  also  associated  with  large 
“convective”  and  radiative  energy  fluxes,  as  in  the  cases 
of  optical  waveguide  glass  deposition  (see,  e.g.,  Roener  and 
Park,  1988)  or  pulverized-coal-flred  electric  power  stations 
(see,  e^.,  Flagan  and  Friedlander,  1978).  Yet,  the  science 
of  “aerosols”  had  its  origins  in  the  understanding/control 
of  environmental  problems  which  were  relatively  unin¬ 
fluenced  by  simultaneous  energy  transfer  (see,  e.g., 
Whytlaw-Gray  and  Patterson,  1932;  Fuchs,  19M;  Fried¬ 
lander,  1977;  Hidy,  1984;  Flagan  and  Seinfeld,  19S^).  For 
example,  flne  particles  suspended  in  isothermal  “fluids” 
undergo  Brownian  motion  due  to  molecular/particle  im¬ 
pacts  (Einstein,  1926)  and  this  random  walk  motion  is,  in 
turn,  associated  with  the  now  well-known  phenomena  of 
particle-particle  “coagulation”  and  particle  Brownian 
diffusion  to  macroscopic  solid  surfaces  (immersed  or  con¬ 
tainment).  The  simultaneous  presence  of  signiflcant  en¬ 
ergy  fluxes  (e.g.,  Fourier  conduction  or  radiative  transport 
through  the  carrier  gas)  raises  the  important  question  of 
if  and  how  the  famiUar  laws  of  particle  “diffusive”  trans¬ 
port  and  coagulation  are  modifled  from  those  normally 
associated  with  isothermal  gases  in  the  absence  of  a  net 
energy  flux. 

Our  recent  research  on  combustion  gases  containing 
suspended  spherical  or  nonspherical  particles,  briefly 
summarized  in  this  paper,  has  provided  many  examples 
in  which  the  nuus-tran^er  effects  of  simultaneous  energy 
transfer  are  quite  dramatic,  because  of  the  onset  of  “new” 
mechanisms  of  particle  transport  in  the  presence  of  gas 
temperature  gradients  and/or  radiative  fluxes.  Actually, 
some  of  these  mechanisms  (conventional  “thermophoresis” 
(section  2.1)  and  “photophoresis”  (section  2.2))  are  dis¬ 
cussed  in  cliesical  aerosol  science  texts,  but,  until  relatively 
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recently,  were  not  thought  to  be  of  engineering  importance. 
Others  (associated  with  asymmetrical  particles  in  the 
presence  of  energy  transfer  (section  2.3))  are  much  less 
familiar  and  are  now  tmder  active  investigation  (Garcia- 
Ybarra  and  Roener,  1989).  Moreover,  some  of  the  phe¬ 
nomena  we  describe  can  be  viewed  as  extreme  cases  of  the 
Soret  effect,  which  describes  the  drift  of  molecular  species 
in  nonisotheimal  disparate  molecular  weight  systems  (see, 
e.g.,  the  review  of  Rmner  (1980)  and  Garcia- Ybarra  and 
Roener  (1989)).  In  each  case  we  briefly  illustrate  the 
consequences  of  simultaneous  energy  transfer  for  deter¬ 
mining  fine  particle  deposition  rates  via  transport  across 
gaseous  boundary  layers  (section  3.1),  fine  particle 
“scavenging  rates”  by  evaporating  droplets  (as  in  “wet 
scrubben”;  Cooper,  1^2)  (section  3.2),  arid  the  coagulation 
dynamics  of  a  population  of  unequal  temperature  radia¬ 
tion-heated  particles  (section  3.3).  Further  details  on  each 
of  these  developing  topics  will  be  found  under  Literature 
Cited,  with  our  present  conclusions  and  recommendations 
sununarized  in  section  4.  To  some  extent  this  paper  ntay 
also  be  viewed  as  a  companion  to  (and  extension  Rosner 
(1987),  which  emphasized  the  non-Fickian  (-Brownian) 
nature  of  small  particle  mass  transport  in  the  presence  of 
convective  energy  transfer,  with  applications  to  electric 
power  generation  and  advanced  materials  processing. 

2.  Examples  of  Heat-Transfer  Effects  on  the 
Dynamics  of  Small  Particles  Suspended  in  Gases 

Consider  the  local  environment  of  a  solid  particle  too 
small  to  appreciably  “sediment”  or  depart  (due  to  its  in¬ 
ertia)  from  the  local  gas  mixture  mass  motion.  In  the 
presence  of  an  energy  flux  it  is  not  difficult  to  anticipate 
that  its  normally  “isotropic”  Brownian  motion  becomes 
“biased”  in  a  direction  related  to  the  direction  of  the  local 
energy  flux.  When  one  averages  over  many  gas  mole¬ 
cule/particle  surface  impacts,  this  “bias”  is  associated  with 
a  net  particle  “drift"  (called  phoresis)  which,  for  isotropic 
spheres,  is  usually  in  the  same  direction  as  the  local  energy 
flux  vector.  The  most  familar  case  of  Fourier  heat  con¬ 
duction  through  the  gas  is  first  briefly  discussed  in  section 
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2.1,  whereas  the  case  of  radiation  energy  flux  through  a 
locally  isothermal  gas  containing  absorbing  particles  is 
taken  up  in  section  2.2.  In  section  2.3  we  consider  some 
ununial  properties  of  asymmetric  particles  (e.g.,  aggregates 
of  two  spheres)  not  in  dermal  equilibrium  with  the  local 
carrier  gas.  In  all  cases  the  resulting  particle  motion  in 
the  dire^on  of  the  net  energy  flux  is  dnonatically  altered, 
giving  rise  to  the  engineering  consequences  discussed  in 
section  3. 

2.1.  Thermophoresis  of  Spherical  Particles  Sus¬ 
pended  in  a  Gas.  A  small  spherical  particle  which  finds 
itself  in  a  radiation-free  ’^uid”  sustaining  a  temperature 
gradient  (and  hence  Fourier  energy  diffusion)  comprises 
the  simplest  and  best-studied  case.  Two  limiting  subcases 
are  easiest  to  visualize,  depending  on  the  ratio  of  the  host 
gas  mean  free  path  1,  to  the  c^acteristic  particle  di¬ 
mension  (say  the  diameter  dp  for  a  sphere).  When  this 
ratio  (the  s&^ed  particle  Knudaen  number)  is  very  laj^e, 
the  particle  behaves  like  a  “giant*  (almost  stationary) 
molecule  in  a  disparate  molecular  weight  gas  mixture 
(Garcia- Ybarra  and  Rosner,  1989).  Waldmann  (1961) 
showed  that  such  a  particle  will  drift  down  the  temperature 
gradient  with  the  quasi-steady  velocity  (relative  to  the  gas) 

Vt  -  (aTD)pl4"F/(*,T)|  (1) 

where  q'V  >  k,|-(grad  7^|  is  the  Fourier  energy  flux  vector 
and  the  '^ermophoretic  diffusivity*  (a^DL  was  found  to 
be 


(ajD),  -  (3/4)i»,|l  +  (T/8)a^|-‘  (2) 


V.  »  (M/p)g  is  the  local  gas  momentum  diffusivity 
(^kinematic  viscosity’*)  and  the  ga8/«oIid  momentum 
accommodation  coefficient.  Equation  2  has  recently  been 
shown  to  be  consistent  with  oburved  “particle-free”  zone 
thicknesses  measured  using  counterflow  laminar  diffusion 
flames  containing  suspended  submicron  TiOg  particles 
(Gomez  and  Rosner  (1992)  and  section  3.4).  In  the  op¬ 
posite  limit  Kn-  «  1  emphasized  in  what  follows  (sections 
2.2  and  2.3),  the  surrounding  gas  can  be  treated  as  a 
“continuous”  fluid  but  one  must  abandon  the  conventional 
“no-slip”  boundary  condition  of  classical  fluid  mechanics 
(see,  eg.,  Rosner,  1986, 1989b).  Indeed,  J.  C.  Maxwell  in 
1879  showed  that  whenever  0.^01  ^  0  the  presence  of  a 
tangential  gas  temperature  gradient  in  the  vicinity  of  such 
a  surface  wfll  be  associated  with  a  local  gas  motion  directed 
from  cold  to  hot,  which  will  act  to  “propel”  the  surface 
toward  the  cooler  gas.  This  “pumping”  action,  which  can 
be  viewed  as  the  “mechanism”  of  thermophoresis  (and 
photophoresis,  section  2.2  below)  in  the  Knp  «  1  limit, 
leads  to  the  result  (see,  e.g..  Brock,  1962) 


la-jD). 


I, 

1  +  -=(A  +  B«-«*/<.) 
a 


-:|  l  +  2^  +  2c,-i| 
a  \  *p  ‘a/ 


(3) 


for  an  isotropic  sphere  thermal  conductivity  kp  in  a  gas 
of  (translational)  thermal  conductivity  k,.  (Remaining 
quantities  are  defined  in  the  Nmnenclature  section). 

It  is  uiteresting  that  (orDIp  in  both  Kn^  limits  is  actually 
particle  size  independent,  with  d.  influencing  “only”  the 
nature  of  the  transition  between  tnese  two  limiting  values 
(see,  e.g..  Figure  1,  constructed  for  silica  particles  in  ca. 
1500  K  “air”  using  a  semiempirical  transition  regime  re¬ 
lation  (Talbot,  1981)).  Equation  3  reveals  that,  unless  Kru, 
»  1,  particles  with  the  smallest  thermal  conductivity  wiO 
experience  the  largest  thermophoretic  effects. 


PARTICLE  DIAMETER  ,  dp,  ( pm) 

Figure  1.  Predicted  particle  thermophoretic  diffusivity  (normalized 
by  the  gas  thermal  diffusivity)  as  a  function  of  particle  size  and 
pressure  level  in  ca.  1S(X)  K  “air*  (after  Roaner  (1980)  and  Fernandez 
de  la  Mora  and  Rosner  (1981)). 


2.2.  Pkotophoresis  of  Spherical  Particles  Sus¬ 
pended  in  a  Gas.  Photophoresis  of  spherical  particles  is 
in  many  respects  similar  to  thermophoresis  b^use  the 
particle  motion  is  induced  by  a  temperature  gradient  in 
the  gas  adjacent  to  the  particle.  In  pkotophoresis,  however, 
the  gradient  arises  from  nonuniform  abst^tion  incident 
radiation  within  the  particle;  Le.,  the  gas  temperature  “far* 
from  the  sphere  can  be  uniform.  Considering  that  radia¬ 
tive  transport  is  often  the  dominant  mode  of  heat  transfer 
in  large-s^e  combustors  (see,  e.g.,  Viskanta  and  Menguc, 
1987)  the  “driving  force”  for  particle  photophoresis  in 
high-temperature  environments  can  be  signiflcant 
In  a  manner  completely  analogous  to  thermophoresis  (cf. 
eq  1),  the  photopboretic  velocity  of  a  particle  can  be  for¬ 
mally  expressed: 

Vp  »  (opD)p{q"R/(fc,7’))  (4) 


where  4"r  >*  the  radiative  energy  flux  vector  and  (apD)p, 
deSned  by  eq  7,  is  the  photophoretic  diffusivity  (Mack- 
owski,  1989;  Castillo  et  al.,  1990).  Prediction  of  (apD)p 
must  take  into  account  the  particular  Knudsen  number 
regime  of  the  sphere,  as  well  as  the  nature  of  radiation 
abrorption  within  the  sphere.  For  monochromatic  radia¬ 
tion  the  latter  is  characterized  by  the  size  parameter  x  ^ 
vdp/X  (where  X  is  the  radiation  wavelength)  and  m,  the 
complex  index  of  refraction  (n  -f  tk).  Particles  having  size 
parameters  much  less  than  unity  will  not  undergo  appre¬ 
ciable  photqphoretic  motion  because  radiation  absorption 
within  the  particle  becomes  uniform  in  this  limit.  Con¬ 
sequently,  in  typical  thermal  radiation  environments 
photophoresis  will  only  be  significant  for  super-micron¬ 
sized  particles,  especially  in  the  “window”  of  particle  sizes 
before  the  onset  of  appreciable  inertial  impaction  effects 
(near-critiral  Stokes  number).  The  analysis  of  photo¬ 
phoresis  under  atmospheric  pressure  (or  high  pressure) 
conditions  can  thus  be  approached  from  a  near-continuum 
(i.e.,  Kn.  «  1)  viewpoint.  In  this  regime  (Mackowski, 
1989) 


(apD)p  =  - 
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1  -I-  -:(A  -I-  Be-^*/*.) 
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Figure  2.  Ratio  of  pfiotophoretic  to  thermophoretic  'diffusion' 
factors  for  equal  Fourier  and  radiative  energy  fluxes  (after  Castillo 
et  al.  (1990)). 

which  contains  the  familiar  Cunningham-rMillikan  cor¬ 
rection  factor  to  the  Stokes  drag  (Buckley  and  Loyolka, 
1989).  Ji  is  the  dimensionless  radiative  asymmetry  factor 
(characterizing  the  degree  of  nonuniformity  of  radiation 
absorption  within  the  particle),  defined  by  (Yalamov  et 
al.,  1976) 

J,  -  (3/(2o®q"H))f  cos  9  sin  9  de  dr  (6) 

Jo  " 

O"'(r,0)  is  the  radiative  absorption  rate  distribution.  In 
^e  limit  x  »  1  most  absorption  occurs  essentially  at  the 
sphere  surface  and  Ji  approaches  its  miriimnm  value,  -0.5. 
Our  evaluation  of  J]  (analogous  to  computing  radiative 
cross  sections)  for  arbitrary  x  and  m  utili^  Lorentz-Mie 
theory  to  calculate  0'"(r,d),  for  which  eq  6  can  be  ana¬ 
lytically  (term  by  term)  integrated.  The  photophoretic 
diffusivity  for  spectrally  distributed  radiation  is  obtained 
through  integration  of  Ji  over  the  wavelength  distribution 
of  the  incident  radiation. 

Presented  in  Figure  2  is  the  ratio  (ap/a-f)  of  the  pho- 
tophoretir  and  thermophoretic  factors  vs  sphere  diameter 
for  both  highly  absorbing  (e.g.,  carbonaceous)  and  weakly 
absorbing  (e.g.,  fly  ash)  particles.  In  either  case  this  ratio 
is  equivalent  to  the  ratio  of  photo-  and  thermophoretic 
velocities  under  equal  radiative  and  Fourier  (conductive) 
heat  flux  conditions.  For  simplicity,  the  illustrative  results 
displayed  pertain  to  a  wavelength-independent  refractive 
index  and  a  blackbody  incident  radiation  spectrum  cor¬ 
responding  to  1800  K.  Ashlike  particles  display  the  in¬ 
teresting  phenomenon  of  "negative”  photophoresis  for 
diameters  less  than  about  25  ton;  i.e.,  the  "drift"  is  directed 
toward  the  radiation  source  (a  consequence  of  radiation 
"focusing”  by  the  particle  to  the  rear  (nonilluminated 
side)).  However,  phoretic  ratios  for  highly  absorbing 
particles  are  positive  throughout  and  several  times  larger 
than  the  weak-absorber  ratios.  Although  absolute  values 
of  this  ratio  do  not  exceed  about  0.1,  the  radiative  flux  to 
heat-exchanger  surfaces  in  large  furnaces  will  often  exceed 
the  "conductive”  flux  by  over  lO-fold.  Under  such  con¬ 
ditions  photophoresis  can  significantly  augment  the  de¬ 
position  rate  of  micron-sized  absorbing  particles  (Castillo 
et  al.,  1990),  as  shown  in  Figure  3. 

2,3.  "Thermal  Phoresis”  of  Asymmetric  Particles. 
In  a  general  sense  photophoresis  is  the  result  of  thermal 
nonequilibrium  between  the  particle  and  surrounding  gas 
(via  radiative  absorption  or  emission)  and  an  asymmetry 
in  the  particle-gas  heat-transfer  process.  In  "conventional” 
photophoresis  (outlined  above)  the  asymmetry  arises  from 
an  anisotropic  distribution  of  radiation  absorption  and 
results  in  a  "space-fixed”  phoretic  motion.  But,  it  is  also 
entirely  possible  for  heat-transfer  as>’mmetry  to  result  from 
an  inherent  asymmetry  of  the  particle,  such  as  non- 
sphericai  shape,  and/or  nonuniform  surface  properties. 


PARTICLE  DIAMETER  ,  dp.  ( ^m) 

Figure  3.  Relative  importance  of  particle  mass  transport  by  the 
mechanisms  of  Brownian  diffusion  (reference  case),  thermophoreais, 
and  photophoresis  for  equal  radiative  and  Fourier  energy  flux  acroaa 
laminar  gaseous  boundary  layers.  Size  dependence  for  absorbing 
particles.  T.  «  1000  K.  -  0.1. 

Phoresis  arising  from  particle  asymmetry  will  be  aligned 
with  respect  to  the  particle,  i.e.,  "body-fix^”.  Such  effects 
have,  indeed,  been  recently  observed  in  a  series  of  inter¬ 
esting  experiments  by  Rohatschek  (1985),  who  concluded 
that  photophoretic  motion  of  irregularly  shaped  particles 
can  be  3  or  4  orders  of  magnitude  greater  than  that  pre¬ 
dicted  for  uniform  property  spherical  particles. 

Few  theoretical  investigations  of  body-fixed  photopho¬ 
resis  (BFP)  have  been  published,  due,  in  part,  to  mathe¬ 
matical  difficulties  associated  with  the  analysis  of  non- 
spherical,  nonuniform  particles.  Rohatschek  (1985) 
presented  a  free-molecular-limit  (Kn^  »  1)  analysis  in 
which  the  particle  was  modeled  as  a  perfect  sphere  but 
with  a  nonuniform  surface  thermal  accommodation  coef¬ 
ficient.  Forces  arising  from  small  differences  in  thermal 
accommodation  were  found  to  be  well  in  excess  of  those 
due  to  nonuniform  particle  radiation  absorption.  We  re¬ 
cently  extended  the  same  model  to  the  near-continuum 
regime  with  similar  conclusions;  i.e.,  resulting  BFP  forces 
greatly  exceed  space-fixed,  conventional  photophoretic 
forces. 

One  interesting  and  singular  aspect  of  this  sphere-based 
model  is  that  body-  and  space-flxed  phoretic  actions  will 
not  produce  a  torque  on  ^e  particle.  Rotation  of  such  a 
sphere  would  therefore  be  due  to  Brownian  action  alone, 
and  the  resulting  BFP  motion  would  be  stochastic  in 
character.  To  a  fixed  observer,  the  long-time  phoretic 
motion  would  thus  resemble  a  "diffusion”  process.  An 
order-of-magnitude  analysis  of  the  "effective”  translational 
diffusion  resulting  from  BFP  leads  to  the  expression 

Dpjett  -  (▼bFp)*/ ^pjot  (7) 

where  Vbfp  is  the  photophoretic  velocity  and  the 
Brownian  rotational  diffusion  coefficient.  Present^  in 
Figure  4  are  near-continuum  model  predictions  for 
(normalized  with  respect  to  the  gas  momentum  diffusivity 
O  vs  sphere  diameter  at  several  ambient  gas  temperatures. 
These  results  pertain  to  a  difrerence  in  thermal  accom¬ 
modation  coefficient  of  only  0.1  between  the  two  halves 
of  the  sphere,  and  particle/gas  thermal  nonequilibrium 
resulting  from  radiative  emission  from  the  sphere  to  a  cool 
environment.  Also  shown  is  the  "normal”  Brownian 
(translational)  diffusion  coefficient  for  a  sphere  of  the  same 
size  under  the  same  ambient  conditions.  Our  results  in¬ 
dicate  that  the  “effective”  diffusion  resulting  from  pho¬ 
tophoresis  can  be  several  orders  of  magnitude  greater  than 
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Figure  4.  'Effective'  particle  diffusion  coefficient  (normalized  by 
gas  momentum  diffusivity)  compared  to  true  Brownian  diffusivity; 
consequences  of  a  0.1  difference  in  thermal  accommodation  coeffi¬ 
cient  across  a  spherical  particle  radiating  to  an  environment  at  1(XX) 
K. 

the  radiation-independent  Brownian  counterpart.  Such 
behavior  can  have  important  consequences  (section  3)  for 
coagulation  and  deposition  rates  of  radiation  absorbing/ 
emitting  “real”  particles  in  high-temperature  environments 
and  could  abo  falsify  particle  siie  inferences  based  on  the 
experimental  method  called  “dynamic  light  scattering”. 

Analysis  of  the  phoretic  behavior  of  asymmetric-shoped 
particles  is  considerably  more  complicated  than  sphere- 
based  models.  As  a  first  step,  we  recently  investigated  the 
behavior  of  an  idealized  two-sphere  aggregate  in  the 
near-continuum  regime  (Mackowski,  1990).  In  this  model 
each  of  the  two  spheres  is  taken  to  be  completely  homo¬ 
geneous  with  respect  to  surface  properties  and  the  ab¬ 
sorption  or  emission  of  radiation.  Heat-transfer  asym¬ 
metry  is  introduced  via  differences  in  properties  (i.e.,  size, 
thermal  conductivity,  and/or  refractive  index)  between  the 
two  spheres.  A  semianal^^ical  spherical  harmonics-based 
technique  is  used  to  solve  the  heat  conduction-  and 
thermal  slip-induced  creeping  flow  equations.  Results  of 
this  analysis  parallel  those  obtained  from  the  above- 
mentioned  spherical  particle/nonuniform  accommodation 
coefficient  investigation,  in  that  radiative  absorption  or 
emission  by  the  asymmetrical  aggregate  results  in  a  BFP 
motion  directed  along  the  aggregate  axis.  In  the  absence 
of  “external  aligning”  forces  ^e  BFP  motion,  when  com¬ 
bined  with  inevitable  Brownian  rotation,  will  result  in  an 
“effective”  Brownian  (translation)  diffusion  which  can  be 
significantly  greater  than  ordinary  Brownian  diffusion. 

This  model  also  predicts  that  conventional  thermo¬ 
phoresis,  due  to  a  temperature  gradient  in  the  bulk  gas, 
will  act  to  align  the  axis  of  the  asymmetrical  aggregate  with 
the  gradient.  In  so  doing  the  BFP  motions  become  es¬ 
sentially  space-fixed  and  contribute  to  (or  diminish)  the 
thermophoretic  velocities — leading  to  a  significant  change 
in  the  predicted  “effective”  aggregate  thermophoretic 
diffusivity,  Similar  conclusions  were  reached  by 

Rohatschek  (1985)  in  that  space-fixed  phoretic  forces,  while 
considerably  smaller  than  body-fix?d  forces,  contribute 
importantly  by  aligning  the  particles  in  a  particular  di¬ 
rection.  The  phoretic  velocity  due  to  aggregate  asymmetry 
may  augment  or  sometimes  oppose  the  ordinary  thermo¬ 
phoretic  velocity;  e.g.,  radiative  absorption  from  a  hot 
environment  can  lead  to  an  overall  decrease  (even  sign 
reversal)  in  («xD)p,,ff.  Incidentally,  our  parallel  Monte 


Carlo  studies  of  particulate  deposit  microstructure  (Tas- 
sopoulos  et  al.,  1989)  reveal  that  the  alignment  of  incoming 
nonspherical  particles  (relative  to  <|"f)  will  lead  to  more 
anisotropic  deposits. 

Currently,  we  are  extending  the  above-mentioned  ag¬ 
gregate  model  to  three  or  more  contacting  spheres.  These 
extensions  should  prove  useful  because,  among  other 
things,  they  will  allow  investigation  of  the  transport  be¬ 
havior  or  particles  having  no  plane  or  axis  of  symmetry 
as  well  as  bridging  the  wide  gap  between  the  properties 
of  small  and  large  “agglomerates”  (Mountain  et  al.,  1986; 
Rosner  et  al.,  1991). 

2.4.  Effect  of  Gas/Particle  Thermal  Nonequilib- 
riom  on  Coagulation  Rates.  The  previous  remarks  have 
interesting  consequences  for  predicting  coagulation  rates 
between  particles  not  in  local  thermal  equilibrium  with  the 
gas.  Several  heat-transfer  mechanisms  act  to  drive  the 
temperature  of  a  particle  above  or  below  the  surrounding 
“transparent”  gas,  e.g.,  radiative  abeorption/emission  and 
heterogeneous  chemical  reaction.  Regardless  of  the 
mechanism,  energy  transfer  from/to  the  particle  will 
perturb  the  adjacent  temperature  field  and,  as  a  result  of 
thermophoresis,  influence  encounter  rates  with  neighboring 
particles.  Our  analysis  (Mackowski  et  al.,  1990a,b;  1991) 
begins  by  investigating  the  thermophoretic  attraction  (or 
repulsion)  between  a  pair  of  particles  separated  by  the 
distance  R.  Let  Q,  (t  =  1, 2)  denote  the  net  energy  transfer 
rate  to  each  particle  t  via  radiative  absorption,  emission, 
etc.  Furthermore,  we  assume  that  the  particles  and  gas 
are  in  a  (quasi-)  steady  state  and  neglect  higher  order 
thermal  interactions  between  the  particles.  From  energy 
conservation,  particle  2  will  experience  a  gas  temperature 
gradient  resulting  from  energy  transfer  to  particle  1  equal 
to 

|dT,/(k|2  =  -0i/{4irR%)  (8) 

where  the  direction  of  the  length  x  is  taken  to  be  from 
particle  1  to  particle  2.  An  equivalent  expression  is,  of 
course,  obtained  for  the  gradient  about  partide  1  as  a  result 
of  partide  2,  with  these  results  applying  in  both  asymptotic 
extremes  of  the  particle  Knudsen  number.  Neglecting 
“hydrodynamic”  interactions  and  inertial  effects,  we  con¬ 
clude  that  due  to  thermophoresis  the  particles  will  move 
toward  (or  away  from)  each  other  at  the  rate 

Vt,i2  =  t(“T^)l$2  +  (aTD)20ll/(4TR*ft|To)  (9) 

The  motion  is  as  if  aa  attractive  (or  repulsive)  force  acted 
between  the  two  particles,  of  magnitude 

-  ''^T.12  (®1 

where  B,  and  B2  are  the  respective  partide  “mobilities” 
in  the  surrounding  gas.  Introducing  this  effective  ther¬ 
mophoretic  force  (proportional  to  R"^),  the  analysis  be¬ 
comes  equivalent  to  that  for  coagulation  between  elec¬ 
trically  charged  particles  (Zebel,  1966).  Consequently,  the 
ratio,  Z,  of  coagulation  rates  including  thermophoresis  to 
that  due  to  Brownian  motion  alone  can  be  expressed  in 
the  familiar  form  (cf.  eq  14) 

Z  =  y/[exp(y>  -  1)  (11) 

where,  now 

y  =  RT,J2iRl2^Rl2^^B^  (12) 

and  R,2  is  the  separation  distance  at  contact  (i.e.,  Oj  -(■  02). 
When  the  particles  are  cooler  than  the  surrounding  gas, 
the  energy  transfer  rates  0,  and  the  force  Fj  ,2  are  negative 
and  the  coagulation  “correction”  factor  Z  exceeds  unity. 
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Figure  5.  Correction  to  the  Brownian  (isothermal)  coagulation  rate 
constant  due  to  energy  exchange  to/from  spherical  particles  of 
unequal  radius  (after  Mackowski  et  al.  (1990,  1991)). 


as  expected.  Particles  hotter  than  the  gas  will  have  a 
Z-factor  less  than  unity. 

Suppose  the  particles  are  out  of  thermal  equilibrium 
with  the  local  gas  due  to  radiative  absorption  or  emission, 
and  radiate  from  (about)  the  local  gas  temperature  to  a 
background  environment  of  temperature  T„v-  Then  the 
above-mentioned  heat  source  for  either  particle  will  be 

=  4»a.2aBl«.T.„v^  -  c.Vl  (13) 

where  c,  and  a,  are  the  particle  emittance  and  absorptance, 
respectively,  estimated  from  Lorentz-Mie  theory  at  the 
corresponding  radiation  emission  and  absorption  tem¬ 
peratures. 

A  plot  illustrating  the  anticipated  effect  of  radiative 
heating  or  cooling  on  particle  coagulation  rates  is  given  in 
Figure  5,  in  which  the  coagulation  factor  Z  is  given  as  a 
function  of  the  two  radii,  Oi  and  02.  The  temperatures  Tq 
and  T,nv  are  1800  and  1000  K  (and  reversed),  and  the 
particle  refractive  index  was  taken  to  be  m  =  1.6  +  0.3i, 
representative  of  carbonaceous  soot.  For  radiation  cooling, 
this  effect  leads  to  preferential  growth  of  large  particles 
via  “scavenging”  of  small  particles.  The  growth  of  radia- 
tively  heated  particles,  on  the  other  hand,  is  essentially 
“frozen"  for  supermicron  particle  radii  (see  section  3.3  and 
Mackowski  et  al.  (1990a,b,  1991)  for  the  startling  conse¬ 
quences  for  PSD  evolution.  Clearly,  similar  PSD  conse¬ 
quences  would  follow  from  the  use  of  the  Zebel-corrected 
Brownian  rate  constant  for  unipolar  charged  aerosols,  a 
topic  outside  of  the  scope  of  the  present  paper). 

3.  Some  Engineering  Consequences 

The  “phoretic"  phenomena  and  “effective"  Brownian 
diffusivities  described  in  section  2  have  important  conse¬ 
quences  in  high-temperature  “dusty  gas"  applications, 
some  of  which  are  briefly  considered  below. 

3.1.  Fine  Particle  Deposition  Rates  on  Heat-Ex¬ 
changer  Surfaces.  At  sufficiently  subcritical  Stokes 
numbers  (see,  e.g.,  Fernandez  de  la  Mora  and  Roener,  1981, 
Konstandopouloe  and  Rosner,  1991)  particle  Brownian 
diffusion  across  a  thin  convective  diffusion  boimdary  layer 
would  normally  govern  local  particle  deposition  rates  in 
nearly  isothermal  systems  (see,  e.g.,  Friedlander,  1977). 
But  in  nonisothermal  forced  convection  systems  with 
nearly  constant  thermophysical  properties,  the  order  of 
magnitude  of  the  ratio  of  the  thermophoretic  deposition 
rate  to  the  ordinary  convective-Brownian  diffusion  rate  is 
set  by  the  dimensionless  ratio 

-At  =  (orlpICr.  -  T,Yr^WutffieJ^r)]/lNuJRe,Sc)]  (14) 


where  (ax)p  is  the  above-mentioned  (dimensionless) 
thermophoretic  factor,  and  the  indicated  iVu-functions  are 
the  corresponding  Nusaelt  numbers  (dimensionless  coef¬ 
ficients)  for  heat  and  mass  transfer.  More  generally, 
Rosner  (1980)  has  shown  that  -By  plays  the  role  of  a  di¬ 
mensionless  “suction"  Peclet  num^r,  i.e.,  u^  Jb^/Dp,  with 
total  deposition  rates  in  the  presence  of  heat  tranter  being 
augmented  by  the  familiar  approximate  factor 

By(e*P<AT>  -  1)  (15) 

For  the  most  commonly  encountered  systems  (with  Dp  « 
O  this  convenient  result  is  known  to  overestimate  particle 
oeixieition  rates  to  cooled  solid  surfaces  by  a  factor  (caUed 
F(outer  sink)  by  Rosner  (1980)  and  Gokoglu  and  Roener 
(1984))  equal  to  the  steady-state  reduction  in  particle  nuum 
fraction  across  the  thermal  boundary  layer  (smaller  than 
unity  but  larger  than  T^fTf  Correlations  of  F(outer  sink) 
can  be  obtained  from  the  results  of  Goren  (1977),  (Okoglu 
and  Roener  (1984, 1986a,b),  Batchelor  and  Shen  U985)  and 
Shen  (1985)  for  laminar  boundary  layers,  and  Roener  and 
Fernandez  de  la  Mora  (1984)  for  turbulent  boundary 
layers.  Noting  that,  frequently,  Nuh/Nu„  ~  (fV/Sc)*/®, 
eq  14  can  be  rewritten  in  the  form 

-Bt  =  \(arD)^/r,]\(T,  -  T,)/TJ(Pr)‘/»(Sc)*/3  (16) 

But,  since  (a^D^fv^  and  (Pr)*/^  are  frequently  of  order 
unity  (cf.  eq  2,  when  Kn^  »  1),  one  sees  that  the  magni¬ 
tude  of  -Bt  is  often  sethy  the  product  (AT/rj(Sc)^/®, 
which  can  be  quite  large  even  for  submicron  p^cles  (e.g., 
for  dp  =  1  Min  at  p  =  1  atm  in  1500  K  “air"  Sc  =  3  x  10®, 
so  that  (Sc)*/®  =  2  X  10^).  Indeed,  in  the  experiments  of 
Rosner  and  Kim  (1984)  (deposition  of  ca.  0.7-Mm  MgO(s) 
particles  from  seeded  combustion  gases  to  Pt  ribbons  at 
T^/T,  =  0.7  and  modest  Reynolds  numbers)  actual 
(thermophoretically  enhanced)  particle  deposition  rates 
were  found  to  be  over  10®  times  those  expected  in  the 
corresponding  isothermal  flow  (see,  also,  Figvire  3).  It  is 
also  significant  that  absolute  soot  particle  deposition  rates 
of  Roener  and  Kim  (1984),  Eisner  and  Rosner  (1985, 1986), 
and,  recently,  Makel  and  Kennedy  (1991),  as  well  as 
“dust-free  zone  thicknesses”  observed  near  laminar  diffu¬ 
sion  flame  “sheets"  (Gomez  and  Rosner,  1992;  see,  also, 
section  3.4  below),  are  consistent  with  kinetic  theory 
predictions  of  (aTD)p  for  single  spheres  (see  eq  2  and 
Waldmann  (1961))  even  though  it  is  well  established  that 
the  particles  in  sudi  flames  are  comprised  of  large  numbers 
of  snuill  quasi-spherical  “primary"  piuticles  bound  together 
(see,  e.g.,  Megaridis  and  Dobbins,  1990).  The  explanation 
for  this  interesting  apparent  paradox  appears  to  be  that, 
while  the  Brownian  diffusivity  and  drag  properties  of  ag¬ 
gregates  are  quite  sensitive  to  aggregate  size  (e.g.,  radius 
of  gyration)  and  morphology,  the  orientation-averaged 
thermophoretic  diffusivity,  ((axDlp),  is  remarkably  in¬ 
sensitive  to  aggregate  size  and  morphology  in  either  the 
free-molecule  or  continuum  regime  (Rosner  et  al.,  1991), 
a  phenomenon  which  facilitates  engineering  calculations 
of  total  deposition  rates  and  makes  “thermophoretic 
particle  sampling”  (see,  e.g.,  Dobbins  and  Megaridis,  1987) 
particularly  attractive. 

It  should  also  be  remarked  that  further  rate  enhance¬ 
ments  are  expected  in  “cold-wall"  systems  with  nonne- 
gligible  mainstream  particle  mass  loading  (as  in  the  de¬ 
position  of  optical  waveguide  “pre-forms"  (Rosner  and 
Park,  1988;  Park  and  Rosner,  1989a-c)).  Moreover,  one 
expects  that  the  microstructure  (and  hence  microstruc¬ 
ture-sensitive  properties  like  the  effective  thermal  con¬ 
ductivity  tensor)  of  such  deposits  will  be  different  from 
those  grown  by  ordinary  convective  diffusion  (see,  e.g.. 
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Tassopouloe  et  al.,  1989).  This  will  influence  the  thermal 
resistance  growth  rate,  of  decisive  interest  in  ‘fouling" 
situations  (Marner,  1989). 

3J2.  Hot  Dusty  Gas  "Cleanup":  Scavenging  of  Fine 
Particles  by  an  Evaporating  Droplet.  Dust-laden  gases 
are  said  to  "scrubbed"  clean  when  they  are  "contacted" 
by  a  spray  whose  droplets  capture  the  suspended  dust  and 
concentrate  the  particulate  solids  in  an  effluent  slurry. 
Most  discussions  and  measurements  of  the  performance 
of  wet  scrubbers  (see,  e.g.,  (hooper,  1982;  Hidy,  1984)  focus 
on  the  inertial  capture  of  supermicron  dust  particles  which 
"impact"  upon  the  (much  larger)  moving  droplets — a 
mechanism  only  weakly  influenced  by  the  energy  fluxes 
necessarily  associated  with  droplet  evaporation.  In  view 
of  the  considerations  of  section  2.1,  it  is  of  interest  to  ask 
to  what  extent  the  heat  transfer  associated  with  droplet 
evaporation  accelerates  the  capture  of  submicron  dust,  i.e., 
"respirable"  dust  particles  too  small  to  be  captured  by  the 
mechanism  of  inertial  impaction.  The  situation  has  some 
similarities  to  that  discussed  in  section  3.1  but  with  one 
important  distinction,  viz.,  whereas  section  3.1  dealt  im¬ 
plicitly  with  an  internally  cooled  solid  surface,  we  are  now 
(in  the  case  of  each  droplet)  dealing  with  an  evaporatively 
cooled  surface  undergoing  outward  vapor  mass  transfer. 
Thus,  whereas  the  temperature  gradient  associated  with 
the  gas  supplying  the  latent  heat  of  evaporation  tends  to 
thermophoretically  draw  small  suspended  particles  into 
the  droplet,  the  radial  outflow  of  vapor  (Stefan  flow)  tends 
to  "convect"  the  particles  away.  The  following  simple 
analysis  establishes  the  conditions  under  which  thermo- 
phoretic  particle  capture  occurs  and  predicts  the  rate  of 
this  capture  per  droplet.  While  a  more  general  analysis 
is  possible,  we  limit  ourselves  here  to  the  simplest  case  of 
negligible  particle  Brownian  diffusion,  small  particle  mass 
loadings  (wp  «  1)  and  negligible  "slip"  between  the  local 
gas  and  the  droplet  (Re'^^  «  1).  For  further  details  and 
discussion,  see  Rosner  and  Mackowski  (1992). 

If  one  considers  a  quasi-steady  evaporating  sphere,  then 
the  mass  flow  rates  m"p(4irr^  and  m"  (4irr^),  and  the  radial 
energy  flow  rate  <}'V(4irH)  are  independent  of  radius,  r, 
where  4"f  (*k|(-grad  D)  is  the  local  Fourier  heat  flux  and 
the  subscripts  p  and  v  refer  to  suspended  particles  and 
vapor,  respectively.  The  Fourier  flux  at  the  vapor/liquid 
interface  supplies  the  required  latent  heat  of  evaporation, 
L  (actually,  eq  17  below  is  valid  under  QS-transient  con¬ 
ditions,  in  whi^  case  the  effective  value  L^i  (=(-4 'V/ hi ")») 
is  the  true  latent  heat  L,  plus  a  contribution  associated 
with  droplet  heatup),  linking  -<)"f  and  m",,  but  -q'V  is 
also  linked  to  the  p^icle  thermophoretic  flux  m  via  eq 
1.  Combining  these  interrelations,  and  eliminating  the 
energy  flux,  leads  to  a  remarkably  simple  and  interesting 
result  (Rosner  and  Mackowski,  1992): 


-fhp  *  fh,«p 


^  1 
Le,  CpT  "  ^1 


(17) 


linking  the  flne  particle  capture  rate  -rhj,  with  the  pre¬ 
vailing  single  droplet  evaporation  rate  rhy-  ui  the  bracketed 
term  appears  the  thermophoretic  diffusivity  ratio  a-rLe^ 
displayed  in  Figure  1,  and  the  more  familiar  Lewis  number 
for  vapor  diffusion  D,/  {k/pCp)^,  which  is  quite  near  unity 
(e.g.,  for  water  vapor  in  air).  It  is  interesting  that  for  the 
bracketed  term  to  be  positive  the  ambient  temperature 
cannot  exceed  a  readily  calculated  threshold  (ca.  900  K  for 
a  water  droplet  collecting  submicron  particles  from  ordi¬ 
nary  combustion  products).  For  process  gas  conditions 
above  this  threshold,  the  Stefan  flow  associated  with 
droplet  evaporation  prevents  the  thermophoretic  sca¬ 
venging  of  dust  particles,  but  much  below  this  threshold 


ambient  temperature  droplets  will  be  "surprisingly"  ef¬ 
fective  in  removing  small  suspended  particles  (in  accord 
with  eq  17),  far  exceeding  the  performance  expected  based 
on  flne  particle  capture  by  Brownian  diffusion  alone.  Also, 
below  this  threshold  gas  temperature,  particles  larger  than 
some  "cutoff  diameter"  (obtained  via  Figure  1,  by  setting 
the  bracketed  term  in  eq  17  equal  to  zero)  will  not  be 
scavenged  in  this  manner.  Incidentally,  by  integrating  both 
sides  of  eq  17  over  the  total  evaporative  lifetime  of  a 
droplet,  the  bracketed  factor,  which  above  plays  the  role 
of  a  discriminant,  can  also  be  shown  to  be  interpretable 
as  the  ratio  of  the  actual  total  number  of  particles  cap¬ 
tured  during  the  evaporative  lifetime  of  an  individual 
motionless  droplet  to  the  number  of  particles  initially 
contained  in  the  volume  required  to  accommodate  the 
initial  droplet  mass  if  the  droplet  density  were  reduced 
to  the  ambient  gas  density.  Further  interesting  corollaries 
of  eq  17  will  be  found  in  Rosner  and  Mackowski  (1992). 

3.3.  Coagulation  Dynamics  of  a  Population  of 
Unequal  Temperature  Spherical  Particles.  To  dem¬ 
onstrate  the  effect  of  gas/particle  thermal  nonequilibrium 
on  the  evolution  of  a  particle  size  distribution  (PSD),  we 
consider  an  aerosol  comprised  of  spherical  particles  (of 
volume  o)  evolving  due  to  coagulation  only.  In  this  case 
the  PSD  n(i;,t)  =  diVp/dc  obeys  the  integro-differential 
equation  (see,  e.g.,  Friedlander,  1977) 

dn(v,t)/dt  =  (1/2)  [  Z<ti,u>  8(v,u)  n{v-u)  n{u}  du  - 
Jo 

n(v)  f  Z(v,u}  b(v,u)  n(u}  du  (18) 
Jo 

where  we  have  introduced  Z,  the  above-mentioned  cor¬ 
rection  factor  for  thermophoretic  scavenging/repulsion  (eq 
11)  and  0{v„vj)  is  the  coagulation  rate  constant  due  to 
ordinary  Brownian  coagulation.  To  account  for  coagulation 
in  the  transition  regime,  we  simply  use  the  "harmonic 
average  approximation"  for  the  collision  kernel  (see,  e.g. 
Pratsinis,  1988;  Rosner  and  Tassopoulos,  1989),  viz.,  0  = 
0TM0c/i0FM  ^c)>  where  uid  dc  the  coagulation 
rate  constants  ("kernels")  in  the  free-molecule  and  con¬ 
tinuum  regimes,  respectively.  The  three  lowest  integral 
moments  of  eq  18  satisfy  the  well-known  DDEs  (see,  e.g., 
Cohen  and  Vaughn,  1971) 

dAf(/d/  =  -(l/2)f  J  Ziu.u}  P(i.',u)  n(i')  n{u)  df  du  (19' 
Jo  “ 

dAf,/d/  =  0  (20' 

dMj/dt  =  I  j  Ziv.u}  P(o,u)  n(f)  n(u)  vu  dt  du  (21) 
Jo  ® 

Here  is  the  ^th  moment  of  the  distribution  deflned  by 

Af»  -  J”o*n(u)  do  =  Afo'^/exp((*®/2)  In^  o,>  (22) 

where  the  second  equality  in  eq  22  is  a  useful  property  of 
"log-normal"  distributions.  Note  that  the  zeroth  and  first 
moments  of  the  PSD  correspond,  physically,  to  the  total 
particle  number  density,  Np,  and  the  total  particle  volume 
fraction,  <t>^  respectively.  In  situations  for  which  the  PSD 
shape  remains  "log-normal”,  eqs  19-21  become  a  closed 
set  of  coupled  DDEs  (see,  also,  Rosner  and  Tassopoulos, 
1991),  used  below  to  predict  the  effects  of  thermophoresis 
on  aerosol  population  evolution. 

In  Figure  6  we  plot  the  PSD  standard  deviation,  <r^,  and 
mean  particle  diameter  (corresponding  to  the  median  of 
the  log-normal  PSD)  vs  dimensionless  time  t  ~  N^  qI  (see 
the  Nomenclature  section)  where  the  evolution  occurs  from 
an  initial  PSD  characterized  by  a^iO)  =  3.0  and  d,(0)  =  0.01 
um.  The  curves  labeled  1  correspond  to  the  isothermal 
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Figure  6.  Energy-trsnsfer  effects  on  the  evolution  of  a  log-normal 
population  of  unequal  temperature  (‘overheated')  apherical  parti¬ 
cles:  time  dependence  of  population  ‘spread*  and  mean  particle 
diameter  (after  Mackorraki  et  aL  (1990a,b,  1991)). 

case  (particle  encounters  due  only  to  Brownian  motion) 
while  curves  labeled  2  correspond  to  the  case  of 
‘overheated’  particles.  Note  that  in  the  isothermal  case 
the  population  spread  (solid  curve  1)  eventually  reaches 
a  constant  value  corresponding  to  a  ‘self-preserving’  dis¬ 
tribution  (see,  e.g.,  Lee,  1983).  In  the  case  with  heat 
transfer  ‘the  larger  the  particles,  the  larger  the  repulsive 
forces’,  so  small  particles  in  the  population  tend  to  ‘catch 
up’  with  the  larger  ones,  resulting  in  a  narrower  PSD  with 
a  smaller  mean  size  (dashed  curve  2).  In  view  of  the 
above-mentioned  fact  that  electrostatic  repulsion  and 
‘thermal  repulsion"  lead  to  interparticle  collision  rate  laws 
of  identical  form  (cf.  eq  11),  this  PSD  ‘sharpening’  (or 
‘focusing’  in  size  space)  should  also  be  observed  in  the 
coagulation  of  unipolar  aerosols.  In  either  case  the  ability 
to  obtain  narrower  distributions  than  those  ordinarily 
associated  with  Brownian  coagulation  may  prove  to  be 
useful  in  the  controlled  synthesis  of  fine  particles  (section 
4).  Cases  involving  ‘underoooled’  particles  (Tp  <  T,)  may 
also  be  of  considerable  practical  interest  (Mackowski  et 
al.,  1990a,b;  1991),  but  they  require  a  more  general 
treatment  due  to  significant  departures  from  the  single¬ 
mode  log-normal  PSD-shape  assumption.  We  have  re¬ 
cently  employed  a  discretized  population  balance 
(‘sectional’)  method  based  on  the  work  of  Hounslow  et 
al.  (1988)  to  examine  the  behavior  of  radiation-cooled 
particle  populations  (Mackowski  et  aL,  1990a,b,  1991).  In 
Figure  7  we  show  the  evolution  of  the  population  median 
particle  volume  vs  dimensionless  time  for  an  initially 
monodisperse  aerosol  coagulating  under  isothermal  (solid 
line)  and  undercooled  (dashed  line)  conditions.  After  a 
short  ‘induction"  period,  after  which  particle  sizes  become 
disparate  enough  so  that  thermophoretic  forces  dominate 
(see,  also.  Figure  5),  the  larger  particles  grow  at  signifl- 
cantly  higher  rates,  thus  appreciably  increasing  the  pop¬ 
ulation  average  size. 

3.4.  Thermophoretic  Gas  Cleaners.  When  the 
‘target”  surface  is  hotter  than  the  dusty  gas  stream  (i.e., 
when  T.  >  T,),  the  incoming  particles  are  thermally 
‘repelled*  and  the  associated  thermophoretic  ‘suction* 
parameter,  -Bj,  changes  sign.  This  implies  dramatically 
reduced  particle  deposition  rates  on  ‘overheated’  surfaces, 
and  in  the  limit  of  negligible  inertia  and  Brownian  diffu¬ 
sion,  this  can  bring  a^ut  a  nearly  ‘dust-free*  zone  in  the 
immediate  vicinity  of  the  surface  (see,  e.g.,  Friedlander, 
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Figure  7.  Energy-transfer  effects  on  the  evolution  of  an  initially 
monodisperse  population  of ‘undercooled*  spherical  partidea;  time 
dependence  of  mean  particle  diameter  (after  Mackowski  et  aL 
(1990a,b)). 

1977;  Gokoglu  and  Rosner,  1986a,b;  Friedlander  et  aL, 
1988;  Strattmann  et  al.,  1988;  Park  and  Rosner,  1987). 
This  can  be  exploited  in  many  ways  as,  for  example,  in  a 
continuous  “thermophoretic  gas  cleaner*.  Tlius,  if  the  gas 
in  the  ‘dust-free*  zone  is  continuously  ‘skimmed  oft*,  it 
can  be  conveyed  to  dust-sensitive  equipment  (e.g.,  a  com¬ 
bustion  turbine)  whereas  the  remainmg  gas  will,  of  course, 
be  dust-enriched.  Some  of  the  properties  of  an  axisym- 
metric  laminar  stagnation  flow  ‘hot-wall*  configuration, 
including  the  effects  of  high  particle  mass  loading  and 
intermediate  Reynolds  number,  are  computed  in  Park  and 
Rosner  (1987, 1989a-c).  In  effect,  of  course,  this  is  the 
same  phenomenon  which  acts  to  ‘shut  ofP  the  further 
coagulation  of  ‘large*  overheated  particles,  as  described 
above  (section  3.3). 

4.  Conclusions  and  Recommendations 

From  the  representative  examples  discussed  in  section 
3  and  the  references  cited  therein,  we  conclude  that  in  the 
presence  of  radiative  and/or  diffusive  (Fourier)  energy 
transfer  small  suspended  particles  in  a  gas  experience 
unbalanced  forces  which  cause  them  to  drift,  usually  (but 
not  always)  in  the  direction  of  the  local  energy  flux.  This 
“plunesiB*  can  cause  dramatic  changes  in  the  rates  at  adiich 
fine  suspended  particles  are  deposited  (section  3.1)  or 
scavenged  (section  3.2),  as  well  as  the  rates  at  wUch 
unequal  temperature  partides  coagulate  among  themselvee 
(section  3.3).  Particularly  interesting  and  not  yet  e^lwed 
in  depth  are  the  phenomena  predicted  here  for  unsym- 
metrical  particles  in  the  presence  of  a  heat  flux  (section 
2.3),  e.g.,  unequal  sphere  agglomerates,  which  should  ex¬ 
hibit  ‘anomdous*  thermophoretic  coefficients  and 
‘effective*  Brownian  diffusivities,  in  addition  to  a  tendency 
to  ‘align*  themselves  in  the  direction  of  the  energy  flux. 
Such  phenomena,  several  of  which  are  being  experimen¬ 
tally  studied  in  our  (HTCRE)  laboratory,  will  inevitably 
have  important  implications  for  the  morphology  of  ag¬ 
glomerates  suspended  in  gases,  the  deposition  (capture) 
rates  of  single  particles  and  agglomerates,  and  the  mi¬ 
crostructure  (hence  thermal  properties)  of  resulting  de¬ 
posits. 

Because  of  the  particle  mass  transport  mechanisms 
described  here,  well-documented  methods  which  have  been 
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successfully  used  to  predict  or  correlate  coagulation  and 
deposition  rates  in  “dusty  gas'  systems  without  appreciable 
heat  transfer  cannot  be  straightforwardly  applied  to  situ¬ 
ations  with  appreciable  radiative  and/or  diffusive  energy 
fluxes.  Conversely,  by  imposing  suitable  energy  fluxes, 
many  of  these  phenomena  can  be  creatively  exploited  to 
improve  the  performance  of  devices  intended  to  produce 
or  separate  fine  particulate  matter  inf  from  flowing  gases. 
We  hope  that  the  present  contribution  will  stimulate  ad¬ 
vances  along  these  lines. 
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Nomenclature 

a  «  spherical  particle  radius 
-Bj  ■  thermophoretic  “suction"  parameter,  eq  14 
B,  *  mobility  of  particle  i  in  the  prevailing  gas,  eq  10 
Cc  *  Cunningham-Millikan  correction  to  Stokes  drag 
<^aom  *  momentum  (tangential)  exchange  coefficient 
c^  *  thermal  exchange  coefficient  (gas/solid) 

Cj,  a  specific  heat  of  carrier  gas 
a  >  diameter 

Fi2  ■  effective  interparticle  force 

J,  »  photophoretic  radiative  asymmetry  factor 

Ji  *  wavelength  integrated  photophoretic  radiative  asymmetry 

factor 

k  »  Fourier  thermal  conductivity 

kg  *  Boltzmann  constant 

k  «  imaginary  part  of  complex  number  m 

Kn  »  Knudsen  number  (based  on  /,  and  particle  diameter) 

L  *  latent  heat  of  vaporization 

Le  «  Lewis  number  (Fick/thermal  diffusivity  ratio) 

1.  *  gas  mean  free  path 

*  kth  moment  of  aerosol  PSD  n(i;,t);  eq  22 
m  *  complex  index  of  refraction  of  particle 
n  *  real  part  of  m 

n  =  PSD  function  (of  u);  n  =  dNj,/dv 
AL  =  total  particle  number  density  (=Afo) 

Nuf,  =  heat-transfer  Nusselt  number 
Nu„  =  mass-transfer  Nusselt  number 
O'"  =  radiative  energy  absorption  rate/volume 
Q,  =  net  heat  source  for  particle  i 
=  Fourier  energy  flux  vector 


“  radiative  energy  flux  vector 
p  ^  pressure 

Pr  *  Prandtl  number  of  gas  mixture  (momen:um-/ther- 
mal-diffusivity  ratio) 

R  =  center-to-center  distance  between  two  spheres 
r  K  radius  measured  from  center  of  spherical  particle  or 
droplet 

Re  *  Reynolds  number 

Sc  “  Schmidt  number,  e/D.fBrownian) 

T,  =  temperature  at  outer  edge  of  thermal  ‘boundary  layer’ 
T.  *  surface  temperature  at  gas/solid  interface 
T.  >  temperature  far  from  surface  (“upstream  infinity') 
Vj  =  particle  thermophoretic  velocity 
Vp  K  particle  photophoretic  velocity 

V  >  particle  volume  (treated  as  a  continuous  variable) 

X  B  radiative  size  parameter,  rdp/X 

X  =  distance  between  particle  1  and  particle  2  (section  2.4) 
y  =  dimensionless  thermophoretic  “potential" 

Z  =  correction  to  the  Brownian  coagulation  rate  constant 

Greek  Symbols 
a  =  absorptivity 

a-r  =  thermophoretic  diffusion  factor 
^mom  ~  momentum  accommodation  factor 
d  =  coagulation  rate  constant,  eq  18 
6  =  boundary  layer  thickness 
<  «  emittance 

>  aerosol  particle  volume  fraction 
M  =  dynamic  viscosity  of  host  gas 
v,  K  momentum  diffusivity  (kinematic  viscosity) 

Wp  =  particle  mass  fraction  in  mixture 
p  >  mass  density 

T  =  dimensionless  time  (proportional  to  hL^t) 

(Tg  =  standard  deviation  of  log-normal  PSD 
(Tg  ^  Stefan-Boltzmann  radiation  constant 

Subscripts 

BFP  =  body-fixed  phoresis 
eff  =  effective 

env  =  environment  (radiative) 
g  =  gas 

h  =  heat  transfer 
i  =  particle  i 

k  =  moment  index  (k  =  0,  1,  2, ...) 
m  ~  mass  transfer 

p  =  associated  with  particle  (dimension) 
rot  =  rotational 

V  =  vapor 

w  =  evaluated  at  the  gas/target  interface 
«>  =  evaluated  far  from  the  target 

Abbreviations 

BFP  =  body-fixed  phoresis 

C  =  continuum 

FM  =  free  molecule 

ODE  =  ordinary  differential  equation 

OWG  =  optical  waveguide 

P  =  photophoretic 

PSD  =  particle  size  distribuvion 

QS  =  quasi-steady 

T  =  thermophoretic 
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Abstract— Aerosol  sampling  from  industrial  environments  (e.g.  combustion  engines)  or  natural 
environments  (e.g.  the  troposphere)  frequently  involves  conveying  the  sample  to  a  downstream 
('sheltered')  instrument  via  an  upstream  tube  or  duct.  While  the  instrument  may  be  capable  of 
characterizing,  say.  the  particle  size  distribution  (PSD)  of  the  aerosol  actually  presented  to  it,  the 
investigator  is,  of  course,  usually  more  interested  in  the  PSD  of  the  aerosol  entering  the  upstream 
sampling  tube.  Invariably,  this  differs  from  that  measured  because  of  several  systematic  phe¬ 
nomena  —perhaps  the  two  most  obvious  of  which  are  particle  size-dependent  losses  to  the  tube  walls 
(i.e.  incomplete  'penetration')  and  PSD  distortion  due  to  suspended  ptarticle-particle  coagulation 
when  the  particle  concentrations  are  sufficiently  high.  We  show  here  how  recent  research  on  the  use 
of 'moment  methods'  to  predict  the  effects  of  size-dependent  walls  loss  and/or  Brownian  coagulation 
in  flow  systems  can  now  be  brought  to  bear  to  conveniently  solve  this  ‘inverse*  problem  by 
numerically  integrating  the  quasi-one-dimensional  coupled  moment  equations  in  the  upstream 
direction,  using  downstream  (measured)  aerosol  properties  in  the  definitions  of  all  dimensionless 
dependent  variables  and  parameters.  Illustrative  'universal*  graphs  are  presented  here  for  the 
sampling  of  log-normally  distributed  'inertialess*  (Brownian)  aerosols  in  long  straight  adiabatic 
ducts  for  both  commonly  encountered  extremes  of  particle  Knudsen  number  Kn^  >  I(free  molecule) 
or  Kn^  I  (continuum),  as  well  as  convenient  rational  approximations  derived  from  the  leading 
terms  of  a  Taylor  series  expansion  of  the  above-mentioned  dimensionless  moment  equations. 


NOMENCLATURE 


A  cross-sectional  area  of  duct 

h  exponent  describing  dependence  of  St„  on  particle  volume  v 
C  dimensionless  coagulation  rate  parameter  (fm  or  c) 

C,  friction  factor  (dimensionless  wall  shear  stress) 

C'l*'  first  order  coefficient  in  Taylor  series  expansion  of  M.fi) 
d,  particle  diameter  corresponding  to  r, 

D  Brownian  diffusion  coefficient  for  particle  of  volume  r 
d.  duct  diameter 

/^(u.  r)  function  of  particle  volumes  appropriate  to  moment  k 
F^  dimensionless  moment  shift  function  (equation  (12)) 

F,  entrance  effect  correction  for  streamwise  averaged  mass  transfer  coefficient 
k  moment  index  (e.g.  It  =  0. 1, 2  are  of  particular  physical  interest) 
kg  Boltzmann  constant;  Table  2 

K  environment/mechanism-dependent  part  of  the  coagulation  rate  constant  p 
Knp  Knudsen  number  based  on  prevailing  gas  mean-free  path  and  particle  diameter 
L  total  length  of  sampling  duct  (upstream  of  detector) 
mass  of  particle  of  volume  r 
W,  It"*  moment  of  n(v, .  .  . )  s  dA/p/dr 
n  size  distribution  function  dfi/p/dr 
Np  total  particle  number  density  (  e  Af,) 

Nu„  mass  transfer  Nusselt  (Sherwood)  number  (Rosner.  1986) 

P  ‘wetted*  perimeter  of  the  duct  of  cross-sectional  area  A(:) 

•  Supported,  in  part,  by  AFOSR  (Grand  89-0223),  DOE-PETC  (Grant  DE-FG-22-9  0PC90099)  and  the  1991 
Yale  HTCRE  Lab  Industrial  Affiliates  (Union  Carbide.  Shell,  SCM  Chemicals,  DuPont,  General  Electric), 
t  Professor  of  Chemical  Engineering,  Yale  Univ.;  Director  HTCRE  Lab. 
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X 

X 

2“' 

Greek  letters 

P 

i 

Pp 

Pt 

Ppp. 

V 

o. 


Reynolds  number  based  on  4/4//*  and  mean  velocity  U 

radius  (of  curvature)  of  centerline  of  a  bend  (e.g.  elbow)  in  duct 

particle  Schmidt  number  litlp\^lD 

mass  transfer  Stanton  numlxr,  Nu./(ReSc)  (Rosner,  1986) 

time 

particle  stopping  time  in  the  prevailing  carrier  gas 

u./tp/v:  dimensionless  particle  stopping  time 

absolute  temperature  (K) 

mean  velocity  of  carrier  gas  in  sampling  duct 

particle  volume  (dummy  variable) 

friction  velocity  (C,/2)'  ^1/ 

volume  of  spherical  particle 

geometric-mean  particle  volume  in  log-normal  PSD  n(v, . . .) 
mean  particle  volume.  of  local  aerosol 

L—z\  distance  measured  upstream  from  the  detector 
position  vector 

distance  measured  downstream  from  the  sampling  tube  inlet 

collision  integral  functions  defined  by  coagulation  rate  law  (equation  (IS)) 


coagulation  rate  constant  in  “mass-action"  law  (equation  (14)) 
scaled  dimensionless  distance  upstream  of  the  aerosol  detector 
particle  volume  function  (  s  M, ) 
intrinsic  density  of  a  particle  (nip/rp) 
dimensionless  k"‘  moment  of  PSD  (Po  =  K  Pi  ~  I) 
dynamic  Newtonian  viscosity  of  the  carrier  gas 
momentum  diffusivity  of  the  carrier  gas, 

(ReSc)~ '  (/./d.)  (in  equation  (40)) 

geometric  standard  deviation  of  the  log-normal  PSD  ii(v, . . . )  (see,  e.g.  Rosner  and  Tassopoulos, 
1989) 


Subscripts  and  superscripts 

C  pertaining  to  interparticle  coagulation 
c  coagulation 
eff  effective 

k  pertaining  to  moment  k  where  k  need  not  be  an  integer  nor  positive 
(ft)  pertaining  to  moment  k 
g  pertaining  to  log-normal  PSD  (geometric) 
gas  pertaining  to  the  carrier  gas 
L  evaluated  at  ;  (downstream  location  of  aerosol  detector) 
max  maximum  value 
meas  measured 
p  particle(s) 

o  evaluated  at  the  sampling  system  inlet  station  (r^O) 

WL  pertaining  to  wall  loss 

'  first  derivative  with  respect  to  { 

*  second  derivative  with  respect  to  ( 

Others 

(*)  normalized  by  the  value  (  )  u  the  aerosoi^etector  (at  z  L) 

(')  mean  value  (e.g.  f  w  0p/Afp,  d  =  {6S/n)''^;  St„  (streamwise  mean)) 

Abbreviations 

c  continuum  (KUp  <  I) 
coag  pertaining  to  (caused  by)  coagulation 
fm  free  molecule  (KOp  p  I) 

GDE  general  dynamic  equation  (population  balance) 
h.o.t.  higher  order  terms 
ODE  ordinary  differential  equation 
PSD  particle  size  (volume)  distribution 
RHS  right  hand  side  (of  equation) 

WL  pertaining  to  (caused  by)  particle  loss  to  walls 


1.  INTRODUCTION,  BACKGROUND  AND  MOTIVATION 

Historically,  the  need  to  accurately  sample  industrial  and  natural  aerosols  has  led  to  many 
advances  in  the  science  and  technology  of  fine  particulate  matter  suspended  in  gases, 
including  the  invention  and  development  of  a  wide  variety  of  clever  instruments  as  well  as 
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ancillary  theoretical  developments  concerning  techniques  for  avoiding  or  minimizing  ‘bias' 
when  extracting  the  sample  and  delivering  it  to  the  instrument  (see.e.g.  Fuchs,  1964;  Davies, 
1966;  Mercer,  1973;  Friedlander,  1977;  Hinds.  1982;  Hidy.  1984;  Vincent,  1989).  Some 
phenomena  will  cause  an  unavoidable  systematic  ‘distortion'  in  the  recorded  particle  size 
distribution  (PSD)  which  must  be  corrected  for — in  particular — the  size-dependent  loss  of 
particles  to  the  walls  of  an  upstream  sampling  tube,  and  any  suspended  particle-particle 
coagulation  that  occurs  between  the  inlet  of  the  sampling  tube  and  the  (downstream) 
monitoring  instrument.  The  present  contribution  provides  rational,  yet  simple,  methods  to 
make  such  corrections,  under  conditions  that,  while  deliberately  idealized,  are  frequently 
encountered.  Indeed,  one  of  our  goals  is  to  produce  a  set  of  generalized  graphs  which  can 
facilitate  making  such  corrections  for  Brownian  aerosols  (in  the  absence  of  appreciable 
inertial-  and  sedimentation-phenomena),  as  well  as  guide  the  design/selection  of  sampling 
systems  which  will  keep  future  corrections  acceptably  small. 

Wall  deposition  losses  in  the  absence  of  interparticle  coagulation  have  been  the  subject  of 
many  early  investigations,  especially  for  the  case  of  dilute,  steady,  fully-developed  laminar 
gas  flow  in  straight,  constant-area  ducts  (see,  e.g.  Fuchs  and  Sutugin,  1971).  This  under¬ 
standing.  in  fact,  led  to  the  development  of  the  so-called  ‘diffusion  battery'  (see,  e.g.  Mercer 
and  Green.  1974)  for  estimating  PSDs  based  on  the  observed  penetration  characteristics  of 
the  aerosol  as  a  function  of  particle  size,  duct  diameter,  length,  and  gas  flow  rate.  At 
suspended  particle  concentrations  high  enough  for  non-negligible  coagulation  (Brownian, 
shear.  .  .  . )  the  situation  becomes  rather  more  complicated  and  the  computational  methods 
used  by  earlier  investigators  to  predict  PSD  distortion  effects  fail  (i.e.  assuming  that  the 
penetration  of  one  particle  size  class  is  not  influenced  by  the  simultaneous  presence  of  the 
others).  Indeed,  the  interesting  interactions  between  coagulation  and  convective-diffusion 
deposition  rates  now  comprise  an  active  field  of  research  (see,  e.g.  Park  and  Rosner,  1989; 
Pratsinis  and  Kim,  1989;  Bai  and  Biswas.  1990;  Rosner  and  Tassopoulos,  1990),  to  which 
the  present  paper  contributes.  The  present  work  is  also  an  extension  and  application  of 
techniques  we  have  recently  developed  for  predicting  total  mass  deposition  rates  from 
polydispersed  aerosol  populations  formed  by  coagulation  (Rosner,  1989;  Rosner  and 
Tassopoulos.  1989).  but  without  the  complication  of  appreciable  coagulation  within  diffu¬ 
sion  boundary  layers  (see.  e.g.  Park  and  Rosner.  1989;  Bai  and  Biswas.  1990). 

Finally,  mention  should  be  made  of  relevant  theoretical  studies  whose  primary  purpose 
was  to  support  interpret  fundamentally-orienteu  experim<*nts  on  the  mechanisms  of  or¬ 
ganic-  and  inorganic-soot  production  in  laminar  flames  (see,  e.g.  Ulrich  and  Subramanian, 
1977;  Ulrich  and  Riehl,  1982;  Dobbins  and  Mulholland.  1984;  Dobbins  and  Megaridis, 
1987;  Frenklach  and  Harris,  1987;  Megaridis  and  Dobbins,  1990;  Megaridis.  1987;  Zachar- 
iah  and  Semerjian.  1989).  Indeed,  our  treatment  below  of  free-molecule  regime  coagulation 
in  aerosol  sampling  system  is  equivalent  to  that  of  Dobbins  and  Mulholland  (1984)  (rather 
than  that  of  Lee  el  al.,  1984).  suitably  augmented  to  include  the  simultaneous  occurrence  of 
particle  size-specific  convective-diffusion  to  the  walls  of  the  sampling  tube. 

Rather  than  explicitly  dealing  with  the  'direct'  problem  of  predicting  the  ‘evolution'  of  the 
penetrating  (and  depositing)  aerosol  as  a  function  of  dimensionless  downstream  distance 
and  coagulation  parameters  (see.  e.g.  Rosner  and  Tassopoulos,  1990,  1991b)  we  show  here 
(section  3.6)  that  by  'backwards'  (upstream-)  integration  of  the  governing  moment  equations 
and  the  use  of  dimensionless  dependent  variables  and  parameters  based  on  downstream 
reference  quantities,  we.  in  effect,  can  provide  direct  solutions  to  the  canonical  'inverse' 
problem  of  aerosol  sampling  theory — i.e.  correction  factors  which  allow  the  investigator  to 
go  immediately  from  PSD-parameters  recorded  at  the  downstream  instrument  to  the  PSD- 
parameters  that  must  have  characterized  the  aerosol  entering  the  inlet*  of  the  sampling 


*  Additional  correction  factors  may  be  necessary  due  to  peculiarities  of  the  inlet  orifice  itself,  and  its  geometric 
surroundings  (see,  e  g.  Vincent.  1989,  Mercer,  1973;  Hungal  and  Willeke,  1990  ;  Belyaeu  and  Levin,  1974;  Zabel, 
1978;  Ivie  ei  al .  1990;  lluebcrt  ft  ill.  1990).  In  considering  the  need  for  such  additional  factors,  however,  the  reader 
should  bear  in  mind  that  we  emphasi/e  here  Brownian  particles  well  below  the  'inertial'  range,  tn  any  case,  if 
necessary,  additional  'inlet'  correction  factors  can  be  superimposed  on  the  present  predictions. 
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tube.  Toward  this  end,  in  section  2  we  spell  out  our  underlying  assumptions/idealizations, 
and  in  section  3  we  state  the  laws  used  here  to  predict  the  evolution  of  the  three  lowest 
integral  moments  of  the  log-normal  shape  aerosol  PSD  for  (quasi-)  steady,  (quasi-)  one¬ 
dimensional  gas  flow  in  a  straight  adiabatic  aerosol  sampling  duct.  General  results 
(‘universal'  graphs  and  formulae)  are  presented  in  section  4,  together  with  a  simple 
numerical  (‘backward-marching')  procedure  for  dealing  with  much  more  complex  cases. 
Some  implications  and  applications  of  the  results  of  section  4  are  provided  and  briefly 
discussed  in  section  5  with  the  help  of  two  specific  numerical  air-monitoring  examples 
(nearly  isothermal) — one  ground-based  and  one  airborne.  Section  6  offers  a  concise  defense 
of  many  useful  idealizing  assumptions,  as  well  as  an  indication  of  how  to  relax  many  of  the 
assumptions  exploited  here  for  particular  sampling  applications  which  are  more  complex, 
such  as  extractive  sampling  from  high  temperature  combustors  (for  further  details  the 
reader  is  referred  to  the  references  cited  in  section  6  and  Rosner  and  Tassopoulos  (1990, 
1991b)).  Our  summary  recommendations,  and  an  indication  of  future  relevant  work  are 
contained  in  section  7,  which  concludes  the  present  paper.  Extensions  based,  in  part,  on 
section  6  will  be  the  subject  of  future  communications  f^rom  this  laboratory. 

2.  BASIC  ASSUMPTIONS 

In  keeping  with  our  goal  of  deriving  reasonably  general  results  applicable  to  a  wide 
variety  of  intrinsically  similar,  commonly  encountered  sampling  situations  we  introduce  the 
following  explicit  simplifying  assumptions  and  idealizations.  These  are  defended  or  shown 
to  be  relaxable  in  section  6. 

Al:  Quasi-one  dimensional  steady  gas  flow  in  a  straight  adiabatic  sampling  duct  with 
constant  effective  diameter  d,„  =  4A/P. 

A2;  Aerosol  PSD  well-represented  by  a  single  mode  continuous  log-normal  function  of 
the  dense  spherical  particle  volume  u. 

A3:  Axial  diffusion  of  particles  negligible  compared  to  axial  convection. 

A4;  Previously  calculated  (or  measured)  fully-developed  perimeter-averaged  mass  trans¬ 
fer  coefficients  adequately  describe  local  convective-diffusion  particle  mass  transport  to  the 
wall,  and  depend  on  a  single  power,  b,  of  the  particle  volume  valid  in  the  high  Schmidt 
number  (Sc  =  v/D^  P  1)  limit. 

AS:  Negligible  interfacial  ‘barrier’  to  incident  particle  capture  by  the  wall,  and  negligible 
reentrainment  of  already  deposited  particles. 

A6:  Brownian  coagulation  occurs  primarily  in  the  core  of  the  sampling  duct  flow 
according  to  the  mass-action  laws  of  free-molecule  transport  (KOp^  1)  or  continuum 
transport  (Knp<^  1). 

As  discussed  in  section  6,  Assumptions  A3,  A4,  and  A6  can  apply  to  either  laminar  or 
turbulent  duct  flow  (perhaps  more  accurately  to  the  latter  than  the  former  (because  of  A4, 
see  sections  5.2  and  6)).  Methods  to  include  the  systematic  consequences  of  ‘entry'  effects, 
non-power-law  deposition,  transition  regime  particle  transport  and  coagulation,  diabatic 
walls,  bends  (in  a  ‘piecewise-straight'  sampling  system)  and  non-spherical  'aggregated'  solid 
particles  are  taken  up  briefly  in  section  6,  and  in  greater  detail  in  Rosner  and  Tassopoulos 
(1990, 1991).  It  should  also  be  remarked  that  in  the  continuum  limit  (KUp  1)  many  of  our 
methods/results  will  apply  to  particle  suspensions  in  flowing  liquid  (hydrosals). 

3.  MOMENT  METHODS  AND  LOG-NORMAL  PSDs 

Rather  than  dealing  directly  with  the  complete  suspended  particle  population  balance  or 
‘general  dynamic  equation’  (GDE)  governing  the  continuous  size  distribution  function 
n(v,  X,  f)  =  dNp/dv  (presumed  to  be  measurable  at  the  location  of  the  downstream  aerosol 
instrument)  the  predictions  made  below  are  based  on  solving  a  closed  set  of  coupled 


Sampling  error  correction  for  tube  flow 


847 


diflerential  equations  for  the  three  lowest  integral  ‘moments'  of  n{v, . . .),  defined  by: 

Mt  =  J  . .  .)dp  (k  =  0,1.2).  (1) 

This  strategy,  and  its  implementation  for  the  case  of  steady  quasi-one  dimensional  duct 
flow,  are  the  subjects  of  sections  3. 1-3.6  below.  Owing  tc  more  extensive  discussions 
elsewhere  (see,  e.g.  the  mass  deposition-oriented  work  of  Rosner,  1989;  Rosner  and 
Tassopoulos,  1989,  1990,  1991)  only  the  most  important  conclusions/results  are  stated 
below  for  completeness. 

3. 1 .  Method  of  moments 

The  use  of  diflerential  equations  for  the  moments  of  n(p, . . .)  to  construct  rational 
approximate  solutions  to  the  GDE  (an  integro-dilTerential  equation)  is  a  subject  that  has 
been  under  intense  development  in  the  past  two  decades,  and  progress  continues  unabated. 
Moment  methods  have  been  especially  successful  in  elucidating  many  aerosol  situations 
involving  coagulation  (see,  e.g.  Cohen  and  Vaughn,  1971;  Lee  et  al.,  1984;  Dobbins  and 
Mulholland,  1984;  Megaridis  and  Dobbins.  1989, 1990;  Pratsinis  and  Kim,  1989;  Frenklach 
and  Harris,  1987, . .  .).  Although  different  authors  use  somewhat  different  methods  to 
achieve  closure,  especially  when  coagulation  occurs  near  or  in  the  ‘free-molecule'  regime, 
almost  all  authors  derive/use  equations  for  Mq,  M,,  and  Mj,  which  are,  indeed,  parti¬ 
cularly  simple  (see  sections  3.4  and  3.5).  The  first  two  moments  also  have  the  merit  of  being 
of  direct  physical  significance,  with  Mo  =  Np  (the  total  particle  number  density)  and 
Mi  =  4>p  (the  total  particle  volume  fraction).  In  the  present  case  of  steady  quasi-one- 
dimensional  flow  in  a  sampling  duct  we  therefore  focus  our  attention  on  the  coupled 
ordinary  differential  equations  governing  the  M^fz),  where  k  =  0, 1, 2,  and  z  is  the  axial 
coordinate  down  the  sampling  tube  (measured  from  the  inlet).  When  use  is  made  of  the 
known  laws  of  particle  deposition  (section  3.4)  we  encounter  terms  involving  additional 
moments  M^.  where  k^0A,2  so  that  it  becomes  necessary  to  interrelate  moments  to 
obtain  a  closed  set  of  integrable  equations.  In  the  present  case  this  is  accomplished  using  the 
log-normal  PSD  shape  approximation,  as  discussed  in  section  3.2,  and  defended  in  section  6 
for  this  class  of  problems. 


3.2.  Implications  of  log-normality 


It  is  now  well  known  that  'coagulation-aged'  populations  of  aerosol  particles  closely 
approach  single-mode  log-normal  size  distributions  (i.e.  Gaussian  in  ln(t;)).  Moreover,  it  is 
easy  to  demonstrate  that  power-law  wall  removal  processes  (see  A4  and  section  3.3)  will  not 
alter  this  situation.  Hence,  for  our  present  purposes  the  single-mode  log-normal  approxima¬ 
tion  is  a  particularly  appropriate  and  convenient  starting  point.  For  such  single-mode  log¬ 
normal  PSDs  it  can  also  be  shown  that  all  moments  are  necessarily  interrelated  as  follows; 

Mt  =  i;Npexp^^y^ln*a,^.  (2) 


where  is  the  median  (geometric  mean)  particle  volume  and  o,  is  the  corresponding 
geometric  standard  deviation  of  the  PSD.  These  PSD  parameters  can  be  related  to  Np. 
and  Mi  via: 


and 


(3) 


I 


> 


(4) 
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but  any  other  desired  moment  2)  can  now  be  obtained  in  terms  of  either  N^,  4>p, 

and  M2,  OT  A/p,  c,  and  a^.  We  will  also  have  frequent  need  for  the  relation: 


r,  =  p-exp 


(-(O'"'”-)- 


(5) 


where  v  s  ^p/A/p  is  the  mean  particle  volume. 

While  not  discussed  further  here,  the  assumption  of  log-normality  in  effect  also  allows  us 
to  escape  the  numerical  problems  often  encountered  in  ‘inverting*  aerosol  population 
measurements  (see,  e.g.  Cooper  and  Wu,  1990).  Thus,  our  “upstream”  integrations  (section 
3.6)  are  perfectly  straightforward,  but  deliberately  limited  to  cases  where  the  required  PSD- 
parameter  ‘correction*  factors  are  well  within  one  decade  of  unity. 


3.3.  Modeling  particle  size-dependent  particle  wall  losses 

For  steady  quasi-one-dimensional  (2-direction)  flow  in  a  duct  of  area  Alz)  and  wetted 
perimeter  P(z)  the  wall  loss  contribution  can  be  written  (see,  e.g.  Rosner  and  Tassopoulos, 
1990) 

St„,(p)p*n(p)dp.  (6) 

A  QZ  jwalllon  A  Jo 

where  Stp,(p)  is  the  appropriate  perimeter-mean  particle  volume-dependent  dimensionless 
mass  transfer  coefficient  (Stanton  number  =  Nu„/(ReSc)),  Rosner  (1986)  at  station  z,  and 
r(p,  2)  is  the  PSD  describing  the  bulk  (mixing-cup  averaged)  suspended  particle-gas 
mixture  at  this  same  streamwise  station.  As  shown  and  exploited  in  our  recent  work  (see,  e.g. 
Rosner,  1989;  Rosner  and  Tassopoulos,  1989)  in  the  submicron  particle  size  range  of 
primary  interest  here  it  is  often  possible  to  accurately  represent  the  v-dependence  of 
St„(D, . . .)  by  simple  power-law: 

(7) 


where,  e.g.  for  Sc  |>  1  turbulent  convective-diffusion  (see,  e.g.  Rosner  (1986,  1990)  equation 
(6.5-1  lb))  Nu„  —  Sc°  St„  ~  Sc"®-’®^  and  therefore: 

and,  for  fully-developed  laminar  flow  Nu„  =  const,  St„  Sc~  so  that: 


b  = 


(9) 


Since  for  dense  spherical  particles  in  the  free-molecule  regime  (Knp>  1),  and 

D  ^  f  ~  for  dense  spherical  particles  in  the  continuum  limit,  the  exponents  b  used  in  our 
calculations  below  are  readily  obtained  and  summarized  in  Table  1.  For  isothermal 
convective-diffusion  they  are  seen  to  fall  between  the  extremes  of  —0.235  and  —0.667. 
Other  mechanisms  of  particle  transport  to  the  sampling  tube  wall  (e.g.  eddy  impaction,  or 
particle  thermophoresis)  and  other  types  of  particles  (e.g.  aggregates  comprised  of  smaller 
diameter  ‘primary*  particles)  will,  of  course,  be  characterized  by  rather  different  6-values,  as 
discussed  elsewhere  (Rosner,  1989;  Rosner  and  Tassopoulos,  1989;  Rosner,  1991). 

Taking  into  account  this  power-law  behavior,  and  considering  the  particular  case  of 
constant  area  adiabatic  duct  flow,  we  find  that  equation  (6)  can  be  simplified  to: 


- -t/  suffi  r- 


(10) 


where,  in  the  equations  below  k  =0,1,2.  In  practice  (see,  e.g.  section  3.5)  we  rewrite 
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Table  I.  Values  of  the  exponent  b  for  fully-developed  convective  diffusion  in 

straight  ducts 


Flow  regime 

Particle  Knudsen  number 

Laminar* 

Turbulent '-smooth  wall 

Free-molecule  (fm) 

-2/3 

-0.469 

Continuum  (c) 

-1/3 

-0.235 

*  In  this  case  Nu.  ^  constant  so  that  St„  ^  Sc" 

t  Based  on  Sc  >  I  transport  by  Brownian  motion  within  the  viscous  sublayer. 


equation  (10)  in  the  rormally  simple  way: 


I-*?’!. 


where  is  the  dimensionless  ‘moment-shift’  function*  defined  by: 

_  I* 

‘  ■  M*  /I* 

where  we  have  introduced  the  dimensionless  moments: 

— - — -In'ff.J. 

Equations  (12)  and  (13)  follow  from  the  PSD-condition  of  log-normality. 


(11) 


(12) 


(13) 


3.4.  Inclusion  of  coagulation  processes  at  high  particle  number'densities 

If  the  coagulation  rate  between  n(r)dt'  particles/volume  of  volume  v±dv/2  and  n(u)du 
suspended  particles/volume  of  volume  u  ±  du/l  is  taken  to  be  of  the  ‘mass-action’  form 
^(u,  r)n(u)n(r)dudt-  then  it  is  known  that  (in  a  well-mixed,  transient  situation): 

I  /t(M,r)-^(M,t)ii(p)  n(w)dt’du.  (14) 

Clearly  /o  =  —  I  and  /,  =  0,  and  it  can  be  shown  that  /j  =  2ui’  (see,  e.g.  Cohen  and  Vaughn 
(1971)).  Following  Dobbins  and  Mulholland  (1984),  for  log-normal  distributions,  once  the 
rate  constant  0(u,  r)  is  specified  the  right-hand  side  of  equation  (14)  can  always  be  calculated 
numerically  and  the  results  expressed  in  the  form: 

where  the  coefficient  K  is  independent  of  particle  size  and  number  density  but  dependent  on 
environmental  conditions  and  perhaps  particle  material  properties.  Values  of  K  and  the  so- 
called  collision  integrals  Z***  for  the  free-molecule  (fm)  coagulation  of  hard  spheres  (Dob¬ 
bins  and  Mulholland,  1984)  are  reproduced  (in  our  notation)  in  the  top  rows  of  Table  2.  In 
the  continuum  (c)  limit  it  is  well  known  that  (dMo/dt)co.,  is  proportional  to  — (Mq 
-I-  1,3)  and  (dM^/dtlc^,  is  proportional  to  (Mf  +  (see,  e.g.  Cohen  and 

Vaughn  (1971)  or  Flagan  and  Seinfeld  (1988))  with  the  proportionality  constant  (hence  K- 
value)  depending  on  the  combination  T/p,„,  When  these  results  are  recast  in  terms  of 
equation  (IS)  using  equation  (S)  we  obtain/use  the  values  of  X,  and  2^*’  shown  in  the 
bottom  row  of  Table  2.  Thus,  in  the  equations  below  we  set  the  Brownian  coagulation 
contribution  to  UdMJdz  equal  to  the  RHS  of  equation  (IS),  and,  for  a,  <  3.32,  used  the 


*  In  our  previous  papers  (Rosner,  1989;  Rosner  and  Tassopoulos,  1989)  particular  attention  was  focused  on  F, 
which  was  shown  to  be  equal  to  the  physically  interesting  ratio  of  the  actual  total  mass  deposition  rate  to 
the  total  mass  deposition  rate  if  all  particles  in  the  population  had  the  mean  volume  f. 


8S0 


D.  E.  RosNEt  and  M.  Tassopoulos 


Table  2.  Summary  of  factors*t{  appearing  in  the  coagulation  rate  term  [equation  (IS)] 


Coagulation  regime* 

K(T....) 

2<oi 

2“' 

2(2l 

Free-molecule*(Knp>>  1) 

-4,y2i-;'‘  expj^^ln'o.j 

0 

8>/2t>i’'‘-exp[^ln»<r.] 

4k,  r 

(  /I  \'l 

r  /i  M 

Continuum  (Knp<t  I) 

-2jl+exp(^-ln>o.j| 

0 

4p'|l+exp|^-ln>o,jJ 

*  Values  taken  from  Dobbins  and  Mulholland  (1984);  see,  also,  Megaridis  and  Dobbins  (1990a)  for  e,  >  3.32. 
t  As  a  consequence  of  the  assumed  conservation  of  volume  when  two  particles  collide  and  coalesce  Z"  ‘  vanishes  for 
all  cases.  (Recall  that  M ,  =  -  particle  volume  fraction.) 

t  Note  that,  according  to  equation  (t  5)  only  the  products  necessarily  have  the  same  dimensions  (in  the  fm  and 
c  cases).  However,  as  is  clear  from  the  entries  above,  the  K-factors  (and  the  corresponding  Z“*  factors)  themselves  do 
not  have  the  same  dimensions  in  the  free-molecule  and  continuum  cases. 

§  For  intermediate  KOp-values  (transition  regime),  see  Rusner  and  Tassopoulos  (1991b)  and  section  6.6. 


values  of  K  and  Z'*'  assembled  in  Table  2.  If  necessary,  this  approach  could  clearly  be 
extended  to  other  coagulation  mechanisms  (see,  section  6,  i.e.  other  rate  constants 
/}(u,  t’, .  .  . )  (cf.  equation  (14))  appropriate  to  the  mechanism  of  interest  (see,  e.g.  Mackowski 
et  ai,  1991). 

3.5.  Coupled  ODEs  for  free-molecule  and  continuum  coagulation;  steady,  qua<i~one~dimen- 
sional  duct  flow  with  simultaneous  wall  losses 

Combining  the  results  of  sections  3. 1-3.4  above,  we  see  that,  subject  to  the  assumptions  of 
section  2,  in  all  cases  the  PSD  moments  (k  =  0, 1, 2)  will  satisfy  the  coupled  ODEs: 

(16) 

where  f » is  the  above-mentioned  moment  shift  function  appropriate  to  the  mass  transfer 
law  (section  3.3)  and  the  appropriate  values  of  /C(r, . . .)  and  Z**’  appearing  in  the 
coagulation  rate  terms  can  be  read  from  Table  2.  All  of  the  conclusions  we  derive  and 
state  below  follow  from  these  coupled  non-linear  ordinary  differential  equations  and  the 
abovementioned  moment  interrelations.  Before  proceeding  it  is  worth  noting  that  the 
coagulation  term  (second  term  on  RHS)  will  not  appear  in  the  ODE  governing  M,  ( = 
i.e.  Z*°’  =  0.  Moreover,  if  instead  of  an  axial  position  variable  z  measured  downstream  from 
the  sampling  duct  inlet  r  =  0  (with  0  <  z  <  L)  we  introduce  a  position  variable  x  measured 
upstream  from  the  aerosol  instrument  at  z  =  L  then  x  =  L— r  and  dx=— dz.  Thus, 
introduction  of  such  a  variable  (to  form  the  corresponding  U(dMJdx}  equations)  will 
merely  change  the  sign  of  all  terms  on  the  RHS  of  equation  (16).  This  sets  the  stage  for  the 
‘backward'  integration  of  equation  (16)  to  find  the  upstream  dependence  of  the  aerosol  PSD 
moments  based  on  a  knowledge  (measurement)  of  their  downstream  (instrument)  values. 


3.6.  Non-dimensionalized  inverse  problem  ODEs 

In  Rosner  and  Tassopoulos  (1990,  1991b)  we  illustrated  the  use  of  equation  (16)  to 
calculate  the  downstream  evolution  of  the  aerosol  PSD  in  terms  of  inlet  values  and 
associated  dimensionless  coagulation  parameters  for  a  variety  of  flow  conditions  (laminar, 
turbulent,  free-molecule  transport,  continuum  transport).  To  solve  the  ‘inverse’  problem  of 
aerosol  sampling  theory,  we  merely  redefine  our  dimensionless  variables  and  parameters 
using  downstream  values  (presum^  measujed)  and  integrate  in  the  upstream  direction. 
Thus,  we  here  introduce  <A  =  </'p/</»p.L.  and  define  the 

‘rescaled'  axial  distance  variable: 


i  =  |^SUi\)f,.Lj(Z.-z). 


(17) 
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The  coupled  ODEs  satisfied  by  N{i, . . .),  •  •  •)  and  . . .)  can  be  obtained  from 

the  dimensionless  ODEs  of  Rosner  and  Tassopoulos  (1990,  1991)  by  making  the  appro¬ 
priate  parameter  replacements  and  sign  changes,  with  the  following  results.  Let  us  intro¬ 
duce  the  coagulation  parameter: 


Cs 


|(0  i/  St„(l\)f..L  fVp.L} 


(18) 


(a  measure  of  the  relative  importance  of  coagulation  and  wall  loss  in  altering  the  PSD 
parameters).  Then,  subject  to  the  ‘initial'  normalization  conditions:  fV(0)=  1,  i^(0)=  I, 
MjfO)  =  1  and  the  assumptions  of  section  2,  the  coupled  dimensionless  moment  equations 
become  the  non-linear  first  order  system: 


where  the  subscript  L  implies  evaluation  of  the  appropriate  function  at  z  =  Z.  (x  =  0).  These 
same  equations  formally  hold  for  either  the  free-molecule  or  continuum  case  provided  the 
appropriate  coagulation  parameters/functions  are  introduced  according  to  Table  2.  Sim¬ 
ilarly,  subject  to  the  assumptions  indicated  in  section  2  and  discussed  in  section  6,  they  hold 
for  either  turbulent  or  laminar  flow  in  a  straight  duct. 


3.7.  Similitude  implications 

Inspection  of  the  above-mentioned  moment  equations  governing  the  present  simplified 
model  reveals  that  the  idealizations  of  section  2  lead  to  correction  factors  of  the  functional 
form: 


^p../^p.L  =  fcto{^„„.  b,  C)  (22) 

0P../0P.L  =  fcti  (<J„,..  b,  C)  (23) 

M2.o/M2.l  =  fct2({„„,  b,  C)  (24) 

from  which  we  can  also  calculate  the  interesting  correction  factors: 

=  fctjl^n,.,,  <r,.L.  b,  C)  (25) 

‘’t.o/Vt.  L  =  fct*  (4,„„  <r,.  L,  b,  C),  (26) 


where  relevant  values  of  the  exponent  b  are  given  in  Table  1,  the  value  of  the  PSD  spread 
parameter  a,,  i,  is  presumed  to  be  measured,  is  the  value  of  the  'scaled'  length  coordinate 
{  (equation  (17))  evaluated  at  the  sampling  tube  inlet  (z=0)  and  C  is  the  relevant 
dimensionless  coagulation  parameter  (which  is  seen  to  scale  linearly  with  Np,  l).  In  section  4 
we  present  the  results  of  numerical  integrations  of  (equations  (19-21))  for  many  combina¬ 
tions  of  physical  interest  to  display  the  explicit  sensitivity  of  the  desired  correction  factors  to 
macroscopic  flow  regime  (laminar,  turbulent)  and  the  Knudsen  regime  of  particle  transport. 
We  also  present  an  extensive  set  of  graphical  results  for  the  commonly  encountered  limiting 
case  of  negligible  Brownian  coagulation  (i.e.  Afp,|,  is  low  enough  to  cause  C  <  1).  Note  that, 
although  all  of  our  graphs  deliberately  go  beyond  the  range  of  values  where  the  leading 
term  of  a  smali-{  Taylor  series  expansion  is  valid  (see  section  4.2),  rather  frequently  if„„- 
values  of  physical  interest  will  be  small  enough  to  allow  the  simple  explicit  formulae  of 
section  4.2,  thus  circumventing  the  need  for  the  above-mentioned  graphs,  or  specific 
numerical  integrations  (sec,  e.g.  the  numerical  examples  discussed  in  section  5). 
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4.  RESULTS 

4.1.  'Universal'  graphs  for  convective-diffusion  wall  losses  without  {and  with)  appreciable 
particle  coagulation 

In  this  section  we  display  and  briefly  comment  upon  the  results  of  our  numerical 
integrations  of  the  coupled  set  of  ODEs  explicitly  given  by  equations  (19H21)  in  the  text. 


Fig.  1.  Effects  of  particle  size  dependent  convective-diffusion  wall  loss,  Knudsen  number  regime, 
and  PSD  spread  on  the  particle  number  density  correction  factor  .Vp  „  l  for  turbulent  gas  flow  in 
straight^  sampling  tubes  of  scaled  length  s  in  the  absence  of  appreciable  Brownian  coagulation. 


Fig.  2.  Eflccis  of  particle  size  dependent  convective-diffusion  wall  loss,  Knudsen  number  regime, 
and  PSD  spread  on  the  geometric  mean  particle  volume  correction  factor  for  turbulent  gas 

flow  in  straight  sampling  tubes  of  sealed  length  s  in  the  absence  of  appreciable  Brownian 

coagulation. 
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subject  to  their  initial  (normalization)  conditions.  Results  are  grouped  as  follows:  Figs  1-6 
pertain  to  FSD-parameter  corrections  in  the  limit  of  negligible  coagulation,  with  the  first 
three  dealing  with  turbulent  carrier  gas  flow.  These  graphs  allow  N^-,  r,-  and  a.-corrections 
to  be  made  for  scaled  sampling  tube  lengths  1  with  either  free-molecule  or  continuum 


Fig.  3.  Effects  of  particle  size  dependent  convective-diffusion  wall  loss.  Knudsen  number  regime, 
and  PSD  spread  on  the  correction  factor  to  the  geometric  PSD-spread.  <r,.,/(T,  for  turbulent  gas 
flow  in  straight  sampling  tubes  of  scaled  length  (  in  the  absence  of  appreciable  Brownian 

coagulation. 


Fig.  4.  Effects  of  particle  size  dependent  wall  loss  by  convectiv;-diffusion.  gas  flow  regime  and  PSD 
spread  parameter  on  the  total  particle  number  density  correction  factor  for  gas  flow  in 

straight  sampling  tubes  of  scaled  length  in  the  absence  of  appreciable  free-molecule  regime 

inlcrparticle  coagulation. 
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increased  size  dependence  of  the  Brownian  diffusivity  D{v)  when  Knp->  cx/.  Figures  3-6  also 
pertain  !o  negligible  coagulation  but  compare  laminar  flow  results  to  those  for  turbulent 
flow  when  the  particle  transport  mechanism  is  free-molecule.  As  discussed  in  sections  5.2 
and  6.1,  for  the  laminar  flow  cases  ‘entrance  effects’  may  dictate  the  need  for  heff-values  not 
explicitly  included  in  these  graphs.  When  4m„<0.5  (say),  these  effects  can  be  dealt  with  by 
using  the  small  ^„,„Taylor  series  formulae  presented  in  section  4.2  (below)  and  Appendix  A. 
For  sm.x=  0(1)  one  expects  b,f{-*b  (Table  1)  and  Figs  4-6  may  be  useful  in  their  present 
form.  Note  that  the  PSD-correction  factors  depart  further  from  unity  for  laminar  flow  cases 
when  compared  at  the  same  value  of  the  scaled  sampling  system  length, 

Figures  7-9  display  our  computed  PSD-correction  factors  in  the  presence  of  appreciable 
interparticle  Brownian  coagulation  effects  (non-negligible  C)  when  the  gas  flow  is  turbulent 
and  the  PSD  spread  parameter  at  the  particle  detector  is  near  the  appropriate  ‘self- 
preserving’  value.  Clearly,  for  any  value  of  these  correction  factors  depart  further  from 
unity  when  the  coagulation  parameter  C  is  increased,  in  either  the  free-molecule  or 
continuum  regimes.  Particularly  interesting  is  the  rapid  rise  in  the  spread  correction  factor 
<7,  „/ct,  L  which  sets  in  at  sufficiently  large  values  (see  Fig.  9).  This  can  be  shown  to  be  a 
consequence  of  equation  (3),  given  the  sensitivity  of  (compared  to  and  Afj)  to 
coagulation  in  the  sampling  system. 

4.2.  Small-i  expansions 

As  will  be  appreciated  from  several  numerical  examples  (sections  5.1  and  5.2)  the 
maximum  value  of  ^  (the  ‘scaled’  distance  measured  upstream  from  the  aerosol  instrument 
to  the  sampling  tube  inlet)  is  frequently  a  small  number,  often  less  than  10~ '.  This  suggests 
that  even  the  above-mentioned  universal  graphs  or  the  numerical  integrations  on  which 
they  are  based  can  be  bypassed  by  using  a  small  ^  Taylor  series  expansion,  the  first  terms  of 
which  can  be  written  down  by  inspection  of  equations  (20,  21).  Thus,  if  we  write: 

=  1  -t-  C',*'(C,  ir,.L,  b)4  +  h.o.t.  (27) 


Fig.  7.  Eflects  of  Brownian  coagulation,  particle  size-dependent  convective-diffusion  wall  loss  and 
panicle  'gas  Knudsen  number  regime  on  the  correction  factor  for  total  particle  number  density 
A’p.„/.Vp  L  for  turbulent  gas  flow  in  straight  sampling  tubes  of  scaled  length  Curves  shown  penain 
to  nearly  self-preserving  PSDs  at  the  detector  (i.e.  =  2.46  and  (<T,,t)c  =  2.3). 


22:7-0 
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Figure  8.  Ejects  of  Brownian  coagulation,  particle  size-dependent  wall  convective-diffusion  loss 
and  particle/gas  Knudsen  number  regime  on  the  correction  factor  to  geometric  mean  particle 
volume  t',.,/i  ,.L'  for  turbulent  gas  flow  in  straight  sampling  tubes  of  scaled  length  Curves  shown 
pertain  to  nearly  self-preserving  PSDs  at  the  detector  (i.e.  ((T,.i.Vm  =  2.46  and  (<r,.L)e  =  2.3). 


Fig.  9.  Effects  of  Brownian  coagulation,  particle  size-dependent  wall  convective-diffusion  loss  and 
particle/gas  Knudsen  number  regime  on  the  correction  factor  for  to  the  PSD  spread.  <r,.,/<r,.L. 
Curves  shown  pertain  to  nearly  self-preserving  PSDs  at  the  detector  (i.e.  {a,  L)fm  =  2.46  and 

=  2.3). 
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we  find  that  the  first-order  coefficientst  C\“  are  explicitly: 

Cr=l+C  (jt  =  0) 

CV’=1  (*  =  1)  (28) 

and 

C\^’  =  exp(-2filn^<T,,L)— 2Cexp(  — In^ffjL)  (<c  =  2). 

Under  circumstances  for  which  assumptions  A I -6  are  reasonable  and  the  ‘higher  order 
terms'  of  equation  (27)  can  be  neglected  t  we  can  therefore  state  the  following  potentially 
useful  rational  explicit  correction  factor-relations: 

1 -K1 -l-C)^„„  (number  density)  (29) 

^P.L 


(independent  of  b  and  Og  i.) 


and,  since  f  s  (<t>plNp)  we  find: 


g..o  .  1  ,  (C\^’-Krr-2CV*)^, 

2lntT,., 


(spread)  (30) 


^  S  1  -  C  (mean  size).  (31) 


The  last  two  equations,  when  combined  with  equation  (S),  immediately  yield  the  corres¬ 
ponding  correction  to  the  geometric  mean  particle  size:  tyo/t-'i.L* 


i.e.  S  1  -[C-l-i(C\2'  -I-  C'.o'  -t-  C',*')]  (32) 

Thus,  it  would  appear  that  when  0.3  (say)  and  the  assumptions  Al-7  arc  simultan¬ 
eously  reasonable  equations  (29)-(3I)  can  provide  acceptable  rational  corrections  to 
measured  aerosol  data  to  account  for  the  systematic  effects  of  upstream  coagulation  and/or 
wall  losses.  It  is  interesting  to  note  that  while  in  the  explicit  forms  stated  they  appear  to  be 
valid  for  either^  laminar  or  turbulent  gas  flow,  and  free-molecule  or  continuum  particle 
transport,  in  practice  they  would  often  incur  unacceptably  large  errors  for  laminar  flow 
systems  due  to  so-called  ‘entrance’-eflects  on  the  mass  transfer  coefficient  (see,  e.g.  sections 
3.2  and  6.1,  and  Rosner  (1986)).  A  numerical  procedure  to  effectively  eliminate  such  errors 
without  introducing  any  new  parameters  has  been  developed  (Rosner  and  Tassopoulos, 
1991b)  but  is  beyond  the  scope  of  the  present  paper.  However,  the  present  results  can  be 
approximately  corrected  for  this  effect  by  simply  replacing  the  ‘fully-developed'  value  of 
Stmlft)  ( =  Nu„/(ReSc))  by  the  streamwise-integrated  average  value  f  (entrance)  •  Nu„  (fully 
developed)/(Re  •  Sc(i;L))  and  by  introducing  an  effective  value  of  the  St„ -exponent  b,  written 
which  varies  between  -4/9  and  —2/3  (Table  1)  according  to  the  value  of  F(entrance). 
This  correction  procedure  for  the  present  laminar  results,  using  Nu„(fully  developed) 
s  3.657  for  a  straight  circular  duct,  is  illustrated  in  section  S.2. 


4.3.  Direct  numerical  integration  for  particular  cases 

The  ‘universality'  of  the  results  derived  and  presented  above,  and  the  relatively  small 
number  of  dimensionless  parameters  on  which  they  depend  (section  3.7),  of  course  follow 
from  our  underlying  idealizations,  made  explicit  in  section  2.  However,  the  basic  moment- 
method  approach  (section  3),  together  with  our  treatment  of  wall  losses  (section  3.3)  and 

t  Closed-form  expressions  are  available  for  the  C/'  cocflkienis  (of  )  (sec  Appendix  A). 

t  These  physical  differences  do.  of  course,  affect  the  results,  but  only  implicitly  via  establishment  of  the 
appropriate  values  of  {...(equation  (3.6-1))  and  the  St.  exponent  h  (sec  Tabic  I). 
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interparticle  coagulation  (section  4,  and  the  primary  references  cited  therein)  can  clearly  be 
used  to  obtain  PSD-parameter  correction  factors  in  much  more  complicated  aerosol 
sampling  duct  situations,  such  as  for  example,  systems  which  are  cooled  near  the  inlet  (to 
prevent  thermal  failure)  and  subsequently  heated  (to  prevent  water  vapour  condensation) 
f^or  extracting  carbonaceous  soot  (smoke)  aerosols  from  gas  turbine  engine  combustors  (see, 
e.g.  Colket  et  al  (1977),  and  section  6).  In  such  cases  the  recommended  procedure  would  be 
to  integrate  the  appropriate  coupled  set  of  ODEs  (of  the  form  of  equation  (16),  but  with  the 
LHS  written  {(Xj A)A{AVM^)ldz}  and  introduction  of  the  ‘upstream’  position  variable  x 
=  L  —  z),  after  incorporating  the  required  special  features  of  the  sampling  system  of 
particular  interest  (see  sections  6.1-6.12). 

5.  IMPLICATIONS  AND  APPLICATIONS 

To  fix  ideas,  become  familiar  with  typical  orders-of-magnitude,  and  illustrate  how  to  use 
the  results  calculated/presented  in  section  4  it  is  useful  to  briefly  consider  at  least  two 
specific  numerical  examples  that  arise  in  conveying  aerosol  samples  to  a  ‘remote’  instrument 
via  a  constant  area  duct.  The  first  (section  5.1)  is  based  on  an  exercise  in  Friedlander’s  1977 
textbook  involving  turbulent  gas  flow  in  a  ground-based  aerosol  sampling  system.  The 
second  example  involving  laminar  gas  flow  is  based  on  conditions  encountered  in  a 
proposed  system  for  sampling  the  Earth’s  atmosphere  from  a  research  aircraft  flying  at 
ca.  8  km  altitude.  In  both  cases  we  examine  log-normal  distributions  of  highly  submicron 
particles  for  which  wall  losses  are  often  non-negligible  even  for  relatively  ‘short’  ducts.  We 
also  define  the  total  concentration  levels  beyond  which  Brownian  coagulation  processes  of 
such  particles  within  the  duct  will  become  non-negligible.  Explicit  applications  to  more 
complex  (e.g.  diabatic)  systems,  as  frequently  encountered  in  sampling  combustion  engines, 
are  beyond  the  scope  of  this  paper  (see  sections  4.3  and  6). 


5. 1 .  Turbulent  flow  in  duct  of  ground-based  air  monitoring  station 

Exercise  5,  Chapter  3  of  Friedlander  (1977)  deals  with  a  ground-based  air  monitoring 
station.  Outside  air  is  delivered  to  the  aerosol  instruments  via  a  2-in  diameter,  12-ft  long  duct 
at  a  mean  velocity  of  10  fts" '.  The  goal  is  to  “calculate  the  correction  factor  that  must  be 
applied  for  submicron  particles  as  a  result  of  diffusion  to  the  wails  of  the  duct,  e.g.  the 
percentage  by  which  the  measured  concentration  must  be  altered  to  provide  the  true 
(ambient)  concentration  as  a  function  of  particle  size”.  To  illustrate  our  procedures  let  us 
imagine  that  the  aerosol  instrument  ‘sees’  a  sub-micron  aerosol  PSD  which  is  log-normal 
with  a  spread  of  <7,  =  2.46  and  a  geometric  mean  diameter  corresponding  to  dp  =  0.02  //m, 
i.e.: 

dp  =  dp  ■  exp  [  -  6  «^,]  =  (0  02  /rm)  •  exp  [  -  i  In^  (2.46)]  =  0.01 75  ftm.  (33) 

The  Brownian  diffusivity  for  a  0.02  /jm  diameter  particle  in,  say,  20'C  air  at  1  atm  is 
(see,  e.g.  Table  2.1,  Friedlander  (1977))  1.34  x  10~*cm^s”'  and  (since  the  momentum 
diffusivity  of  air  is  1.50  x  10"'cm^s"'  under  these  same  conditions)  we  also  find  that 
Sc(Cp  l)  s  v/Dp(rp,L)  =  1.12  x  10^  for  the  mean  particle  diameter.  The  mass  transfer  coeffi¬ 
cient  St„  appropriate  to  this  case  can  be  estimated  once  the  Reynolds  number  l/d,/v 
characterizing  this  air  flow  is  determined.  Converting  to  metric  units  we  find  that 
d,  =  5.08  cm,  U  =  3.05  x  10^  cm  s"  ‘  and,  therefore.  Re  =  1.03  x  10*,  which  is  well  into  the 
turbulent  range.  For  a  smooth  wall  straight  duct  and  particles  in  this  size  range  (see  section 
6.7)  St„  can  therefore  be  estimated  from,  say,  equation  (6.5-1  lb)  of  Rosner  (1986),  i.e.: 

St„*0.0889^^y^^Sc-°''°*  (SOI),  (34) 

where  C,(Re)  is  the  friction  factor  (non-dimensional  time-averaged  shear  stress  at  the 
hydraulically  smooth  duct  wall)  corresponding  to  the  prevailing  Reynolds  number.  This  can 
be  estimated  as  0.079  Re"  so  that  C,  =  0.783  x  10"^,(Cf/2)*'^  =0.626  x  10"'  and,  hence. 
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S(„(Pp.|.)  %  3.97  X  10~‘  (via  equation  (34)).  We  are  now  in  a  position  to  calculate  the 
decisive  dimensionless  ‘scaled'  length  of  the  sampling  tube,  i.e.: 

{m„=4^^jst„(i5p.L)fo.L  (35) 

'fm..=  4^^ySt„(tjp.L)expj^^^^-?^ln^ff,,Lj=  1.51  x  10-\ 

where  we  have  used  the  fact  that  L  =  3.66  x  10^  cm,  L/d^-12  and  (see  Table  1) 
- 0.469  for  particles  which  are  small  compared  to  the  prevailing  gas-free-path 
(0.065  ^m).  In  this  case,  since  small,  the  Taylor  series  results  of  section  4.3  will  be 
useful,  i.e.; 

^Sl+(1  +  Q„)i„„,  (36) 

^p.L 


so  that  if  coagulation  is  negligible  (i.e.  (see  below))  then  Np.o/Np.L=  1  +  1-51 

X  I0~^  =  1.0015.  Thus,  under  such  conditions  total  number  densities  in  the  ambient 
atmosphere  would  only  be  about  0.15%  above  those  present  at  the  detector.  Based  on  the 
definition  of  the  dimensionless  coagulation  parameter  Cr„  it  is  also  possible  to  calculate  the 
level  of  particle  number  density.  A/p  ^  at  which  coagulation  effects  and  wall  loss  effects  on 
Np  „/A/p  L  would  be  comparable — i.e.  conditions  corresponding  to  Cf„  =  1.  In  the  present 
case; 


(-2a>.L)(iX,„)Np., 

4/d,L/SUCp.L)fo.L' 


(37) 


Setting  C,„=l  and  evaluating  the  factors  (-2}°J.l)  and  from  their  respective 
definitions  (see  Table  2,  first  row)  we  find  that  Cr„=l  when  Np,L  S  0(1.1  x  10’  particles 
cm*^),  a  level  that  could  be  experienced  near  pollutant  sources.  However,  on  the  assump* 
tion  that  particle  concentrations  are  frequently  much  lower  than  this  we  proceed  to 
calculate  the  (small)  corrections  needed  to  arrive  at  the  ambient  PSD  parameters  a,.,  and 
Vp.p.  These  can  be  obtained  from  equations  (30)  and  (32),  respectively.  In  this  case;  C‘®’  =  1, 
CV'=  1,  C'i^>  =  exp[-2(-0.469)ln^(2.46)]  =  2.138  so  that; 


and; 


<’i.O 

<^f.L 


(2.138- 1)(1.51  X  10-’) 
2  In  (2.46) 


1.00095 


exp(-i(2.138  -  I)(1.51  x  10"’))  =  0.9997, 

i.e.  when  dp  L  =  0.02^m  and  ap,L  =  2.46  the  systematic  corrections  are  modest  indeed. 
However,  the  present  formulation  has  the  merit  that  it  allows  rapid  calculations  of 
corrections  under  other  conditions  of  dp,L,  U, . . .  for  which  the  corresponding 

corrections  due  to  wail  loss  and/or  coagulation  may  be  appreciable,  including  cases 
requiring  the  graphs  of  sections  4.1  and  4.2  rather  than  the  small  {,„,,exi>ansions  of  section 
4.3  and  Appendix  A. 


5.2.  Laminar  flow  in  sampling  tube  upstream  of  an  airborne  condensation  nucleus  counter 

As  our  second  specific  example  consider  the  following  set  of  conditions  typical  of  an 
airborne  condensation  nucleus  particle  counter  (CNC)  system  for  flight  at  an  altitude  of 
25,000  ft.  Suppose  the  cabin  instrument  package  sees  an  aerosol  with  a  mean  particle 
diameter  0.014 /im  and  a  spread  parameter  of  2.46  in  air  at  a  volume  flow  rate  of 
24 cm’ s'*,  with  the  pressure  and  temperature  levels  established  at  0.4  atm  and  244  K, 
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respectively.  If  the  insulated  upstream  tube  that  conveys  the  aerosol  sample  to  the 
instrument  has  a  length  of  1 2 1  cm  and  an  internal  diameter  of  0.953  cm  then  we  seek  to 
estimate  the  correction  factors  IVp.o//Vp.L<  ‘^i.o/o'i.L  and  d,  ,/d,  ^  to  be  applied  to  data 
recorded  at  the  instrument  in  order  to  estimate  the  corresponding  aerosol  properties  in  the 
local  troposphere. 

If  we  tentatively  assumed  ‘fully-developed*  particle  mass  transfer  in  this  L/d,  =  128 
sampling  tube  system  (see,  however,  section  6.1  and  below)  we  would  proceed  as  follows; 
Under  the  stated  conditions  the  Reynolds  number  Ud,/v  can  be  calculated  to  be  only 
1.165  X  10^  so  the  gas  flow  will  be  laminar.  The  prevailing  mean-free-path  under  the  above- 
mentioned  conditions  is  about  0.2  Aim  so  that  Knp(iip  i,)=  14.4  and  most  particles  will 
experience  conditions  appropriate  to  free-molecule  (fm)  flow.  In  particular,  the  mean  size 
particle  will  be  characterized  by  a  Brownian  diffusivity  of  about  0.93  x  10’^  cm^s'‘, 
corresponding  to  a  Schmidt  number  Sc  =  v/Dp  of  about  2.96  x  10^.  For  fully-developed 
steady  laminar  flow  in  a  straight  circular  duct  it  is  well  known  that  Nu„  =  3.657  (Nusselt)  so 
that: 


St„(i',.) 


Nu„  3.657 

Re  Sc~(1.165x  10^)(2.96x  10^) 


1.06  X  10“*. 


It  would  follow  that; 


E 

Smax 


St. 


.(ftlexp  — ^^In^ff,  =0. 


852 X  10' 


(38) 


(39) 


since  (Table  1)6  =  —2/3.  If  the  coagulation  parameter  C,„«l  (see  below)  then  this  would 
correspond  to  A^p.o/^p.t  =  ^  085  due  to  particle  size-dependent  wall  loss  (based,  again,  on 
equation  (29)).  Noting  that  exp  { —  2(  — 2/3) In^ (2.46)}  =  2.947  we  find  (equation  (30)) 
that  the  remaining  correction  factors  are:  o/ffg.t  109  and  d,  ,/d,  L  *  0.973  (equation 
(32)).  Returning  to  the  possibility  of  coagulation  effects,  we  can  calculate  what  particle 
number  density  Np  l  would  be  required  to  cause  %  10" say.  Since  this  turns  out  to  be 
cfl  3x  10*  cm"^  in  the  present  case  (a  value  very  far  in  excess  of  ambient  atmosphere 
particle  number  densities  (0(10"^  cm”^)  at  such  altitudes)  the  neglect  of  coagulation  in 
such  sampling  tubes  is  seen  to  be  self-consistent. 

Returning  to  the  question  of  the  validity  of  our  tentative  assumption  of ‘fully-developed' 
particle  mass  transfer  (i.e.  A4,  negligible  ‘entrance  eflcct’)  we  note  that  while  the  viscous  flow 
itself  would  become  fully-developed  (r/d* -independent)  beyond  a  L/d»  of  only  about 
0.2  Re  =  23,  the  value  of  St„  will  not  approach  the  fully-developed  asymptote  until  Z./d»  is 
about  0.2ReSc  =  6900.  Since  the  actual  L/d».  of  the  system  is  128  we  see  a  posteriori  that  A4 
is  not  justified  under  such  laminar  flow  conditions  (see  section  6.1).  All  is  not  lost,  however, 
since  there  is  good  reason  to  believe  that  if  is  recalculated  using  the  actual  length- 
averaged  transfer  coefficient,  and  an  effective  value  of  h  appropriate  to  the  magnitude  of  the 
entrance  effect,  then  our  previous  formulas  (or  graphs)  will  provide  a  useful  first  approxima¬ 
tion.  The  ‘entrance-effect’  correction  factor  F,  to  the  fully-developed  value  of  Nu^  or  St,,  is 
known  to  depend  upon  the  dimensionless  distance  variable: 


1 

Re  Sc 


dw 


(40) 


For  fully-developed  laminar  flow  in  a  straight  circular  duct  a  convenient  approximation  to 
F(entrance)  is  (see,  e.g.  Rosner,  1986) 

F,  St  [1-K7.60C„)- *'"]*'*,  (41) 


which,  in  the  present  case,  gives  F,  =  3.29  corresponding  to  St„(t;L)  =  3.49  x  10" 
the  appropriate  “effective”  value  of  h  obtained  from 

:  [■•Sl-t-i”-'--']-: 


Using 


(42) 
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(see.  e.g.  Rosner  and  Tassopoulos,  1991)  we  then  recalculate  the  values  of  this  time 
obtaining  ^„„=0.30,  corresponding  to  Np  o/Np  ^s  1.3.  Recomputing  the  remaining  cor¬ 
rection  factors  now  gives  ffg.o/fft.L=  1^8  and  ^  *  0-95.  Based  on  the  present  work 

this  would  provide  the  best  current  estimate  of  the  required  correction  factors  in  this 
laminar  flow  aerosol  sampling  system.  More  rigorous  calculations  (now  underway)  will 
be  necessary  to  test  the  accuracy  of  this  plausible,  proposed  correction  procedure 
(for  relaxing  A4). 

In  the  light  of  these  two  simple  numerical  examples,  both  involving  highly  submicron, 
dense  {ca  Pp  =  1.83  gem ~^)  spherical  particles  in  straight  nearly  adiabatic  sampling  ducts 
but  one  with  laminar  flow  (section  S.2)  and  one  with  turbulent  flow  (section  S.l),  it  is  now 
appropriate  to  briefly  discuss  many  of  the  abovementioned  assumptions,  indicating,  where 
possible,  methods  to  generalize  the  treatment  to  include  less  ‘idealized'  cases  often  encoun¬ 
tered  in  sampling  from  industrially  important  devices  (e.g.  combustors  for  propulsion  or 
power  generation). 


6.  DISCUSSION  OF  ASSUMPTIONS  AND  GENERALIZATIONS 

It  is  prudent  to  consider  the  domain  of  validity  of  the  underlying  assumptions  of  section  2 
and  useful  to  indicate  here  how  readily  some  of  these  assumptions  can  be  relaxed,  thereby 
opening  the  door  to  the  quantitative  treatment  of  much  more  complex  sampling  duct 
situations. 


6.1.  Entry  effects 

Even  apart  from  development  of  the  viscous  gas  flow  itself,  the  mass  transfer  boundary 
layers  within  the  sampling  tube  require  a  ‘development  length’  before  St„  (equation  (6)) 
becomes  essentially  independent  of  zjd^.  For  laminar  mass  transfer  boundary  layers  the 
required  number  of  duct  diameters  is  approximately: 

=0.2(ReSc),  (43) 

which,  for  the  example  of  section  5.2,  is  about  6900,  far  in  excess  of  the  actual  Ljd^  of  the 
sampling  system  (128)  or  the  L/d»  required  for  development  of  the  host  gas  flow  itself 
(ffl  23).  Thus,  unless  much  smaller  particles  are  of  interest,  Assumption  4  will  not  be 
defensible  for  many  laminar  flow  sampling  systems.  A  rational  ‘correction’  for  this  can, 
however,  be  made  when  coagulation  is  negligible  (using  the  actual  length-averaged  St„- 
value  in  the  definition  of  s  and  modifying  the  effective  exponent  b,  as  discussed  in  section 
5.2)  or  specific  results  can  be  generated  for  laminar  sampling  systems  by  accurately 
including  such  mass  transfer  ‘entrance  effects’  (Rosner  and  Tassopwulos,  1991a).  For  turbu¬ 
lent,  high  Sc  systems  the  situation  is  much  improved  and  (Lld^)„^  will  ordinarily  be  less 
than  about  30,  a  condition  which  is  frequently  met*,  and  certainly  met  in  the  ground-based 
air  monitoring  application  discussed  in  section  5.1. 

6.2.  Variable  =  4A(z)/P(z) 

If  the  effective  duct  diameter  d,„  is  not  constant,  but  sufficiently  slowly  varying,  then  the 
quasi-one-dimensional  equations  (equation  (16))  remain  valid,  and  such  geometric  effects 
can  clearly  be  incorporated  in  the  numerical  integrations.  In  many  situations  d,ft  will  be 
piecewise-constant,  in  which  case  our  analysis  would  apply  to  each  segment  of  the  system 
but  with  the  eflcctive  lengths  of  the  segments  modified  to  allow  for  the  inevitable  ‘entrance- 
effects'  near  each  juncture  (see.  e.g.  section  6.1). 


•  Colkcl  I’l  III  (1977)  dcsiribc  a  diabutic  system  with  L/d,  %  2200  for  sampling  from  gas  turbine  combustors. 
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6.3.  Bends 

Because  of  the  secondary  flows  and  turbulence  modulation  associated  with  bends  in  a 
“piecewise  straight”  sampling  system  St„(z/d«,  Re,  Sc)  will  be  modified  in  accord  with  the 
curvature  d^/R  and  included  angle  of  the  bend.  Recent  (vapor)  mass  transfer  data  of  this 
type  (see,  e.g.  Sparrow  and  Chrysler  (1986)  and  Ohadi  and  Sparrow  (1989,  1990))  are  now 
being  mobilized  and  recast  so  that  it  may  be  conveniently  used  in  the  present  type  of  quasi- 
one  dimensional,  Sc>  1  analysis  (Rosner  et  al,  1991). 


6.4.  Heated  (cooled)  ducts 

Inadvertent  or  deliberate  heat  transfer  can  dramatically  alter  the  rates  of  fine  particle  loss 
to  the  wall,  and  the  particle  size-dependence  of  this  rate  (cf.  the  St„  exponent  b)  as  a  result  of 
particle  thermophoresis.  The  laws  of  particle  thermophoresis  (drift  down  the  transverse 
temperature  gradient)  have  been  under  extensive  development  in  the  last  IS  years  (see,  e.g. 
the  recent  review  of  Rosner  et  ai,  1990)  and  it  has  been  shown  that  for  submicron  particles 
deliberately  heated  ducts  can  be  used  to  sharply  reduce  wall  deposition  rates  [see, 
e.g.  Gokoglu  and  Rosner  (1986)  (before  the  gas  temperature  ‘catches  up'  with  the  wail 
temperature].  While  beyond  the  scope  of  the  present  paper,  it  can  be  shown  that  in  some 
cases  (e.g.  constant  wall  heat  flux  (addition))  the  present  results  can  be  used  as  a  first 
approximation  with  a  suitably  reduced  St„-value  and  a  modified  (more  negative)  b^rr'^alue 
(Rosner  and  Tassopoulos,  1991b).  In  some  cases  wall  cooling  may  be  necessary  near  the  inlet 
to  prevent  local  probe  failure  (as  in  extracting  samples  from  high  temperature  combustors 
(Colket  et  al.  (1977),  arcjets,  etc.).  This  can  cause  appreciable  local  particle  losses  to  the  wall 
unless  the  wall  cooling  is  accomplished  by  foreign  gas  transpiration  (see,  e.g.  Gokoglu  and 
Rosner,  1985).  Sometimes  downstream  heating  is  also  used  to  prevent  the  condensation  of  a 
co-present  vapor,  as  in  sampling  the  combustion  products  of  a  gas  turbine  combustor. 


6.5.  Single  mode  log-normality 

Aerosols  formed  from  several  distinct  mechanisms  may  exhibit  a  multi-modal  PSD,  each 
mode  of  which  can  be  represented  by  a  log-normal  distribution.  If  (as  is  often  the  case) 
further  coagulation  can  be  neglected  wi*hin  the  aerosol  sampling  tube  then  our  present 
results  can  be  applied  without  modification  to  each  mode.  In  the  presence  of  appreciable 
coagulation  this  is  no  longer  possible  and  results  would  have  to  be  obtained  using  an 
extended  coagulation  rate  theory  along  the  lines  of  Megaridis  and  Dobbins  (1990a)  to  deal 
with  particular  cases  of  interest. 


6.6.  Kiiff-interpolation 

If  an  appreciable  portion  of  the  aerosol  is  in  the  transition  regime  of  Knp  =  0(1),  then 
neither  of  the  above-mentioned  asymptotic  regime  results  (Knp  P  1  called  free-molecule,  or 
Knp  ^  1  called  continuum)  would  provide  an  accurate  description.  This  influences  both  the 
wall  loss  behavior  (via  its  effect  on  the  Brownian  diffusion  law  D(v))  and  the  coagulation 
rate  behavior  (via  the  collision  rate  constant  which  will  differ  from  either  of  the  previously 
considered  asymptotes  and  fi^).  A  general  treatment  of  this  regime  is  possible,  which  will 
introduce  an  additional  parameter  like  the  Knudsen  number  based  on  the  median  diameter 
(6t>p.i,/]i)’'^.  In  the  absence  of  coagulation  the  principal  effects  are  to  modify  St„(t\)  and  the 
effective  values  of  the  St„-exponent  b.  In  the  presence  of  coagulation,  use  can  be  made  of  the 
simple  and  successful  so-called  'harmonic-average'  approximation  for  the  coagulation  rate 
constant,  i.e.  fi  which  can  be  introduced  to  generalize  the  moment-method 

treatment  of  section  3.4  (Rosner  and  Tassopoulos,  1991b)  (for  an  alternative  but  closely 
related  approach,  sec  Pratsinis  and  Kim  (1989)  and  Biswas  et  al.  (1989)). 
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6.7.  Eddy-impaction 

When  the  sampling  duct  flow  is  turbulent  (Re  >  2.3  x  10^,  say)  and  a  significant  portion 
of  the  aerosol  is  associated  with  dimensionless  particle  stopping  times  t*  (=  ultjv)  greater 
than  about  16.5  (Sc)‘°  then  time-averaged  wall  losses  will  be  enhanced  by  the  mech¬ 
anism  of ‘eddy-impaction'  (see,  e.g.  the  brief  summary  of  available  correlations  contained  in 
Rosner  and  Tassopoulos  (1989)),  which  also  modifies  the  effective  value  of  the  exponent  b 
for  that  portion  of  the  aerosol.  Only  when  eddy-impaction  is  the  dominant  mechanism  and 
St„  -  (fp  )^  for  most  of  the  prevailing  particles  can  our  present  formulation  be  retained,  but 
with  6  =-(-  1.333  (see  Rosner  and  Tassopoulos,  (1989)).  However,  for  the  illustrative 
problem  discussed  in  section  S.l  we  find  that  eddy-impaction  would  not  set  in  until  dp 
=  3.8  ;<m,  a  size  beyond  which  there  are  few  particles  indeed.  Thus,  the  previous  treatment  is 
self-consistent  in  this  respect.  Other  effects  associated  with  the  mass  of  the  suspended 
particles  (e.g.  sedimentation  and/or  inertial  impaction  in  bends)  are  likewise  expected  to  be 
negligible. 

6.8.  Turbulent  coagulation 

In  principle,  new  mechanisms  of  suspended  particle-particle  coagulation /coalescence 
occur  in  turbulent  flows  (Friedlander,  1977)  which  could  be  incorporated  within  our 
moment  method  calculation  of  coagulation  via  the  appropriate  fi(u,  v, . . .)  functions. 
However,  even  in  turbulent  flows  sufficiently  small  suspended  particles  (the  ones  likely  to  be 
numerous  enough)  will  still  coagulate  according  to  the  Brownian  laws  discussed  in  section 
2,  especially  in  initially  well-mixed  (homogeneous)  submicron  aerosol  situations  at  the  inlet. 
In  many  cases  even  such  small  particles  may  not  be  numerous  enough  to  warrant  the 
inclusion  of  coagulation  within  the  sampling  tube  (see  the  examples  of  sections  5.1  and  5.2). 


C.9.  Axial  dispersion 

For  radially  nonuniform  velocity  steady  flow  in  a  duct  the  effective  axial  diffusivity  is 
known  to  be  larger  than  the  true  Brownian  diffusivity  by  the  Taylor-contribution  (see, 
e  g.  Butt,  1980;  Denbigh  and  Turner,  1971;  Probstein,  1989) 

!Re  ■  Sc 

(laminar) 

5  f  j  Ud^  (turbulent) 


so  that  Assumption  3  (section  2)  is  equivalent  to  satisfying  the  inequality; 


(45) 


This  criterion  is  readily  satisfied  in  the  turbulent  flow  example  (section  5.1)  but  is  only 
marginally  satisfied  in  the  laminar  flow  example  (section  5.2).  In  both  examples  the  effective 
axial  diffusivity  in  completely  dominated  by  the  fluid-dynamic  (Taylor)  contribution.  These 
results  again  suggest  (cf.  section  6.1)  that  the  idealizations  exploited  in  the  present  work  are 
more  appropriate  to  turbulent  flow  sampling  systems  (section  5.1)  than  laminar  flow 
sampling  systems  (section  5.2),  and  that  many  laminar  systems  will  warrant  a  more  accurate 
treatment  than  that  indicated  in  section  5.2,  even  including  the  above-mentioned  'entrance- 
effect'  correction  (section  6.1). 


6.10.  Aggregate  deposition/ coagulation 

For  non-coalcsced  aggregates  comprised  of  primary  particles  of  volume  t'l  each  depos¬ 
iting  'particle'  will  have  the  total  volume  r  =  nr,,  where  n  is  some  integer  (see,  e.g.  Megaridis 
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and  Dobbins,  1990b).  Such  aggregated  particles  will  have  orientation-averaged  Brownian 
dilTusivities  which  will  depend  not  only  upon  v  but  also  on  the  particular  morphology 
(arrangement  of  the  primary  particles).  In  each  case  we  can  write: 

^_g)nSt^_  glnSt„  /glnP\ 

5lnp  5lnSc  \<5lnu /' 

Thus,  when  coagulation  is  negligible,  a  knowledge  of  St„(t;),  where  t;  =  nt;,  and  of  the 
appropriate  value  of  din  D/d  In  v  for  the  morphology  is  question  will  also  allow  useful 
predictions  to  be  made  for  aggregated  particles  (see,  e.g.  Rosner,  1991;  Rosner  el  al.,  1 99 1 ).  A 
treatment  of  suspended  aggregate-aggregate  coagulation  analogous  to  that  outlined  in 
section  3.4  is  apparently  not  yet  available. 


6.1 1.  Non-power  law  capture 

If  the  Schmidt  number  dependence  of  St„  deviates  from  a  simple  power  law  and/or  D{v) 
departs  from  a  simple  power-law  then  the  value  of  the  exponent  b  will  not  be  a  constant 
(cf.  equation  (46))  and  corrections  to  the  convenient  computational  procedures  exploited 
here  will  be  necessary.  In  Rosner  (1989)  we  showed  that  these  corrections  depend  upon 
(d^lnSt„/dln^D)j:  and  the  local  spread  of  the  PSD.  This  provides  either  a  testable  criterion 
for  the  validity  of  power-law  behavior  in  any  particular  case  or  the  basis  of  a  systematic 
correction  procedure  for  the  calculation  of  the  factor  (cf.  equation  (12))  (see,  e.g.  Rosner 
and  Tassopoulous  (1991b)).  For  straight  duct  high  Schmidt  number  turbulent  convective- 
diffusion  transport  in  either  the  free-molecule  (fm)  or  continuum  (c)  limits  these  corrections 
vanish  and  the  exponents  (h-values)  presented  in  the  second  row  of  Table  1  are  true 
constants.  However,  it  should  be  realized  that  no  single  mechanism  of  particle  transport 
(hence  St„(t’)-law)  is  strictly  valid  over  the  entire  range  of  particle  sizes  (see,  e.g.  Rosner  and 
Tassopolous,  1989).  F'.\rthermore,  even  when  a  single  deposition  mechanism  dominates, 
departures  from  simple  power-law  behavior  may  be  associated  with  fluid-dynamic  non¬ 
idealities,  such  as  developing-(section  6.1)  and/or  secondary-flows  (cf.  section  6.3) 

6.12.  Vapor  growth/ particle  evaporation 

In  the  present  analysis  we  assumed  that  the  suspended  particles  are  stable  with  respect  to 
evaporation  or  growth  by  vapor  condensation  in  the  prevailing  carrier  gas  environment 
(A7).  In  “diabatic”  sampling  tube  cases  (often  encountered  in  extractive  sampling  from 
hostile  environments,  such  as  gas  turbine  engine  combustors)  this  may  not  be  valid  and  the 
analysis  of  section  3  would  then  have  to  be  generalized  by  including  the  contribution  to  the 
l/dMt/dz  due  to  vapor  growth  (see,  e.g.  Friedlander,  1977;  Megaridis,  1987;  Megaridis  and 
Dobbins,  1989)  coupled  with  balance  equations  for  the  vapor  concentration  and  the  bulk 
(mixing  cup-averaged)  gas  temperature.  The  present  class  of  moment  methods  would  also 
make  specific  calculations  of  this  type  feasible,  however,  the  assumption  of  log-rormality 
could  prove  overly  restrictive  in  many  such  instances  (see,  e.g.  Landgrebe  and  Pratsinis, 
1989;  Landgrebe,  1989;  Biswas  et  al.,  1987;  Mackowski  et  al.,  1991). 


7.  CONCLUSIONS,  RECOMMENDATIONS  AND  FUTURE  WORK 

A  simple,  rational,  moment-based  method  is  presented  for  correcting  aerosol  sampling 
data  for  the  systematic  upstream  effects  of  particle  size-dependent  wall  losses  and,  if 
necessary,  particle-particle  coagulation.  ‘Universal'  graphs  and  rational  formulae  are 
developed  and  provided  for  the  commonly  encountered  cases  of  log-normally  distributed 
aerosols  in  straight  adiabatic  ducts  of  constant  cross-section/‘wetted'  perimeter  for  the 
limiting  cases  of  free-molecule  particle  transport  (Knp>  1)  and  continuum  particle  trans¬ 
port  (Knp<|  1).  In  many  cases  the  leading  terms  of  a  Taylor  series  expansion  in  the 
dimensionless  distance  (proportional  to  (L/d„)St„(v))  will  be  adequate  to  provide  the 
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desired  correction  factors;  Np.„/Np,me„.  and  as  well  as 

r,  „/ig  n,,.,  and  In  "’ore  general  cases  backward  numerical  integration  of  the 

relevant  coupled  moment  equations  is  recommended,  allowing  several  of  the  idealizations 
discussed  in  sections  2  and  6  to  be  simultaneously  relaxed.  In  addition  to  providing  a  simple 
yet  rational  basis  for  correcting  aerosol  sampling  instrument  data  (i.e.  a  direct  solution  to 
the  canonical  ‘inverse'  problem  of  aerosol  sampling  theory),  the  present  results  can  also  be 
used  to  provide  preliminary  design  information  for  a  prospective  aerosol  sampling  system 
to  circumvent  the  future  need  for  uncomfortably  large  correction  factors  (predicted  or 
laboriously  measured)  in  the  ambient  environment  of  interest. 

The  present  approach  opens  the  door  to  several  generalizations  which  will  be  of  practical 
interest  (see,  also,  Rosner  and  Tassopoulos,  1990, 1991b) — especially  the  inclusion  of  addi¬ 
tional  deposition  and/or  coagulation  mechanisms,  effects  of  sampling  duct  bends,  variable 
cross-sectional  area,  diabatic  walls,  departures  from  single-mode  log-normality,  and  the 
presence  of  suspended  particulate  aggregates.  Several  of  these  extensions,  beyond  the  scope 
of  this  introductory  presentation  but  of  importance  in  sampling  from  combustion  turbines, 
rockets,  and  stationary  fossil  energy  power  plants  (necessarily  only  briefly  discussed  in 
section  6),  will  be  the  subject  of  future  communications  from  this  laboratory. 
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APPENDIX  A:  C'/'  COEFFICIENTS 

In  this  Appendix  we  provide  the  second-order  terms,  of  the  smal)-{  expansions  (see  equation  (27)).  The  wall- 
loss  (WL)  terms,  clearly  the  same  for  free-molecule  and  continuum  coagulation,  are; 

where  Nj,,  and  are  the  first-order  expansion  coefficients*  defined  in  equations  (28).  Next  we  provide  the 
coefficients  associated  with  coagulation’.  For  free-molecule  coagulation  we  find: 

0;.ca)=o 

M”  m  r  ) 

M,.ca)  =  C  ^ - ^ ^ -  >. 

The  corresponding  coefficients  for  continuum  coagulation  are: 

Tr-  ^  2[l-t-exp(iln’(T,  ,.)];v;  +  iexp(ln^e,,L)(-2^p-t-M'j  +  N;) 

l+exp(ilnXJ - 

i;.c(L}=o 


2  2^1  +expQln'e,.Ljj-t-^exp(ln'(;,  L)(-20,-hMi-(-N;) 
expdn'o,  ,.)-[l  -t-expdn'ej  L)] 


These  coefficients  can  prove  useful  to  evaluate  correction  factors  when  exceeds,  say  O.Z  but  does  not  exceed 
about  0.6.  For  values  of  {...much  above  0.6  direct  use  should  be  made  of  the  figures  of  section  4.2  or  specific 
numerical  integrations  in  accord  with  section  4.3. 
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•  In  the  present  notation  C,*’ »  N,  C',"  s  l.  CY*  ■  W'j  l. 
t  Note  that  CV"  m  N; ,*^(1.)  N;,c(U  etc. 


